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OCCURRENCE OF URINARY ACID MUCOPOLYSACCHARIDES IN 
THE HURLER SYNDROM E* 


By ALBERT DorRFMAN AND ANDREW E. Lorinczt 


DEPARTMENT OF PEDIATRICS, UNIVERSITY OF CHICAGO, AND LA RABIDA JACKSON PARK SANITARIUM, 
CHICAGO, ILLINOIS 


Communicated by Lowell T. Coggeshall, April 3, 1957 


The presence of abnormal amounts of a chemically unidentified substance in 
tissues of patients with the Hurler syndrome (lipochondrodystrophy gargoylism 
dysostosis multiplex) has been recognized.! Brante? recovered a substance with 
characteristics of an acid mucopolysaccharide from tissues obtained at autopsy 
(liver, meninges) of patients with this disorder. He suggested that this material, 
which was metachromatic and contained 3.9 per cent sulfur, 27 per cent hexosamine, 
and 26 per cent glucuronic acid, was similar to or identical with chondroitinsulfurie 
acid. Dawson* and Uzman‘ have also suggested that the “‘storage material” in 
affected tissues is a complex polysaccharide. Since acid mucopolysaccharides are 
not known to be synthesized in the parenchymal cells affected in this disease (e.g., 
liver, spleen, nervous system) it seemed possible that widespread deposition in 
tissues results from an increase in circulating mucopolysaccharides. Accordingly, 
an attempt was made to determine whether urine of afflicted individuals contains 
unusual amounts or types of mucopolysaccharides. This communication describes 
the isolation and identification of two acid mucopolysaccharides from the urine 
or a six-year-old Negro girl who has the typical grotesque skeletal and facial de- 
formities, hepatosplenomegaly, and corneal clouding characteristic of this dis- 
order. Preliminary studies on a second, male, patient with the Hurler syndrome 
indicated the presence of similar materials. 

Experimental.—After dialysis against cold running tap water, 48-hour urine 
samples containing thymol were concentrated in vacuo at 50° C. to '/i volume. 
Sodium hydroxide was added to make a 2 per cent solution, and, after 24 hours at 
room temperature, dialysis against cold running tap water was repeated. Ten 
milligrams of one-time crystallized trypsin was added and the mixture dialyzed 
against 0.1 M, pH 8.0, phosphate buffer at 38° C., with external stirring for 4-5 
days. Protein was precipitated by the addition of '/; volume of 40 per cent tri- 
chloroacetic acid and was removed by centrifugation after standing at 4° C. for 1 
hour. After dialysis against distilled water, the solution was concentrated to ap- 
proximately '/;) volume, and the ‘‘crude mucopolysaccharide” fraction was pre- 
cipitated by the addition of 4 volumes of 95 per cent ethanol in the presence of 2 
per cent sodium acetate. The precipitate which formed after 48 hours at 4° C. 
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was collected by centrifugation, washed with absolute ethanol and ether, and 
dried in vacuo in a desiccator containing calcium chloride and paraffin. 

These ‘‘crude mucopolysaccharide” fractions were purified by slab electrophoresis 
as described by Schiller et ai.6 One hundred and fifty to 200 mg. of “crude muco- 
polysaccharide” per 48-hour sample was isolated from the urine of this patient, as 
compared with 20-50 mg. per 48-hour sample obtained from the urine of normal 
individuals of comparable age and sex. The “crude mucopolysaccharide’”’ ob- 
tained from normals, when subjected to slab electrophoresis, yielded insufficient 
material for further identification. Recently, Di Ferrante and Rich® isolated from 
23 liters of pooled normal human urine 140 mg. of an acid mucopolysaccharide which 
was indistinguishable from chondroitinsulfurie acid—A (CSA-A). 

Slab electrophoresis of the “crude mucopolysaccharide” fraction derived from 
the urine of the patient with the Hurler syndrome indicated two peaks. Most of 
the material was present in the rapidly migrating peak, which had a mobility com- 
parable to that observed for CSA-A and chondroitinsulfurie acid—B (CSA-B, 
8-heparin). Similarity of behavior to CSA-A and CSA-B was also noted by paper 
electrophoresis, with respect both to mobility and to metachromatic staining with 
toluidine blue. 

The analyses (Table 1) of the more rapidly migrating peak are consistent with an 


TABLE 1 
COMPOSITION OF MUCOPOLYSACCHARIDES ISOLATED FROM URINE* 


Uronic Ester 
Hexosaminet Nitrogent Acid § Sulfur! N-Acetyl# 


“Crude fraction,’’ normal 1.00 3.42 0.71 
“Crude fraction,’’ Hurler syndrome 1.00 2.30 0.76 
‘Purified fraction,’’ rapidly migrat- 

ing peak, Hurler syndrome 1.00 1.32 0.55 1.37 
“Purified fraction,’’ slowly migrat- 

ing peak, Hurler syndrome 1.00 1.46 1.40 : 0.87 

* All analyses on each fraction were performed on the same solution, with the exception of the N-acetyl deter- 
mination on the rapidly migrating peak. Analyses are expressed as molar ratios with respect to hexosamine. 

+ Hexosamine was determined by the method of Boas (J. Biol. Chem., 204, 553, 1953), omitting the Dowex 


treatment. 
t Nitrogen was determined by a micro-Kjeldahl procedure. 
§ Uronic acid was determined by the method of Dische (J. Biol. Chem., 167, 189, 1947) 
! Ester sulfur was determined by the method of L. Anderson (Acta Chem. Scand., 7, 689, 1953). 
# N-acetyl determination by chromic oxidation according to the procedure of E. Wiesenberger (Mikrochim. 


Acta, 33, 51, 1948). 


acid mucopolysaccharide containing one sulfate group per disaccharide repeating 
unit. The values for hexosamine by the method used’ are generally about 10 per 
cent low. Some contamination by protein is present, as indicated by the slightly 
elevated nitrogen value and by the demonstration of free amino acids in hydroly- 
zates subjected to paper chromatography. The chromatograms were developed 
with tertiary amy! alcohol-formic acid-water (8:2:3), and the amino acids were 
demonstrated by the ninhydrin color reaction. 

The low uronic acid value obtained by the carbazole method® is characteristic of 
CSA-B. As has been found for CSA-B,*:? this material is resistant to the action 
of testicular hyaluronidase, as demonstrated by the method of Mathews, Roseman, 
and Dorfman. The optical rotation was [a]? = —68°, calculated on the basis 
of hexosamine content of the solution. Galactosamine was demonstrated by the 
method of Stoffyn and Jeanloz.!' 
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CSA-B has been isolated previously from skin,*: ' gastric mucosa,'* and beef 
lung.'* It differs from CSA-A in that it contains L-iduronic acid: '* instead of p- 
glucuronic acid. The presence of L-iduronic acid in this material has been demon- 
strated by paper chromatography following mild acid hydrolysis." The absence of 
sialic acid was indicated by a negative direct Ehrlich reaction." Preliminary 
studies'* indicate that this material has antithrombic activity similar to that of 
CSA-B'® obtained from beef lung. 

The more slowly migrating peak had an electrophoretic mobility somewhat less 
than CSA-A and CSA-B but greater than that of hyaluronic acid. Analyses of this 
material (Table 1) were consistent with an acid mucopolysaccharide containing one 
sulfate group per disaccharide unit. The high uronic acid value obtained by the 
carbazole method is of interest in view of the known high values which are ob- 
tained with heparin.’ This material was also resistant to the action of testicular 
hyaluronidase. Glucosamine was demonstrated by the method of Stoffyn and 
Jeanloz.'!_ The optical rotation was [a] 7° = +39°, based on hexosamine content. 
These data suggest that this substance is similar to that isolated by Jorpes and 
Gardell” from ox liver, which they named “heparin monosulfurie acid.” Meyer?! 
has recently indicated the existence of a similar substance in aorta and amyloid. 

Discussion.—The data presented indicate that a marked quantitative deviation 
in the excretion of acid mucopolysaccharides occurs in the Hurler syndrome. Al- 
though the two mucopolysaccharides identified occur in normal connective tissues, 
they have not been previously isolated from urine. Available evidence indicates 
that normal urine contains only small amounts of CSA-A, in contrast to the finding 
in the Hurler syndrome of relatively large amounts of CSA-B (6-heparin) and small 
amounts of a second mucopolysaccharide, tentatively identified as heparin mono- 
sulfuric acid. 

The site of origin of these polysaccharides is unknown, but it is unlikely that they 
are formed in all of the parenchymal cells in which abnormal accumulation occurs. 
More reasonable is the assumption that the fundamental defect resides in the con- 
nective tissues. The nature of this defect cannot be ascertained in the absence of 
more adequate knowledge of the metabolism of these mucopolysaccharides. In 
any case, it is probable that the distortion of tissues in this heritable disorder is re- 
lated to an abnormality of mucopolysaccharide metabolism. The Hurler syn- 
drome thus may be considered the first known example of a human disorder involv- 
ing such a metabolic anomaly. 


The authors are indebted to Julio Ludowieg, Louise K. Roth, and Minoru Mayeda 
for technical assistance. 


*This work was supported by grants from the Chicago Heart Association, the Helen Hay 
Whitney Foundation, and the National Heart Institute of the United States Public Health 
Service (H-311). 

+ Postdoctoral Fellow of the Arthritis and Rheumatism Foundation. 
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ASYMMETRIC DISTRIBUTION OF C% IN SUGARS FORMED 
DURING PHOTOSYNTHESIS 


By Martin Gipss* anp Orro KANDLERT 
DEPARTMENT OF BIOLOGY, BROOKHAVEN NATIONAL LABORATORY, UPTON, NEW YORK 


Communicated by D. D. van Slyke, March 12, 1957 


The concept of the conversion of CO, to labeled sugars during photosynthesis as 
reported by Calvin and coworkers! ? is the following: (1) carbon dioxide fixation 
involving the carboxylation of a symmetrically labeled two-carbon piece derived 
from ribulose diphosphate to yield phosphoglyceric acid predominantly labeled in 
the carboxyl carbon (CH,*OPO;H.-C*HOH-C**OOH), (2) a cyclic transketolase- 
transaldolase sequence involving fructose-6-phosphate, sedoheptulose-7-phosphate, 
and pentulose-5-phosphate to produce the symmetrically labeled CO, acceptor, and 
(3) conversion of the phosphoglyceric acid via the Embden-Meyerhof sequence to 
yield a hexose sugar predominantly labeled in carbon atoms 3, 4 (C-3, C-4) with 
carbon atoms 1, 6 equal to 2, 5 (C*-C*-C**-C**-C*-C*).*-* In earlier papers the 
hexose degradation data were obtained by the Lactobacillus casei” procedure, which 
yields pairs of carbon atoms. In the present report, we have investigated the dis- 
tribution of tracer in sugars produced during photosynthesis in CQO, using the 
Leuconostoc mesenteroides* * degradation method, which permits a determination 
of the C'* content of the individual carbon atoms. The results of these experi- 
ments and their implications for the pathway of carbon dioxide to carbohydrate 
during photosynthesis are reported in this communication. 

Methods and Materials.—The Chlorella was grown as described by Kandler," 
except that the light source was fluorescent lamps with an intensity of approxi- 
mately 1,000 foot-candles. The organisms were harvested by centrifugation, 
washed twice with water, and suspended in distilled water or potassium phosphate 
buffer. The tobacco, sunflower, and Canna leaves, including petioles, were re- 
moved from mature greenhouse stock plants. 

The 10- and 60-second experiments with algae were carried out in a “lollipop” 
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apparatus, as described by Calvin.!_ Conventional Warburg manometry techniques 
were employed for the longer-time algal experiments involving the effect of glucose 
on photosynthesis. In the leaf experiments, the photosynthesis apparatus was a 
desiccator type of chamber. In all experiments the light source was Sylvania 
Birdseye Floodlites (150 watts, 120 volts). 

In general, the photosynthesis was terminated by immersing the test substances 
in boiling 80 per cent ethanol. The 80 per cent ethanol extract and two 20 per cent 
ethanol washings were concentrated under vacuum to approximately 1 ml. The 
residue was analyzed by two-dimensional paper chromatography (solvents: water- 
saturated phenol and butanol—propionic acid; paper: oxalic acid-washed What- 
man No. 4). The area of the chromatogram occupied by the sucrose was eluted, 
and the eluate was hydrolyzed by means of invertase and rechromatographed using 
water-saturated phenol as solvent to separate the monosaccharides. The residue 
remaining after the 80 per cent ethanol extraction was treated twice with cold 20 
per cent ethanol, then washed several times with hot water. 

Following the washing, the starch glucose was derived by hydrolyzing the cells 
with 1 N HCl for 45 minutes. Paper chromatography indicated that essentially all 
the tracer resided in glucose. The hydrolyzate was brought to pH 6-7 with 1 NV 
KOH. 

Glucose was degraded by fermentation with Leuconostoc mesenteroides, which has 
been shown by Gunsalus and Gibbs to produce CO, from carbon atom 1, ethanol 
from carbon atoms 2 and 3, and lactate from carbon atoms 4, 5,and6. The ethanol 
was oxidized to acetic acid by heating for 2 hours at 90° C. with 0.5 gm. of potas- 
sium dichromate in 4 N sulfurie acid. The activity of the carboxy! carbon of the 
acetic acid was obtained by the method of Phares. Methyl carbon activity was 
determined by difference between the activity of the whole acetic acid molecule 
obtained by a persulfate oxidation and the carboxyl carbon determination. Lac- 
tate was subjected to an oxidation with chromium trioxide to yield CO, and acetic 
acid. In order to avoid the addition of carrier ethanol or lactate, about 400-500 
umoles of glucose were fermented. 

Radioactivity measurements were made with a Bernstein-Ballentein methane 
flow counter. All samples were converted to barium carbonate and counted at 
infinite thickness. Activity is expressed in Tables 1-3 as millimicrocuries per 
milligram of carbon (muC/mgC). 

Results and Discussion.—All the glucoses formed during the short periods of 
photosynthesis in C“O, (Table 1) possess an asymmetrical distribution of C™. 
The ratio of carbon atom 3 to carbon atom | (C-1, aldehyde carbon) and/or C-2 is 
less than that of C-4 to C-5 and/or C-6. This type of labeling pattern in the glu- 
cose moiety of sucrose and starch is similar to that reported in our earlier communi- 
cation!! for the glucose phosphate esters (monophosphate, uridine diphosphate 
glucose, and an unidentified glucose phosphate). 

The labeling pattern of the three glucose phosphate esters of our previous com- 
munication can be compared with the data tabulated in the first line of Table 1, 
since the four glucose units were derived from the same experiment. An average 
of the four glucose degradations in which the activity of the various carbon atoms 
are expressed in percentage of the total glucose molecule is as follows: C-1(7.9). 
C-2(7.6), C-3(33.8), C-4(42.5), C-5(3.5), C-6(4.7). 
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TABLE 1 


DIsTRIBUTION oF C4 IN GLUCOSE 
-——T RACER ConTENT oF GLUCOSE CaRBON ATOMS-—— 
Lieut INTENSITY GLUCOSE (MyC/meaC) 
PLANT (Foor-CaNnDLEs) Time Source 3 
Chlorella* 4,000 10 sec. Starch : : ; .90 
Chlorellat 4,060 60 sec. 
Chlorellat 700 45 min. q ‘ 2 
Tobacco$ 4,000 50 sec. j : . 1 
Tobacco$ 100 180 sec. , : 3 
Sunflower § 70 15 min. ; f ; .2¢ : 
Canna 2,000 24hrs. Sucrose : j 5.08 5.08 5.12 


* Thirty milliliters of a Chlorella suspension (1.5 ml. packed cells per 100 ml. water) was illuminated 5 minutes in 
a nitrogen atmosphere before the introduction of 10 wmoles of NaHCO; containing 135 microcuries. 

+ Thirty milliliters of a Chlorella suspension (1.5 ml. packed cells per 100 ml. water) was incubated with 10 
numoles of NaHCO;, containing 135 microcuries, in the dark for 5 minutes before photosynthesis occurred. 

t Five milliliters of Chlorella suspension (1.0 ml. packed cells per 100 ml. M/30 potassium phosphate buffer, pH 
5.6) was illuminated in a Warburg vessel. Introduction of NaHCO; (50 uwmoles, 40 microcuries) at same time light 
was applied. 

§ The leaves were illuminated 5 minutes in an air atmosphere before the introduction of 10 umoles of C'O2 con- 
taining approximately 100 microcuries. 


00 
4 
6 
5 
29 


| 
12 
8 


TABLE 2* 


Errect or GLUCOSE FEEDING ON PosrtT10N LABELING IN GLUCOSE 
oF SucRoOsE FORMED DURING PHOTOSYNTHESIS IN CQO, 
TRACER Content or Giucose Carson AtomMs—————— 

TIME (mpC/meC) 
(Mrn- CO. = ———————Witnovut Giucose————__—._ ——__—_-WiTn Giucose————— 
UTES) (MOLES) 1 2 3 4 5 6 1 2 3 4 5 

15 10 3.16 2.92 4.00 5.21 1.25 1.67 8.08 8.34 13.3 13.8 2.46 

45 50 2.49 2.16 2.52 3.52 1.12 1.46 8.42 8.33 14.9 18.3 3.58 
180 200 4.34 4.31 4.25 5.12 3.98 4.36 3.39 2.53 4.92 6.23 1.30 


* The Chlorella were starved by shaking the cells in water at room temperature for 16 hours in the dark. 
centrifuging, they were suspended in M/30 potassium phosphate buffer, pH 5.6 (1.0 ml. of packed cells per 100 ml. 
buffer). The experiment was carried out in six 150-ml. Warburg vessels, one set of three containing 5 ml. of the cell 
suspension, the other, 5 ml. of suspension containing 55 w~moles/ml, glucose. The side arms carried either 10, 50, 
or 200 umoles of NaHCO; containing 2.75 microcuries. After the vessels were shaken 10 minutes in the dark, 
the NaHCO; was tipped in and illumination was begun. An illumination of 700 foot-candles was used, which 
resulted in glucose respiration (5-fold over endogenous) and photosynthesis being at the compensation point. After 
the appropriate time, the contents of the vessels were pipetted into cold centrifuge tubes containing ice water. 
After the cells were centrifuged and the supernatant fluid discarded, they were extracted twice with hot water. 
The glucose moiety of the sucrose was isolated for degradation from the hot-water washings, as described under 
BR sae and Materials.’’ The starch glucose was obtained for degradation as described under ‘‘Methods and 
Materials.” 


TABLE 3* 


Errect oF GLUCOSE FEEDING ON PosITION LABELING IN STARCH 
FORMED DURING PHOTOSYNTHESIS IN C!4Q, 


——-———_ TRACER ConTent or GLucosE CARBON ATOMS 
TIME MuC/meC) 
(M1n- CO: —_——————-Without Glucose——-——__.~ _——~——-———-With Glucose-————_----_~ 
UTES) (uMOLES) 1 2 3 4 5 6 2 3 4 5 6 
15 10 1.03 0.71 1.71 2.08 0.37 0.53 1.61 1.59 2.88 3.08 0.60 0.82 
45 50 ve ati as ze by: ies 3.00 2.93 5.88 7.15 1.52 1.57 
180 200 2.21 2.32 2.42 2.83 2.31 2.11 4.45 3.60 8.42 10.42 2.44 2.64 


* See note to Table 2. 


Our earlier communication also reported degradation data for two hexose phos- 
phate esters corresponding to the glucose data listed in the second line of Table 1. 
The average values in percentage of the total glucose molecule for the three deg- 
radations are C-1(5.5), C-2(5.6), C-3(41.0), C-4(46.0), C-5(0.8), C-6(1.2). 

The isotope-distribution data recorded in Table 1 were obtained with Chlorella 
and leaves permitted to photosynthesize in CO, for short time periods. Since the 
findings of Kandler’ and Simonis and Grube" have indicated that exogenous 
glucose may influence the photosynthetic process, an experiment was carried out 
in which the rate of glucose respiration and photosynthesis were at the compensa- 
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tion point, with the desire that the exogenous C!*-glucose would enter the photo- 
synthetic cycle, dilute the compound or compounds which are closely connected 
with CO, reduction, and thereby prolong the period of asymmetric distribution of 
label in the newly formed carbohydrate. The isotope-distribution data tabulated 
in Tables 2 and 3 strongly indicate that a substrate of respiration can, at least 
at the compensation point, enter the photosynthetic cycle and be converted into 
compounds involved in the pathway between carbon dioxide and the hexose 
molecule. Thus, even following a photosynthesis period of 180 minutes, the 
Chlorella which received an exogenous C!*-glucose supply still possessed sucrose 
and starch whose glucose moieties exhibited a marked asymmetry in isotope dis- 
tribution, whereas the control algae (minus glucose) contained a uniformly labeled 
hexose unit. 

At present, three mechanisms can be suggested to account for a glucose molecule 
possessing a labeling pattern in which the ratio of C-3 to C-1 and/or C-2 is less 
than that of C-4 to C-5 and/or C-6. 

One pathway, which was suggested to us by Dr. Paul Marks, we have titled the 
‘ exchange and dilution pathway.” This “pathway” involves a slight modifica- 
tion of the photosynthetic cycle published recently by Bassham, Barker, Calvin, 
and Quarck.'* It may be summarized in the following manner (the values as- 
signed to the carbon atoms of the various compounds are arbitrary and are ex- 
pressed as per cent of total activity; they are based on the data of Calvin and 
co-workers) : 


et) (2 C4 


| | | 
O=C* O=C* O=—C” 


ae + (150 ae ~ E100 (0 
; l 


dilution exchange 


(‘100 


| 
| 
C10 aie ; Cm , > (20 
1 2 


| | | hie | 
C 10 (20 C 20 (20 C2 


bs be» ( 


1. The conversion of the aymmetrically labeled ribulose diphosphate and carbon 
dioxide into two moles of phosphoglyceric acid, followed by a conversion to a 
mole of glyceraldehyde-3-phosphate and one of dihydroxyacetone phosphate 
catalyzed by the carboxylating enzyme, phosphoglyceric acid transphosphorylase, 
triose phosphate dehydrogenase, and triose phosphate isomerase. Carbon dioxide 
has been assigned the highest value, since the ribulose diphosphate has presumably 
been diluted by endogenous material. 

2. Dilution of the isotopic dihydroxyacetone phosphate by a pool of inactive 
dihydroxyacetone phosphate derived from the breakdown of endogenous fat. 
Dilution prior to condensation would produce an asymmetrically labeled fructose 
phosphate molecule. 

3. A rapid exchange catalyzed by transketolase between the “active glycolalde- 
hyde” moieties of fructose-6-phosphate and pentulose phosphate. The result of 
the exchange is an increase in tracer C-1 and C-2 of the hexose phosphate molecule. 
Thus a combination of dilution by endogenous dihydroxyacetone phosphate and 
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an exchange catalyzed by transketolase could produce a hexose unit whose isotope 
distribution is in accord with the obtained data. The weakest point in this hy- 
pothesis is the source of dihydroxyacetone phosphate in the starved algal cells. An 
investigation involving a kinetic study on the specific activity of the various inter- 
mediates would shed light on this point. 

A second pathway for forming an asymmetrically labeled hexose molecule is 
based principally on our observation that C-4 has the highest activity of the hexose 
carbon atoms. Since the C!O, has the highest specific activity of all components 
of the photosynthetic reactions, a direct conversion of the CO, to C-4 of the hexose 
unit would keep this carbon atom of the sugar molecule with the highest specific 
activity until a uniformly labeled molecule was synthesized. 

The hypothesis is proposed that pentulose diphosphate is cleaved between C-3 
and C-4. The diose phosphate piece (C-4 and C-5) is condensed with a reduced 
form of CO, to yield glyceraldehyde-3-phosphate: 


(130 Cx» Cc* 
l | | 
0=C* O0=C” O=C* 
l | | 


(50 4+ (150 <> C50 =e (C50 


--}- + | 
cw OH C™ cu 
| | | 


C0 C0 C110 
G10 es 


The upper portion of the diphosphate yields dihydroxyacetone phosphate. A re- 
versal of the Embden-Meyerhof pathway would produce an asymmetrically labeled 
hexose unit. After this point the photosynthetic cycle would be similar to that 
described previously. !* 

The third mechanism is a compromise between the two proposals. It would pro- 
duce an asymmetrically labeled hexose molecule in the following manner: (1) 
addition of CO, to C-4 of the pentose, (2) a reduction of this primary addition com- 
pound by the reducing power generated by light, (3) a rearrangement in which the 
new carbon atom becomes C-4 of the hexose, and (4) a lack of reducing power re- 
sulting in the primary addition product being split into two molecules of PGA. 

It is evident that the basic difference between the three schemes is the fate of 
the carbon dioxide in the initial steps. The first and third proposes a carboxylation 
catalyzed by the ribulose diphosphate carboxylating enzyme, followed by a reduc- 
tion either as a whole unit or after splitting into two molecules of PGA, while the 
second speculates on a direct reduction of the CO, subsequent to an addition. Fur- 
ther studies will be required to elucidate the significance of this asymmetric type of 
label in the hexose molecule. 

Summary.—Photosynthesis was carried out with Chlorella and the leaves of 
higher plants using C'O. under various conditions. The glucose isolated from 
starch and sucrose produced during short periods of photosynthesis possessed an 
asymmetric distribution of tracer. This pattern of label was also obtained during 
longer periods of photosynthesis in the presence of a supply of exogenous glucose. 
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The distribution of tracer is not that which is expected if two equal triose phos- 
phates combined to yield a hexose phosphate. 


The authors are indebted to Dr. R. C. Fuller for discussion, helpful suggestions, 
and assistance in carrying out the short-time experiments. They are also grateful 
to Dr. Jerome Schiff and Dr. G. Kandler for laboratory assistance. 
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MULTIPLICATION OF ANIMAL CELLS IN SUSPENSION 
MEASURED BY COLONY COUNTS* 


Ry R. B. L. Gwarkin, J. E. Titi, G. F. Wuarrmore, L. Stmurnovircn, 
AND A, F. GRAHAM 


CONNAUGHT MEDICAL RESEARCH LABORATORIES AND DEPARTMENT OF MICROBIOLOGY 
SCHOOL OF HYGIENE, UNIVERSITY OF TORONTO, TORONTO, ONTARIO, CANADA 


Communicated by A. Lwoff, April 11, 1957 


During the past few years the development of new techniques for the cultivation 
of an:mal cells in vitro has facilitated the quantitative study of many aspects of 
cell biology. At present the most commonly used method of propagating cell 
strains is based on the ability of cells to multiply while attached to a glass surface. 
The cells may be subcultured by removing them from the surface into suspension 
and then distributing them into other vessels, where they will again adhere to the 
glass and populate the surface. This procedure has been developed by Earle and 
his associates! ? inte the so-called quantitative replicate culture technique and 
applied to a variety of studies with animal cells. Despite the technical advance 
represented by this method, there are, nevertheless, serious experimental limitations 
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inherent in the use of glass surfaces for cultivating large cell populations. Perhaps 
the most obvious of these is the problem of removing representative samples at 
will during the growth of a cell population. In addition, subculture requires the 
removal of cells from the surface, with consequent risk of cell injury. 

In order to obviate these restrictions, efforts have been made by several investi- 
gators to propagate animal cells in suspension.*~* Previous communications from 
this laboratory have described a technique that accomplishes this purpose.’ It 
was shown that, for several strains, the growth of animal cell populations in sus- 
pension, estimated by enumeration of cells in a counting chamber, is similar to 
that of bacteria, in that lag, logarithmic, and stationary phases were observed. 
Cells could be kept continuously in the logarithmic phase of multiplication by 
subculturing at frequent intervals. 

While such a system would seem to be ideally suited to biochemical studies, its 
full usefulness cannot be realized without some knowledge of the number of viable 
cells in the cultures. To obtain this information, we have used the clone technique 
recently developed by Puck and his associates.* * In the present paper it will be 
shown that cells grown in suspension readily form colonies on a glass surface and 
that, for two strains at least, they will plate with a high and constant efficiency 
during the logarithmic phase of growth. 

Two strains of cells were used in this work: (1) Earle’s L strain, a clone derived 
originally from a culture of adult mouse connective-tissue cells, and (2) a clone, 
selected in this laboratory, from an altered monkey kidney cell population provided 
by Dr. R. C. Parker.!! | The method for cultivating cells in suspension, and the 
medium, have been described in detail elsewhere.’ In brief, large test tubes con- 
taining the cell suspensions were rotated at 60 rpm around their long axes at 37° C. 
The medium contained 20 per cent horse serum in solution CMRL-1066,!! which is 
a defined mixture of salts, amino acids, vitamins, coenzymes, and other metabolites. 
Cells of both strains multiplied readily under these conditions, with generation times 
of 20-25 hours in the logarithmic phase of growth. Stock cultures were kept con- 
tinuously in the logarithmic phase by frequent subculture. 

When cells from a logarithmic-phase culture were placed in a stationary container, 
macroscopically visible colonies appeared on the glass surface after 6-10 days’ 
incubation. A linear relationship was obtained between the number of cells plated 
(10-500) and the number of colonies formed. The procedure for plating the cells 
was similar to that used by Puck, Marcus, and Cieciura® for suspensions of cells 
from trypsin-treated monolayers. For the L strain, daily enumeration of the cells 
in ten randomly selected colonies showed that multiplication of the individual cells 
began 12 hours after plating and continued logarithmically, with a doubling time 
of 18 hours, until the colonial populations reached 500 cells or more. Since we found 
in these experiments that both the plating efficiency, i.e., the percentage of plated 
cells that developed into colonies, and the rate of cell increase within colonies varied 
greatly with the particular lot of horse serum employed, a selected serum was used 
in the remainder of the work. 

In order to determine the increase in a cell population by the colony-count method 
a rotating culture containing about 5 X 10* cells per milliliter was prepared from a 
logarithmic-phase stock suspension. At intervals samples were removed and the 
cell densities estimated with a hemocytometer.'!? Serial dilutions of the sample 
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were made in solution 1066 at room temperature. One hundred cells were plated 
in each of six Petri dishes (6-cm. diameter) containing 5 ml. of complete medium 
(solution 1066 supplemented with horse serum) previously equilibrated at 37° C. 
in an atmosphere of 5 per cent CO,. After incubation of the plates at 37° C. for 
8-10 days in the 5 per cent CO, atmosphere, the colonies were washed with buffered 
saline, fixed with Bouin’s solution, treated with Giemsa’s stain, and scored. The 
colonies of L cells on a representative plate are shown in Figure 1. 


Fig. 1.—Colonies of L-strain cells obtained by plating 100 cells from a 
suspension in the logarithmic phase of multiplication. The colonies were 
fixed and treated with Giemsa’s stain. 

Figure 2 shows the growth curves obtained for L cells by hemocytometer count 
(upper curve) and by plate count (lower curve). Each point on the lower curve 
represents the average of six plates, or 400-600 colonies. The method of least 
squares was used to fit the exponential portion of this curve and indicated that the 
standard deviation of the points from the line was about 12 per cent. A similar 
set of curves is shown in Figure 2 for a suspension of the monkey kidney cells. 
The average plate counts deviate from the linear part of the curve by 8 per cent. 
The errors, which include day-to-day variation and sampling errors, are comparable 
to those obtained in colony counts with bacteria. 

It will be seen that the plating efficiencies in both experiments remained constant 
at about 75 per cent throughout the logarithmic phase of growth. During the 
stationary phase, however, the plating efficiencies may decrease significantly (Fig. 
3). This decrease was very marked in some other experiments, indicating that 
the cells may undergo an extensive physiological change during this period. 
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Fic. 2.-—Multiplication of L-strain cells in suspension. Upper curve: hemocytometer count; 
lower curve: colony count. 


It is of interest now to inquire whether the 75 per cent plating efficiency observed 
during logarithmic growth represents all the viable cells in the population or 
whether some inadequacy in the plating procedure or in the medium permits only 
75 per cent of the cells to form colonies. In an effort to answer this question, we 
have utilized the feeder-layer technique developed by Puck e¢ al.+® Table 1 
shows the results of several experiments in which monkey kidney cells from loga- 
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Fic. 3.—Multiplication of altered monkey kidney cells in suspension. Upper curve: hemocy- 
tometer count; lower curve: colony count. 


rithmic-phase suspensions were plated over feeder layers of giant cells obtained from 
X-irradiated cultures of the same cell strain. In each case the plating efficiency 
was higher in the presence of feeder cells than in their absence and was frequently 
in the neighborhood of 100 per cent. These observations indicate that, during the 
logarithmic phase of multiplication, nearly all the cells possess the ability to form 
colonies under the improved plating conditions provided by the presence of feeder 
cells. It should be mentioned, however, that for many practical purposes, e.g., 
studies on virus and radiation sensitivity of cells during the period of logarithmic 
increase, a plating efficiency of 75 per cent achieved without the use of feeder cells is 
adequate, 





BIOCHEMISTRY: GWATKIN ET AL. Proc. N. A. 8. 


TABLE 1 
PLATING EFFICIENCY OF ALTERED MONKEY KIDNEY CELLS 
Without With 
Feeder Cells Feeder Cells* 
Experiment (Per Cent) (Per Cent) 
1 58 98 
2 58 92 
3 62 102 
4 66 90 


*5 X 10¢ irradiated cells were added to each Petri dish to form the feeder layers. 


It is evident that the culture systems described here have the essential properties 
that will allow a highly versatile and precise experimental approach to many 
problems of animal cell physiology and biochemistry. They provide populations of 
individual cells that multiply in suspension, from which representative samples 
may be taken at will. During the logarithmic phase the number of cells that will 
form colonies on a glass surface without feeder cells is a constant and high propor- 
tion of the total population, and the amounts of RNA, DNA, and protein per cell 
are constant.’ In these respects there is a marked similarity between these popula- 
tions of animal cells and many bacterial populations. 

Summary.—When animal cells of two different strains were withdrawn from 
populations growing in suspension, they multiplied, without an appreciable lag, to 
form colonies. While the percentage of cells that formed colonies varied widely 
when they were plated from the stationary phase of growth, 75 per cent of the cells 
gave to colonies when plated at any time during the logarithmic period. Almost 
100 per cent of the cells in the logarithmic phase grew into colonies if they were 


plated over feeder layers of X-irradiated cells, showing that all, or nearly all, of the 
cells were capable of multiplying if provided with the proper nutritional environ- 
ment. The cell-culture systems described should facilitate a variety of studies on 
the physiology and biochemistry of animal cells. 


We wish to acknowledge our gratitude to Dr. R. C. Parker for a great deal of 
helpful discussion and advice during the course of this work. 
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Health Service of the National Institutes of Health of the United States, the National Cancer 
Institute of Canada, and the W. B. Boyd Memorial Fund. 
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12 Tn our previous work on the growth of animal cells in suspension (see nn. 6 and 7), cell densities 
were determined in a Petroff-Hauser bacterial counter rather than with a hemocytometer. Subse- 
quent comparison of the two types of counter showed that estimates made with the hemocytometer 
are correct, while those determined with the bacterial counter were high. Cell densities reported 
in previous publications may be corrected by multiplying them by the factor 0.65. 


ON THE DISTRIBUTION OF PURINE AND PYRIMIDINE BASES IN 
THE NUCLEIC ACID OF TOBACCO MOSAIC VIRUS 


By Roger G. Harrt* 
VIRUS LABORATORY, UNIVERSITY OF CALIFORNIA, BERKELEY, CALIFORNIA 


Communicated by W. M. Stanley, March 22, 1957 


Reddi and Knight! have recently investigated the products of ribonuclease di- 
gestion of the nucleic acids of several strains of tobacco mosaic virus. After the 
enzyme had been allowed to act on the isolated ribonucleic acids, the trichloroacetic 
acid—precipitable fraction (hereafter called ‘‘core’’) of each digest was analyzed to 
determine the molar ratios of purine and pyrimidine bases. As might be expected 
from the specificity of pancreatic ribonuclease,’ the core material was found to con- 
tain several times as many purines as pyrimidines, whereas the ratio of purines to 
pyrimidines is close to unity in the whole nucleic acids. And just as no significant 
differences have been found in the over-all composition of the nucleic acids of the 
strains of tobacco mosaic virus,’ no significant differences were observed in the base 
compositions of the cores. 

In the present study an attempt is made to determine from theoretical consider- 
ations whether the similar compositions of the cores are indicative of some common 
arrangement of purines and pyrimidines in the nucleic acids. In order to do this, a 
hypothetical nucleic acid is considered to undergo the treatments described by 
Reddi and Knight and to serve as a control in their experiments. The control ma- 
terial is assumed to be a single long-chain nucleotide polymer with the characteris- 
tic ribose-to-phosphate linkages of a ribonucleic acid and to have the same relative 
numbers of adenylic, cytidylic, guanylic, and uridylic acid residues as do the nucleic 
acids of the viral strains. However, the sequence of these residues along the hy- 
pothetical polymer is assumed to be completely random, and for this reason it 
serves as a control; if the core compositions for the viral-strain nucleic acids do not 
differ from the core composition which would be obtained for the random polymer, 
then they are not necessarily indicative of any specific nucleotide sequence common 
to the strains. The validity of this test can be appreciated if it is remembered that 
the random arrangement is not any particular sequence of nucleotide residues but is 
a mathematical composite resembling a majority of the vast number of possible 
sequences. If the nucleic acids of the viral strains have core compositions consist- 
ent with such a model, it can be concluded that they are typical by this criterion, 
but not that the strains are identical or even similar with respect to their particular 
nucleotide sequences. 

The question of how closely a viral nucleic acid resembles a random polymer has 
another, more general significance. Recent experiments‘ have indicated that a 
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molecule of viral ribonucleic acid having as few as 1,000 nucleotide residues can by 
itself bring about mosaic infection in plants and the production of complete nucleo- 
protein virus particles. Thus whatever additional information the cell may require 
in order to begin synthesizing the abnormal or viral materials must be somehow 
coded into a sequence of no more than 1,000 nucleotides. The number of different 
possible sequences of this length is of the order of 10®?, a number which dwarfs 
any conceivable figure for the number of different kinds of molecules actually in- 
volved in life-processes. But if there is any simple regularity in the nucleotide 
sequence, such as a repetition of one pattern involving only a few residues, then the 
possibility for information storage is severely limited. One indication of the 
presence of such a simple repeat in viral nucleic acid would be a core composition 
which differed significantly from that of a random sequence. 

The Random Nucleic Acid.—Consider a single, unbranched polynucleotide chain 
of A purines and B pyrimidines in a random sequence. The chain may be con- 
veniently thought of as a series of inter-pyrimidine segments containing 0, 1, or 
more purines, with each such segment separated from its neighbors (or “neighbor,” 
in case of an end segment) by single pyrimidines. If all A purines are arranged in 
the B + 1 segments in a random manner, the probability that a given purine will 
lie in any one of the segments is the same for all segments and is therefore equal to 
1/(B + 1). If A and B are both large, the variable ny, representing the number 
of segments having just N purines, will be distributed according to the Poisson law: 

Be-®R,® 


tin = 


aes where Ro = 4 (1) 

Ribonuclease Digestion and Precipitation of the Core.—Since pancreatic ribonu- 
clease appears to cleave the nucleic acid chain only at one side of each pyrimidine 
residue (hydrolyzing the phosphoester bond at the ribose 5’ position), the fragments 
after digestion will consist of the original purine segments (note that some are of 
zero length), each attached to a terminal pyrimidine.’ 

At this point it is necessary to make some assumption as to the selective action of 
6 per cent trichloroacetic acid (as used by Reddi and Knight) in precipitating the 
various fragments in the digest. For algebraic simplicity and for want of any 
empirical information on the solubilities of polynucleotides under these conditions, 
the following will be assumed: (1) that no fragments containing less than some 
critical number M of nucleotide residues (i.e., M — 1 purines) are precipitated; 
(2) that all fragments containing more than M residues are precipitated; and (3) 
that some fraction 0 < f < 1 of the fragments containing just M residues are pre- 
cipitated. According to these conditions, the fraction of the complete digest (in 
terms of nucleotide residues) which is precipitable is 


fnu.aM + > ny(N +1) 
oe ae en ; 





The molar ratio of purines to pyrimidines in the precipitate is 


Iny.i(M —1)+ z. nnN 
Re N=M : 





fnuaa + ba Nn 
N=M 
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The expressions for F and R in the case of a random nucleic acid are obtained by 
substituting ny from equation (1) into equations (2) and (3) and simplifying: 


e-*, T{MRM—' —-, Ry"(N +1) | 
Pecukutaid ; 4 
Ro + | te el (4) 


Cy N 
fRoM-\(M — 1) + (M1)! YO a - 
Wie a eli e a eee Y at ts (5) 
meoitiu—-y sy & 
vem N! 
where Ao, the ratio of purines to pyrimidines in the undigested polymer, is equal to 
1.19 for tobacco mosaic virus nucleic acid. 

Although F and FR appear to depend on the two parameters f and M, they can be 
uniquely determined by a single parameter P = M + f. This is true because, for 
any such number P, there is only one combination of an integer (i.e., /) and a frac- 
tion (i.e., f) which when added together will yield that number. Therefore, for 
each assumed value of P, the corresponding values of F and R can be computed 
and by repeating the calculations for several P values, a plot can be obtained in 
which R is given as a unique function of F. Thus, although the exact nature of the 
selective precipitation is not known, it is possible to predict for the random nucleic 
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Fic. 1.—A comparison of the ribonuclease digest of tobacco mosaic virus nucleic acid (“experi- 
mental’’) with that of a randomly arranged ribonucleotide polymer of the same over-all base com- 
position. The curve indicates the fraction, F, of the latter digest which would be precipitated 
by trichloroacetic acid and the ratio, R, of purines to pyrimidines in the precipitate for various 
values of the precipitation parameter, P. The tendency of the experimental points to fall farther 
from the curve with increasing digestion time probably reflects a “spontaneous’’ cleavage of 
nucleic acid at sites other than those known to be attacked by ribonuclease. 
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Fic. 2.—The ribonuclease digestion of tobacco mosaic virus nucleic acid. The curve 
represents the fraction, F’, of the original material which can be precipitated by trichloro- 
acetic acid after various periods, ¢, of digestion. Note that there is a rapid initial loss, 
followed by a slow loss, which is believed to be due to the spontaneous breakdown of 
digest fragments. The extrapolated value at ¢ = 0 is presumed to be that which would 
be obtained if the digestion went entirely according to the known specificity of ribo- 
nuclease. 





acid with Ry = 1.19 a relation between the fraction of the digest precipitated and 
the ratio of purines to pyrimidines in the precipitate. 

Comparison of Theoretical and Experimental Results——The theoretical curve of 
R plotted against F for the random nucleic acid is compared in Figure 1 with the 
experimental points determined by Reddi and Knight for tobacco mosaic virus 
nucleic acid. Some rather large discrepancies will be noted: These may be at- 
tributed to one or more of the following: (1) a distribution of purines between suc- 
cessive pyrimidines which does not follow the Poisson law; (2) a breakdown of 
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nucleic acid in some manner other than according to the known specific action of 
pancreatic ribonuclease; (3) a manner of selective precipitation greatly different 
from that which has been assumed here. Evidence will now be given which indi- 
cates that the discrepancies between theory and experiment may be due primarily 
to the second of the possible causes listed. 

The ‘‘Spontaneous” Breakdown of Nucleic Acid.—In Figure 2 the fraction of the 
viral nucleic acid digest precipitable by trichloroacetic acid is plotted as a function 
of digestion time. It will be seen that there was at first a rapid loss of precipitable 
material which was essentially completed within 6 hours. This was followed by 
a slower loss, which continued for the full period of observation without any appar- 
ent deceleration. Analogous trends are seen in Figure 3, where the ratio of purines 
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Fic. 3.—The ribonuclease digestion of tobacco mosaic virus nucleic acid. The curve 


represents the ratio, R, of purines to pyrimidines in the material precipitable by tri- 
chloroacetic acid after various periods, ¢, of digestion. As in Fig. 2, there is a rapid 
initial change, followed by a slow change. The extrapolated value at ¢ = 0 is again pre- 
sumed to be that which would be obtained in the absence of the spontaneous breakdown. 


to pyrimidines is plotted against digestion time. Both curves suggest strongly that 
more than one mechanism was acting during the digestion period. The investiga- 
tors reported also that when the nucleic acid was incubated in the absence of any 
known enzyme, there was a slow but appreciable loss of trichloroacetic acid—pre- 
cipitable material. (The ‘‘spontaneous” breakdown of similar preparations has 
also been noted by Cohen and Stanley,’ Northrup and Sinsheimer,® and others.) 
These results would suggest that in Figure 2 only the rapid, initial loss of precipi- 
table material is due to the specific action of the ribonuclease and that the slower 
loss must be attributed to some other manner of breakdown, With this assumption. 
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the curves in Figures 2 and 3 have been extrapolated to zero digestion time in such 
a way that the trend due to the slow breakdown is continued and the trend due to 
the ribonuclease digestion is ignored. The extrapolated values are then presumed 
to be the values of F and R which would be obtained for the core material if the 
breakdown were due entirely to the known specific action of ribonuclease. In Fig- 
ure 1 it will be seen that the point corresponding to these values lies very close to 
the theoretical curve for the random nucleic acid. The extrapolated point almost 
coincides with the calculated point for P = 3, which would indicate that all the 
digest fragments consisting of 4 or more nucleotide residues, and none smaller, are 
precipitated as core. 

As the digestion period is increased, the experimental points fall farther and 
farther from the theoretical curve. If it is assumed that the discrepancies between 
these points and the calculated point for P = 3 are due to a gradual degradation of 
the theoretical digest fragments, the points can then be used to test various hy- 
potheses as to the nature of the slower degradation. In Figure 4 the data are com- 
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Fic. 4.—Plot similar to Fig. 1. Here the experimental points 
from Fig. 1 are compared with curves calculated for three different 
hypotheses as to the manner in which the fragments after ribo- 
nuclease digestion break down spontaneously. Hypothesis 3, which 
comes closest to fitting the data, postulates that by some unknown 
mechanism only the terminal pyrimidines are removed from the 
digest fragments. 


pared with the curves calculated from three such hypotheses, which are as follows: 
(1) random breakage of phosphoester bonds—that is, within any specified interval 
of time the probability of a given bond being broken is the same for all such bonds; 
(2) random breakage (as defined above) of only those bonds connecting the ter- 
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minal residues to the rest of the chain; and (3) random breakage of only those 
bonds connecting the terminal pyrimidines to the rest of the chain. While none of 
these curves are in good agreement with the point, it will be noted that hypothesis 
3 comes closest to fitting the data. This observation suggests that there is some 
specificity involved in the slower breakdown, though no chemical explanation for 
removal of terminal pyrimidines from the digest fragments can be offered at this 
time. 

Conclusions.—The experimental results for the digestion of tobacco mosaic virus 
nucleic acid appear to be consistent with those which would be obtained with a 
randomly arranged ribonucleotide polymer of the same over-all base composition 
as the viral nucleic acid. This finding does not imply the lack of a definite sequen- 
tial arrangement of residues in the viral nucleic acid. It does suggest, however, 
that the arrangement consists of something more complex than the mere repetition 
of one simple pattern involving only a few residues. For this reason, ribonuclease 
digestion appears to be quite limited as a method for comparing the structure of 
nucleic acids. However, this limitation may be due mostly to the lack of good 
methods for separating digest fragments. If it were possible, for example, to ana- 
lyze only those fragments containing 7 or more residues, the method might become 
very useful. From equation (1) it can be calculated that a nucleic acid of molecular 
weight 10° would be expected to yield only 1 or 2 heptanucleotides. Such frag- 
ments would then be more characteristic of the molecules from which they were ob- 
tained and might provide a means of distinguishing or identifying ribonucleic acids. 

Summary.—The pancreatic ribonuclease digest of tobacco mosaic virus nucleic 
acid is compared with that of a hypothetical nucleic acid in which the same relative 
numbers of purine and pyrimidine residues are arranged in a long random se- 
quence. Specifically, the values determined by Reddi and Knight for the fraction 
of the viral nucleic acid digest precipitable by trichloroacetic acid and the ratio of 
purines to pyrimidines in the precipitate are compared with the corresponding values 
that would be expected for the random nucleic acid. The agreement is excellent if 
allowance is made for a slow side effect in the viral nucleic acid digestion and if it is 
assumed that just those digest fragments containing four or more nucleotide resi- 
dues are precipitated by trichloroacetic acid. The experimental data are also used 
to test some hypotheses as to the nature of the slow side effect accompanying ribo- 
nuclease digestion. 


I am grateful to Drs. K. K. Reddi and C. A. Knight for allowing me to use their 
results before publication. 
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MOLECULAR BASIS OF BIOLOGICAL STABILITY TO 
HIGH TEMPERATURES* 
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WEST LAFAYETTE, INDIANA 


Communicated by A. K. Balls, April 1, 1957 


One of the most plausible explanations for the ability of thermophilic organisms 
to grow at temperatures that usually are lethal to other organisms is that they 
possess relatively heat-stable proteins. Substantial support for this hypothesis 
recently has been furnished by Koffler and his associates," ? who showed that 
cytoplasmic proteins obtained from thermophilic bacteria are strikingly more heat- 
stable than those from mesophiles. Their experiments extend previous reports*~7 
that the activity of certain enzymes from thermophiles is more heat-stable than that 
of similar preparations from mesophilic cells. Unfortunately, the available data, 
except Campbell’s,® were obtained with impure enzymes and thus do not enable one 
to distinguish between inherent thermostability and external stabilization. Camp- 
bell’s experiments, however, were performed with purified enzymes. He showed 
that crystalline amylase isolated from cultures of Bacillus coagulans grown at 55° 
C. is markedly more heat-stable than a similar isolate from cultures grown at 
35° C. Though one could argue that even crystalline proteins may contain sta- 
bilizing contaminants in minute yet effective concentrations, the most reasonable ex- 
planation for this finding is that the observed difference in relative thermostability 
is due to molecular differences and not to external substances. Campbell’s observa- 
tion confirms our own conclusions regarding the relative thermostability of flagella 
from thermophiles. These experiments were presented in preliminary accounts 
before® § and are described more fully in this paper. 

Bacterial flagella are helical organelles made up of fibrous proteins’—" recently 
named “‘flagellins.”!! The characteristic appearance of flagella as observable in 
the electron microscope can be destroyed by acid,® ™ ™ 1® heat,” § urea,” '” and 
acetamide.? Apparently flagella ‘‘dissociate” into their component units or multi- 
ples thereof; this seems true for the action of acid below pH 3-4* !* ™ and is prob- 
ably true also for the effect of the other agents mentioned. Suspensions of flagella 
are viscous, but their viscosity decreases markedly as the flagella dissociate. Such 
a decrease therefore can be used as a measure of stability. The advantage of using 
flagella as a test system is that highly purified preparations from many different 
organisms can be made rather easily by the same method involving fractional 
centrifugation (with or without fractional precipitation with ammonium sulfate). 
Moreover, because flagella are organized aggregates of protein molecules, they lend 
themselves to a study of intermolecular as well as intramolecular stability. 


MATERIALS AND METHODS 


Organisms.—The mesophilic bacteria used in this study were the following: 
Proteus vulgaris Purdue M, Escherichia coli ATCC 26, Serratia marcescens 8 UK, 
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Bacillus megaterium KM, and B. subtilis 712. The following thermophilic strains 
of the genus Bacillus were used: B. stearothermophilus NCA 2184, B. sp. Purdue 
CD, B. sp. Texas 11330, and B. sp. Nebraska 39. We are grateful to various 
colleagues in other laboratories for making some of these cultures available to us. 

Isolation of Flagella—The above organisms were usually grown in a medium con- 
sisting of the following (in gm/liter): gelysate, 10; trypticase, 5; yeast extract, 1; 
dextrose, 2.5; sodium chloride, 5; dibasic potassium phosphate, 2.5. The medium, 
in 125-ml. quantities, was placed in 500-ml. Erlenmeyer flasks, sterilized, and, after 
cooling to the appropriate growth temperature, inoculated with 5 ml. of a broth 
culture. The cultures were incubated until they had reached the end of the loga- 
rithmic growth phase or the early part of the maximum stationary phase (8-12 
hours). Mesophiles were grown at 20° C., thermophiles at 60° C. Aeration was 
increased by agitation on reciprocating shaking machines. 

The cells were harvested by centrifugation at 3,900  g for 30 minutes (100-ml. 
bottles; 4,600 rpm; International Refrigerated Centrifuge, Model PR-1). Usually 
6 liters of culture was handled in one experiment. The cells were resuspended in 
400-600 ml. of 0.9 per cent saline solution and washed by centrifugation. After 
the cells had been washed, they were resuspended in 200-300 ml. of cold saline 
solution and placed in a 500-ml. ground-glass-stoppered bottle that contained ca. 
100 gm. of frozen saline solution. The bottle was placed on a shaking machine 
(“Miracle Paint Rejuvenator’’) and shaken at a rate of 620 1'/,-inch strokes/min 
for 20 minutes. This step broke off most of the flagella from the cell bodies. The 
shaken suspension was then centrifuged twice at 3,900 * g for 45 minutes at 0°- 
4° C. The centrifugate was discarded each time. The supernatant liquid was 
then brought to 0.1 saturation with ammonium sulfate and allowed to stand for 
8-12 hours. The precipitate was sedimented at 3,900 < g for 45 minutes and then 
discarded. The supernatant liquid was brought to 0.5 saturation with ammonium 
sulfate and allowed to stand for 8-12 hours. The precipitate was concentrated by 
centrifugation at 3,900  g for 45 minutes and then centrifuged at 25,000 X g for 
5 minutes (19,000 rpm; high-speed attachment of the International Refrigerated 
Centrifuge). The supernatant liquid was discarded. The surface of the pellet 
and the inside walls of the tube were washed with distilled water. By discarding 
the 0-0.1 and 0.5—1.0 fractions, one eliminates a large portion of nonflagellar ma- 
terials. The precipitate was then resuspended in 40 ml. of distilled water and 
centrifuged twice at 25,000 x g for 30 minutes; the centrifugate was discarded 
each time. The supernatant liquid was centrifuged twice at 59,000 X g for 35 
minutes (30,000 rpm; No. 40 rotor, Spinco Model L Preparative Ultracentri- 
fuge). The purified flagella were stored at 0°-4° C. until used. A second method 
of growing cells and isolating flagella was used for some of the experiments.'* How- 
ever, inasmuch as the procedural variations did not seem to influence the results, 
details regarding these variations are omitted. 

Determination of Purity —To determine dry weights, the samples were placed in 
tared glass or aluminum-foil dishes, dried at 110° C. for 4 hours, and weighed on 
a microchemical balance. Total nitrogen was determined according to the method 
of Willits and Ogg.'* Carbohydrate (as hexose) was determined by the anthrone 
method of Loewus.” Nucleic acid was determined according to the method of 
Warburg and Christian;*! a correction for nonspecific dispersion of radiant energy 
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was made as described by Weibull.’ As mentioned previously, flagella dissociate 
below pH 3-4. Flagella thus treated are no longer centrifugable at 59,000 x g 
for 30 minutes. It was found that the sediment obtained under such conditions 
contains the bulk of nonflagellar impurities.* © The absence of such material is 
perhaps the most informative single criterion for the purity of isolates. In this 
study the routine procedure was exposure of the flagella to pH 2 at 26° C. for 30 
minutes, followed by centrifugation at 59,000 * g for 30 minutes (“pH 2 insoluble 
material”). The pellet, if any, was transferred to a tared glass dish and weighed 
as described. Usually the preparations used were estimated to be approximately 
95-99 per cent pure. The ranges of values encountered for preparations from the 
eight organisms employed in this study were as follows: total Kjeldahl nitrogen, 
15.8-16.4 per cent; carbohydrate, 0.1-0.6 per cent; nucleic acid, 0.1-0.5 per cent; 
‘pH 2 insoluble material,’’ 0.2—0.7 per cent. 

Determination of Viscosity.—To test the stability of flagella to heat, 2 ml. of a 
suspension was placed in a viscometer of the Ostwald type. The initial viscosity 
was determined at 20° + 0.05° C. after a period of equilibration of 1 minute. The 
instrument was then immersed in a water bath at temperatures varying from 40° 
to 80° C. Measurement of exposure time was begun after an equilibration period 
of 1 minute. After an interval of time, the viscometer was removed, quickly cooled 
to 20° C., and the flow time of the suspension was redetermined. The heating 
then was continued in increments of 5 or 10 minutes until the total exposure time 
to a given temperature was 30 minutes. To test the stability of flagella to various 
chemical] agents, 1 ml. of suspension was placed in a viscometer, and the instrument 
was placed in a water bath at 20° C. One milliliter of the solution (at 20° C.) 
containing twice the desired final concentration of the compound under study was 
added to the suspension of flagella. The flow time of the mixture was determined 
at intervals over the course of an hour or longer. 

The results were expressed in terms of relative viscosity, nre1 = n/no; specific 
viscosity, sp = Nrei — 1; or reduced viscosity, tre1 = 1 nsp»/C, where 7 is the viscosity 
of the solvent plus flagella, no is the viscosity of the solvent, and C is the concentra- 
tion of flagella in mg/ml. The ratio n/no was assumed to be equal to the ratio of the 
flow time of the flagellar suspension to that of the solvent. The relative viscosity 
of flagellar suspensions in the concentrations employed was not significantly 
affected by the velocity gradient, changes in pH from 4.4 to 9.1, and small con- 
centrations of sodium chloride (0-0.15 1). 

RESULTS 

Relative Thermostability of Thermophile Flagella——The effect of heat on the sta- 
bility of flagella (usually 0.5 mg/ml) was determined. Decrease in the viscosity 
of flagellar suspensions at various temperatures was used as a measure of thermal 
destruction.2 The products of this “disintegration” (dissociation) are nondialyz- 
able and are precipitable by ammonium sulfate, but they cannot be sedimented 
at the RCF sufficient to remove whole flagella. Figure 1 shows the effect of heat 
on the flagella from E. coli, a mesophile. Figure 2 illustrates the stability of 
flagellar from B. sp. 11330, a thermophile. From experiments similar to these, 
the data summarized in Table 1 were obtained. Table 1 shows the decrease in 
the specific viscosity of flagellar suspensions after exposure to various tempera- 
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Fic. 1.—Effect of various temperatures on the reduced viscosity 
of suspensions of flagella from the mesophile Escherichia colt. 
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Fia. 2.—Effect of various temperatures on the reduced viscosity of 
suspensions of flagella from the thermophile Bacillus sp. 11330. 
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tures for 20 minutes. Disintegration of flagella from mesophiles becomes marked 
at temperatures higher than 50° C., while flagella from thermophiles remain 
intact at temperatures as high as 70° C. This is in agreement with the ex- 
perience of McCoy,?? who showed that the H-agglutinogen of Clostridium thermo- 
saccharolyticum, a thermophile, is much more heat-stable than that of C. butyri- 
cum, a mesophile. In general, the heat stability of the flagella used in the 
present study corresponds to the ability of the organisms from which they were 
isolated to grow at high temperatures. For example, the strain of B. subtilis used, 
for which the maximum growth temperature (52° C.) is between that of the 
mesophiles and that of the thermophiles, produces flagella which are intermediate 
in heat stability. 

Inability To Demonstrate Factors That Protect Flagella from Destruction by Heat.- 
The question whether the relative heat stability of thermophile flagella is an in- 
herent property or is provided by protective materials?’ was considered in a num- 
ber of experiments. In the first, a mixture of equal amounts of flagella from P. 
vulgaris and B. sp. CD was exposed to 60° C. for varying lengths of time. The 
mixture behaved essentially as one would have expected on the basis that the 


TABLE 1* 
THERMOSTABILITY OF FLAGELLA FROM MESOPHILIC 
AND THERMOPHILIC BACTERIA 
Per Cent Initia Speciric Viscosiry arrer 20 Min. at °C. 
50 6 65 70 75 


FLAGELLA FROM 
P. vulgaris 99 67 18 
Y. cola 100 56 36 
B. megaterium 100 72 28 


B. subtilis 98 79 31 fe cf 
B. stearothermophilus 2184 te : ay ie 43 18 
B. sp. CD i a 99 94 41 
B. sp. 11330 100 ; 105 Sa 96 79 5 
B. sp. 39 100 ; 112 ; 94 56 14 


Pi + a were taken every 5-10 minutes over a period of 30 minutes. Viscosity measurements were made 
at 2 . 

flagella from the mesophile disintegrated while the flagella from the thermophile 
remained intact. From kinetic data the following values for the specific viscosity 
of the flagellar suspensions are given as examples: at zero time—0.059 (meso- 
phile), 0.061 (thermophile), 0.120 (mixture calculated), and 0.115 (mixture ob- 
served); at 30 minutes—0.010 (mesophile), 0.059 (thermophile), 0.069 (mixture 
calculated), and 0.061 (mixture observed). <A similar situation was observed 
when a mixture of flagella from B. subtilis and B. sp. CD was exposed to 65° C. 
If the relative instability of the mesophile flagella were due to external labilizing 
factors,”® such factors could have been released when the mesophile flagella disin- 
tegrated and might have affected the flagella from the thermophile. Similarly, 
the presence of protective materials” in the thermophile, if diffusible, might have 
left its mark on the stability of the mesophile flagella. The possibility exists, 
of course, that any labilizing or stabilizing factor may be specific for a given prep- 
aration of flagella under examination, may be firmly bound or unable to reach the 
necessary sites, and thus may escape detection by this method. However, the 
simplest explanation for the previous and the following experiments is that the 
relative stability of flagella depends upon the molecular architecture of the flagella 
themselves rather than on extraneous materials. 
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Cell-free extracts of B. sp. CD were prepared by oscillation of a suspension of 
this organism in distilled water for 90 minutes in a Raytheon 9-ke. Magneto- 
striction Oscillator. The debris was centrifuged at 25,000 X g for 1 hour. The 
clear supernatant liquid, containing 40 per cent of the total dry weight of the 
cells, was concentrated to dryness under vacuum at 23° C. and stored at 0°-4° C. 
until needed. The cell-free extract and flagella from B. subtilis in concentrations 
of 0.5 mg/ml were mixed, and the stability of the flagella under these conditions 
was compared to their stability in the absence of thermophile material. The 
cytoplasmic constituents isolated from the thermophile exerted no appreciable 
protective effect. 

Flagella from B. sp. CD were disintegrated by exposure to 85° C. for 15 minutes, 
to pH 2 for 30 minutes at 26° C., or to sonic oscillation (Brush Development Com- 
pany Hypersonic-Transducer, Model Bu-101-A, 15 minutes). It was hoped that 
a protective factor contained in the thermophile flagella would be released but not 
destroyed by one or more of these treatments. The preparations disintegrated by 
heat or acid (after neutralization) were mixed in equal amounts (usually 0.5 mg/ml) 
with flagella from B. subtilis. The sonically disintegrated thermophile flagella were 
mixed with flagella of P. vulgaris (0.5, 1.0, 1.5, and 2.0 mg/ml of extract with 1.0 
mg/ml of flagella; also, 3.0 and 6.0 mg/ml of extract with 3.0 mg/ml of flagella). 
No significant stabilization of the flagella from B. subtilis and P. vulgaris was ob- 
served at 65° and 55° C., respectively. Also, the effect of calcium ions, previ- 
ously implicated in the stabilization of proteins,?4—* was tested (equal volumes of 
0.05, 0.1, 0.5, and 1.0 M solutions of calcium chloride and suspensions of flagella 
from B. subtilis, 0.5 mg/ml; 65° C.). No stabilizing effect could be demonstrated.” 

The possibility of adsorbing a hypothetical protective factor was tested by de- 
termining whether thermophile flagella treated with various adsorbing agents would 
become less heat-stable.27. Two milliliters of a suspension of flagella from B. sp. CD 
(3.0 mg/ml) was incubated for 15 minutes at 26° C. in the presence of Dowex 50 
(H form; 20 and 40 mg/ml), Dowex 2 (OH form; 20 and 40 mg/ml), or glass beads 
(3 mm, 5 beads/ml); the suspensions were shaken gently. The adsorbing agents 
were removed by centrifugation, and the flagella in the supernatant liquid were 
examined for heat stability at 60° C. Even after treatment with exchange resins 
or beads, the flagella remained unaffected by that temperature. 

Enzymatic removal of a possible protective material was attempted by exposure 
of flagella from B. sp. CD to the action of trypsin, ribonuclease (RNAse), and 
desoxyribonuclease (DNAse). Trypsin was assayed by a method described by 
Crewther,”> and RNAse and DNAse by the method of Roth.*® To eliminate the 
possibility that the flagella might be inhibitory, the activity of these enzymes 
on standard substrates (gelatin, yeast RNA, and sperm DNA) was determined 
in the presence and in the absence of flagella. The presence of the flagella had no 
effect on these enzymes. Sufficient enzyme was added to flagellar preparations to 
hydrolyze a weight of substrate that was from 20 to 40 times greater than the total 
amount of flagella present. After incubation with trypsin, RNAse, and DNAse 
at various temperatures and for varying lengths of time, the thermophile flagella 
were just as thermostable at 65° C. for 20 minutes as they had been before enzy- 
matic treatment. 

The Unlikelihood of Disulfide Linkages Being Responsible for the Relative Stability 
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of Thermophile Flagella.—The possibility that thermophile flagelia are more stable 
because of more effective disulfide bonding did not seem likely, because of the fact 
that the cysteine-cystine contents of mesophile and thermophile flagella are low 
and are similar regardless of the source (0.54-0.75 per cent).‘‘ Moreover, the 
rapid dissociation of flagella below pH 3-4 speaks against the existence of inter- 
molecular disulfide bridges. Nevertheless, the effect of sodium thioglycolate, an 
agent known to break —S—S— bridges, on the stability of flagella was deter- 
mined. Suspensions of flagella from B. stearothermophilus 2184 (ner, 1.23) and B. 
sp. CD (ne:, 1.30) were mixed with thioglycolate in the presence of urea (final 
concentrations: 0.2 M thioglycolate and 5 M urea in M/60 phosphate buffer, 
pH 8). Under the same conditions a 5 per cent solution of ovalbumin was appre- 
ciably denatured (after 24 hours at 26° C. the solution was diluted fivefold and the 
pH adjusted to 4.8; over 25 per cent of the ovalbumin was precipitated, while only 
5 per cent of it precipitated when thioglycolate was absent). Flagellar suspensions 
were incubated for varying lengths of time, and their stabilities were determined by 
viscosity measurements. Thioglycolate did not affect the viscosity of these sus- 
pensions, even after contact for 110 hours. 

Evidence for More ‘Effective’ Hydrogen Bonding in Thermophile Flagella.—If 
the relative thermostability of flagella from thermophilic bacteria is due to hy- 
drogen bonding, 2* *° they should be more resistant than mesophile flagella to com- 
pounds regarded as hydrogen-bond breakers, such as urea and acetamide. The 
effect of urea on the viscosity of a suspension of flagella from B. subtilis is shown in 
Figure 3. To disintegrate thermophile flagella to the same extent, one needs to 
increase the concentration of urea to 9 MV; at a level of 6 1/7, urea does not exert any 
observable effects on thermophile flagella. This is shown in Table 2, in which the 
relative stabilities of four mesophile and four thermophile preparations are com- 
pared. Disintegration by urea, like that by heat, is accompanied by the loss of 
the characteristic flagellar structure; similarly, the products are nondialyzable and 
cannot be sedimented at the RCF at which whole flagella can be centrifuged. 
The data indicate that the thermophile flagella are not only more heat-stable than 
are the mesophile flagella but also more resistant to the action of urea. They are 
also more resistant to 10 M acetamide, as is shown in Table 2. The experiments 
with acetamide were performed in a manner similar to those previously discussed, 
but for the sake of brevity only end-point data are presented. 

Inability To Demonstrate Factors That Protect Flagella from Disintegration by Urea 
or Labilize Them to Its Action.—Evidence for the existence of possible stabilizing 
factors in the thermophile flagella from B. stearothermophilus 2184 or labilizing 
factors in the mesophile flagella from P. vulgaris was sought in experiments similar 
in design to those previously performed. 

Cell-free extracts were prepared as before, except that the cells were exposed to 
oscillation for 30 instead of 90 minutes and the debris was removed by three con- 
secutive centrifugations of 60 minutes each. The ash of whole cells was obtained 
by the method of Bates et al.*! The addition of cytoplasmic constituents from 
the thermophile (2.5 mg/ml) to the mesophile flagella (2.0 mg/ml) or from the 
mesophile (2.4 mg/ml) to the thermophile flagella (1.2 mg/ml) did not markedly 
alter the stability of the flagellar preparations to 5 M urea. No stabilization or 
labilization was brought about when thermophile ash was added to mesophile 
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Fic. 3.—Effect of various concentrations of urea on the reduced 
viscosity of suspensions of flagella from the mesophile Bacillus sub- 
tilis. 


TABLE 2* 
STABILITY OF FLAGELLA TO UREA, ACETAMIDE, AND 
Soprum DoprEcyLsuLFraTE (SDS) 
Per Cent Initia, Repucep Viscosiry in PrResENcE or— — 
0.0035 M 
10 M Acet- Dodecyl- 
6 M Urea amide sulfate 
-——-After Minutes = -————— 
60 30 60 
67 


FLAGELLA FROM 5 30 
. vulgaris 4 
. colt 16 ; 
. megaterium 23 7 60 
. subtilis 83 72 4: 21 F 63 
. stearothermophilus 2184 100 100 100 100 36 103 
. sp. CD 104 108 111 114 100 
. sp. 11330 100 100 100 100 \)- 100 
. sp. 39 103 103 104 107 2 ste 


l 
z 


* Samples were taken every 5 minutes over a period of 60 minutes. The data presented above were abstracted 
from experiments similar in execution to those shown in Fig. 3. 
flagella or when mesophile ash was added to thermophile flagella (all in concen- 
trations of 2.0 mg/ml). 

The products of flagella disintegrated by heat (thermophile flagella at 80° C. 
for 30 minutes; mesophile flagella at 60° C. for 30 minutes), acid (prepared as be- 
fore), and oscillation (26 minutes’ treatment in a 9-ke. Raytheon Magnetostriction 
Oscillator) did not significantly change the behavior of whole flagella obtained from 
cells of the opposite stability type. The following concentrations were used in 
these experiments (mg/ml): whole and disintegrated flagella from B.  stearo- 
thermophilus 2184, 1.2; from P. vulgaris, 2.0. 

Treatment of flagella with Dowex 50 (H form; 2 gm. moist resin mixed with 5.0 
mg. of thermophile flagella [in 4 ml] or 2.4 mg mesophile flagella [in 3 ml]; 26° 
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C. for 5 minutes), Dowex 2 (OH form, 2 gm. moist resin mixed with 5.0 mg. flagella 
fin 4 ml., thermophile; in 5 ml., mesophile]; 26° C. for 5 minutes), or small beads 
(5.0 mg. of Scotch-lite glass beads mixed with 5.0 mg. of thermophile or 5.0 mg. of 
mesophile flagella; 26° C. for 15 minutes) had no significant effect on the behavior 
of the flagella tested in the presence of 5 M urea. 

Exposure of thermophile and mesophile flagella to trypsin, RNAse, and DNAse 
left the manner in which the thermophile and mesophile flagella responded to urea 
essentially unaltered. The amount of flagella used per experiment was 4.0 mg/ml 
of thermophile flagella and 3.0 mg/ml of mesophile flagella. The enzymes were 
used in levels sufficient under the conditions of the experiments to hydrolyze a 
weight of substrate greater than the weight of flagella present. Incubation was 
at 20° C. for 10 minutes (RNAse and DNAse) or 20 minutes (trypsin). 

Mixtures of thermophile and mesophile flagella behaved in the presence of 5 M 
urea as one could have predicted on the basis of the behavior of the individual 
preparations (cf. Fig. 4). 

Differences in the Titration Curves for Mesophile and Thermophile Flagella.—Wei- 
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Fic. 4.—Effect of 5 M urea on the specific viscosity of mixed 
suspensions of flagella from the mesophile Proteus vulgaris and the 
thermophile Bacillus stearothermophilus 2184. 
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bull examined the titrimetric behavior of flagella from P. vulgaris and found that 
after disintegration of flagella at pH 2.5 the solution of flagellin acted as if additional 
titratable groups had become available.** In the intact flagellum such groups 
may not be titratable, because they are sterically inaccessible or combined in link- 
ages. However, the difference in titration curves may also be due to changes in 
the pKa values of the ionizable groups after the flagellum has broken apart. In any 
case, it became of interest to determine whether the relative stability of thermo- 
phile flagella might be due to interactions in which these groups are involved. 
Titration curves were prepared according to the method described by Weibull,** 
and those for the mesophile flagella were similar to that described for P. vulgaris. 
The apparent increase in the number of titratable groups after dissociation of the 
flagella by acid cannot be correlated with the relative thermostability of the flagella. 
On the other hand, there is a striking difference in the total number of H+ bound 
at pH 2.5 and 10.5. This is shown in Table 3.** 


TABLE 3* 


NUMBER OF TITRATABLE GROUPS IN FLAGELLIN FROM 
MESOPHILIC AND THERMOPHILIC BACTERIA 


H* Bounp/105 GM H * Bounp/Mo.e 
MIN. ——- —- ———At pH—— —-——---— - 
FLAGELLIN FROM Mot. Wr. 2.5 10.5 2.5 


P. vulgaris 19,600 +120 —63 +24 
S. marcescens 20 , 800 +111 —53 +23 
B. subtilis 18, 200 +102 —48 
B. stearothermophilus 2184 13,900 + 58 —40 
B. sp. CD 15,200 + 90 —32 
B. sp. 11330 16,700 + 67 —28 
* The values presented in the table were taken from complete titration curves prepared as done by Weibull 


(Acta Chem. Scand., 4 260, 1950). The minimal molecular weights were calculated on the assumption that each 
flagellin molecule contains one N-terminal alanine (see Koffler et al., Arch. Biochem. and Biophys., 64 509, 1956). 


These data can be interpreted in at least two ways. (1) Mesophile flagella 
possess approximately twice as many acidic and basic groups as do thermophile 
flagella. This means that the flagellin molecules constituting mesophile flagella 
potentially have a greater opportunity for electrostatic attraction and repulsion 
than do those of which thermophile flagella are composed. Greater electrostatic 
attraction is not likely in the light of our other observations, but it is possible that 
the net stability of thermophile flagella is the consequence not only of more effective 
bonding between flagellin molecules but also of lesser repulsion between them. 
(2) Thermophile flagella possess the same number of titratable groups as do meso- 
phile flagella, but these groups are more strongly acidic and basic in the former. 
For example, the carboxylate anions in thermophile flagella may be stabilized in 
such a way (probably through hydrogen bonds) as to lower the effective pKa of 
these carboxy! groups (interaction of the ammonium cation might increase the pKa 
of the amino groups).* Thus the total number of H*+ bound at pH 2.5 may be 
considerably less than the maximum acid binding capacity. Comparison of the 
maximum number of titratable groups at pH 1 or pH 1.5 therefore may reveal that 
both mesophile and thermophile flagella possess the same number. 

The differences between the titration curves of mesophile and thermophile 
flagella need further study before one can distinguish definitely between the above 
possibilities. In any case, though, it is clear that the amphoteric properties of 
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mesophile and thermophile flagella differ remarkably and that these differences 
probably constitute useful clues regarding their architecture. 

The Relative Stability of Thermophile Flagella to Sodium Dodecylsulfate.—Thermo- 
phile flagella apparently are more resistant to the action of sodium dodecylsulfate 
(SDS) than are mesophile flagella. For example, 3.5 X 10~* M of SDS in M/60 
phosphate buffer (pH 8) caused the partial disintegration of mesophile flagella but 
had no effect on thermophile flagella (cf. Table 2); to bring about the same degree 
of dissociation of thermophile flagella, the concentration of SDS had to be raised 
to8.7 X 10-* M. 

The mechanism by which anionic detergents denature proteins is not completely 
understood, but it is thought that they first react with basic groups.”* ‘The greater 
resistance of thermophile flagella to the action of SDS, therefore, may be merely an 
expression of the fact that they possess fewer accessible basic groups than do meso- 
phile flagella. On the other hand, in the above experiment the observed dissoci- 
ation of the flagella did not occur until the molar concentration of SDS was greater 
than the molar concentration of titratable basic groups (for instance, the flagella 
from B. subtilis used were 6.25 X 10~* M with respect to titratable basic groups, 
and the flagella from B. stearothermophilus were 5.2 X 10-4 M). Inasmuch as the 
nonpolar moiety of SDS might be expected to have an affinity for the nonpolar 
residues of the flagella,® it is possible that the dissociation of the flagella is caused 
by a disruption of hydrophobic bonds. This may be true especially for the thermo- 
phile flagella, the observed dissociation of which was brought about by concentra- 
tions of SDS considerably in excess of those that are likely to be bound electro- 


statically. The greater stability of thermophile flagella to SDS therefore may be 
a result of more “effective”? hydrophobic bonding, but. this interpretation needs to 
be examined further. 


DISCUSSION 

The remarkable ability of thermophile organisms to grow at high temperatures 
may be conferred by one or both of the following features: (1) the relative heat 
stability of their cellular components, a stability that may be either inherent or 
provided by protective materials, and (2) the rapid resynthesis of essential molecules 
damaged by heat.‘ Various types of substances likely to be present in the cyto- 
plasm of organisms have been found to protect proteins from disorganization by 
heat,?*> among them calcium and other metallic ions,?*~** sodium desoxyribonu- 
cleate,® certain proteins,** and cellular particles.“ Certain impure, relatively 
heat-stable enzyme preparations have been shown to lose their thermostability on 
purification.** These findings make it imperative that studies regarding the nature 
of biological thermostability include a diligent search for protective materials. 

Koffler and Gale showed that the striking relative heat stability of cytoplasmic 
proteins from thermophiles was not significantly diminished when “small par- 
ticles” and nucleic acids were removed.!| However, one could contend that a 
small amount of these materials remained, sufficient for protection. Perhaps more 
convincing were their experiments, in which cytoplasmic proteins from mesophiles 
and thermophiles were mixed and the heat stability of the mixture could be ac- 
counted for by the stabilities of the individual preparations. Campbell’ showed 
that an amylase isolated from a facultatively thermophilic bacterium grown at 
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55° C. remained relatively heat-stable even when purified. to the point of crystal- 
lization; a similar preparation obtained from cultures grown at 35° C. was not 
heat-stable. This work constitutes a valuable contribution on two accounts. 
First, the use of crystalline enzymes increases the likelihood that intrinsic stability 
rather than extrinsic stabilization was observed. Second, the use of the same 
organism for the production of relatively heat-stable and relatively heat-labile amyl- 
ase decreases the possibility that the correlation between the stability of the en- 
zyme and the growth temperature of the organism is fortuitous. His observations 
confirmed our earlier preliminary account,® more fully developed in this paper, re- 
garding the heat stability of purified flagella. 

We have shown that flagella from thermophiles are much more stable than 
flagella from mesophiles. Though it is impossible to state categorically that the 
observed relative stabilities are not brought about by protective or labilizing 
materials, the evidence strongly supports such a view. A persistent search for 
the existence of such substances has yielded no positive information. 

The following facts speak against the importance of —S—S— linkages in holding 
flagellin molecules together in the flagellum. (1) Flagella contain only small 
amounts of cysteine-cystine, regardless of whether they have been produced by 
mesophiles or by thermophiles; each molecule of flagellin contains probably only 
one molecule of cysteine or cystine.'* (2) Flagella dissociate under relatively mild 
conditions, for example, below pH 3-4.* 14% 16 (3) Thioglycolate in a con- 
centration known to denature other proteins does not affect the stability of flagella; 
however, in the native flagellum, the —S—S— bridges may not be accessible. 

Since thermophile flagella are more resistant than mesophile flagella to urea and 
acetamide, agents thought to break hydrogen bonds,'* ® it seems that more “ef- 
fective” hydrogen bonding is involved in the structural stability of thermophile 
flagella. ‘Effectiveness’ may depend on the number, strength, or location of the 
bonds (or on a combination of all three factors). 

The observation that thermophile flagella might contain only half the number 
of titratable groups found in mesophile flagella suggests another explanation for 
the relative stability of the former. Flagellin molecules in the flagellum may face 
sites carrying charges of like sign. Because of the resulting electrostatic repulsion, 
such molecules would have a tendency to dissociate. The main difference between 
mesophile and thermophile flagella therefore may reside not in differences in bond- 
ing but in a greater or lesser tendency of the constituent flagellin molecules to repel 
each other. 

Thermophile flagella are also more resistant than are mesophile flagella to sodium 
dodecylsulfate, an anionic detergent. This may be because thermophile flagella 
possess fewer basic groups to which the SDS can be anchored than do mesophile 
flagella. However, the amount of SDS necessary to cause the observed disintegra- 
tion of flagella is in excess of that theoretically needed to neutralize the available 
basic groups, and the dissociation of flagellin molecules actually may come about 
through the breakage of hydrophobic bonds by the detergent. Hydrophobic bond- 
ing is thought to be due to the tendency of relatively nonpolar amino acid residues 
to avoid the aqueous phase and to adhere to each other.* Such bonding may play 
an important role in holding the flagellin molecules together in the flagellum. The 
nonpolar portion of SDS is more likely to stick to the nonpolar residues of the protein 
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than to enter the aqueous phase. This may result in the rupture of existing hy- 
drophobic bonds. The relatively greater resistance of thermophile flagella to SDS 
therefore may be due to stronger, more numerous, or more strategically located hy- 


drophobic bonds than occur in mesophile flagella. 

In addition to the explanations mentioned, one has to consider the possibility 
that the stability of thermophile flagella is provided by novel types of bonding. 
No information regarding this has been discovered so far. 


SUMMARY 


Flagella from thermophilic bacteria are much more stable than similar struc- 
tures from mesophilic bacteria. Since flagella are organized aggregates of fibrous 
proteins, this relative stability probably is a reflection of the forces holding these 
molecules together. The relative stability of flagella appears to be an intrinsic 
property of these structures, since no stabilizing or labilizing materials could be 
demonstrated. Thermophile flagella are also more stable than mesophile flagella 
in the presence of urea and acetamide, agents regarded as hydrogen-bond breakers, 
and of sodium dodecylsulfate, an anionic detergent. These findings are inter- 
preted to mean that the relative stability of thermophile flagella is due to more 
numerous, stronger, or more strategically located hydrogen and/or hydrophobic 
bonds than are present in mesophile flagella. The latter seem to contain twice as 
many titratable groups as do the former. This suggests that under certain cir- 
cumstances they might be under greater stress, created by the repulsion between 


charged groups having like signs. 


* When this work was started, it was supported in part through a contract with the Office of 


Naval Research. Since 1954 it has been aided by a grant from the United States Public Health 


Service. 

+ This paper is dedicated to my colleague and mentor, Professor P. A. Tetrault, whose interest 
in thermophilic bacteria proved contagious, as an expression of my deep gratitude for his friend- 
ship. 
t Public Health Service Predoctoral Fellow (National Microbiological Institute), 1952-1954. 
Present address: Lilly Research Laboratories, Eli Lilly and Company, Indianapolis, Indiana. 
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It has been noted! that watery solutions of fluorescent dyes profoundly alter their 
optical behavior on freezing: their fluorescent emission disappears. On further 
lowering of the temperature, light emission reappears, but the emitted light differs 
from the earlier fluorescent light in having a longer wave length and coming from an 
excitation having a longer lifetime, answering, thus, to the formal definition of 
‘“‘nhosphorescence.”’ This change in properties is not limited to dyes. Certain 
proteins, amino acids, coenzymes, and steroid hormones, for instance, show anal- 
ogous behavior. 

With regard to its possible biological implications, this phenomenon deserves 
further study. Two questions suggest themselves: (1) Does the phosphorescent 
light originate from a triplet excitation or from a crystal phosphor? (2) Is this 
behavior a universal property of fluorescent substances? The experiments reported 
in this paper suggest answers to these questions. 

There are several criteria which can be utilized in studies of luminescence to 
characterize an emission as arising from a metastable triplet level. 


1. Triplet emission is shifted toward the longer wave lengths relative to the 
fluorescent emission. 
2. Triplet emissions have a much longer lifetime than singlet emissions. 
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Triplet-state emissions obey first-order decay kinetics. 

Triplet-state emissions are usually only slightly temperature-dependent, i.e., 
activation energy is negligible. 

The molecule in the triplet state, having an unpaired electron spin, is paramag- 
netic. 

Triplet emissions induced by polarized exciting light are slightly polarized or 
nonpolarized. 

The molecule in the triplet state, being a chemical species in its own right with 
its own set of properties, will exhibit its own characteristic absorption spectrum. 


For the work reported in this paper we have used Criteria 1—4 to aid in character- 
izing the electronic states involved. 


METHODS 


CRITERION 1: EMISSION SPECTRA 
A. VISUAL 


A knowledge of the fluorescent spectrum for a particular compound permits one 
to observe visually whether a subsequent delayed emission has undergone a wave- 
length shift. For example, tyrosine fluorescence resides exclusively in the ultra- 
violet, but when the compound is frozen in the presence of glucose, it exhibits, upon 
proper excitation, a pale blue afterglow. We recognize this shift of the emission 
toward longer wave lengths as presumptive evidence for the occurrence in the 
molecule of an electronic transition different from that responsible for the ultra- 
violet fluorescence. In this study most of our observations of wave-length shift 
have been visual. 

B. INsTRUMENTAL 

1. Beckman DK-1 Spectrophotometer.—We have occasionally measured fluores- 
cent and phosphorescent spectra with the DK-1 spectrophotometer by a technique 
described by Gemmill.2. Spectra obtained by this method are subject to consider- 
able error unless the spectral response of the prism-photomultiplier assembly has 
been calibrated against a standard lamp. 

2. Leiss Monochromator—Red-sensitive Phototube Assembly.—Emission spectra 
have been obtained with a Leiss fused quartz monochromator coupled with a Du- 
Mont red-sensitive photomultiplier tube, type K 1292, sensitive to 1.3 uw. The 
photomultiplier signal was piped into a Beckman Model V micromicroammeter, 
which was reversed to take a negative input signal. Spectral emission intensity was 
measured as a linear response of the ammeter-galvanometer deflections. The 
monochromator-photomultiplier-ammeter assembly was calibrated for spectral 
distribution intensity in relative units against a standard lamp calibrated for us by 
the National Bureau of Standards. Monochromator slit widths were 0.5 mm. 
Samples under study were contained in flat-faced Pyrex cuvettes of 4-mm. inside 
width and were frozen in liquid nitrogen contained in an unsilvered Pyrex Dewar 
flask. The frozen system was illuminated from one side with a Hanovia Compact 
xenon are, and the emission from the opposite side was introduced into the entrance 
slit of the Leiss monochromator. Judicious use of liquid and glass filters permitted 
proper separation of exciting and emitted light. 
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CRITERION 2: LIFETIMES 


Since fluorescent lifetimes are of the order of 10-* second in duration, any emission 
which is apparent after removal of the exciting source, especially when accompanied 
by a shift toward longer wave lengths relative to the fluorescence emission, is in- 
dicative of a metastable-state involvement. Increased lifetimes are obvious visually 
and may be quantitated from kinetic studies. 


CRITERION 3: DECAY KINETICS 


The Beckman DK-1 automatic recording spectrophotometer was used for kinetic 
studies. With the instrument set to record energy, the wave-length drum was set 
to an appropriate wave length, and the slits were set to maximum width, namely, 
2.0mm. Aqueous solutions of the systems under study were placed in 25 X 200- 
mm. Pyrex test tubes and immersed initially in a methanol-dry ice bath for prior 
freezing. Systems studied at 194° K. were permitted to remain in this bath for 
a minimum of 15 minutes to attain temperature equilibrium. Samples to be studied 
at 77° K., following freezing in the dry-ice bath, were then transferred to Dewar 
vacuum flasks containing liquid nitrogen. Samples were allowed to remain here 
for a minimum of 15 minutes to insure temperature equilibration, following which - 
they were placed for 5 seconds in the unfiltered beam of a Hanovia high-pressure 
mercury lamp. The glowing sample was then immediately transferred into the 
lamp-housing compartment of the DK-1 and positioned at the light-entrance port. 
Following an immediate deflection response of the potentiometer, the magnitude of 
which depended on the intensity of the afterglow, the decay was tracked auto- 
matically, with the recorder chart set to roll at 12 inches per minute. Energy 
recorded in this manner on the DK-1 was read off the graph in relative linear units. 
The curves so obtained were examined for fit to kinetic equations. 

In making kinetic studies with the DK-1, it should be noted that the instrument 
itself has a finite decay time. For example, in recording on energy with a light 
source inducing a 100 per cent energy reading on the potentiometer, shutting off 
the light causes the potentiometer to “decay” to 0 per cent energy. By running 
the chart at 12 inches per minute, we were able to obtain an instrument ‘decay 
curve.” For our instrument this ‘decay curve” was found to fit a first-order kinetic 
equation and had a rate constant of approximately 1. Kinetic data obtained with 
the DK-1, therefore, represent mixtures of two rate constants, that of the instru- 
ment and that of the sample being studied. We express our thanks to Dr. William 
Arnold for deriving an equation for us which permitted the two rate constants to 
be resolved from mixed kinetic data. Using this equation we have checked the 
accuracy of our experimentally determined rate constants. The reported rate con- 
stants are accurate to within +5 per cent; most within *2 per cent usually on 
the low side. 

Ina communication from the Beckman Instrument Company concerning the use 
of the DK-1 spectrophotometer for rapid kinetic studies they pointed out that one 
should be able to determine the time constant (reciprocal of the rate constant) for 
a phosphorescent decay if it is not less than 10 per cent of the time constant of the 
DK-1. In our studies no phosphorescences were observed with rate constants 
exceeding that of our DK-1. 
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It should be emphasized that in the text figures recording our kinetic results, the 
ordinates represent phosphorescent emission intensities in relative energy units. 
No significance is to be attached to the relative vertical positions of the curves on 
the logarithmic scale. The positions were selected for arbitrary convenience in 
plotting the data. 

A statement accompanying the figures to the effect that the emission for a par- 
ticular rate study was unfiltered means that the RCA 1P28 photomultiplier tube on 
the DK-1 was exposed directly (the prism system was not traversed) to the glowing 
sample under study. This technique was possible only at the temperature of dry 
ice, for at lower temperatures nonspecific emissions interfered. The technique also 
required extra precautions to exclude outside light from the photomultiplier tube. 


CRITERION 4: ACTIVATION ENERGY 


Activation energies were measured in the usual way, i.e., from the slope obtained 
by plotting the log of the rate constant against the reciprocal of the absolute tem- 
perature. The temperatures used were those of dry ice (194° K.) and liquid 
nitrogen (77° K.). 


RESULTS 


We have obtained and studied delayed emissions for the following compounds 
(source given): bullock Jens, intact and homogenized; tryptophan and tyrosine, 
H & M Chemical Company; adenine HCl and adenosine, Schwarz Laboratories, 
Ine.; adenylic acid (muscle), adenosine diphosphate, and adenosine triphosphate, 
Pabst; desoxynucleic acid (sperm) and ribonucleic acid, Nutritional Biochemical 
Corporation; progesterone and alpha-estradiol, Ciba; quinine sulfate, Merck; 
quinidine sulfate, Fisher; acridine orange and acridine yellow, National Aniline 
Division, Allied Chemical and Dye Corporation. We will now consider the results 
for the individual materials separately. 


1. BULLOCK LENS 


For the study of the optical properties of structural proteins, the lens, being 
devoid of pigments, offers a unique material.* Lenses taken from different animals 
—bullock, rabbit, and various fishes—show a brilliant fluorescence when illuminated 
with ultraviolet light. When fresh bullock lenses were examined for fluorescence in 
the dark at room temperature, an intense blue emission was observed which had a 
lifetime shorter than 10~* second, since no light was observed in our phosphoroscope 
which had this resolution time. If the lenses were permitted to attain tempera- 
ture equilibrium at 194° K. (dry-ice temperature), an even more intense emission 
was exhibited, which persisted, with progressive fading, for about 20 seconds after 
the exciting light was removed. Such a long-lived afterglow and the shift from 
fluorescence to phosphorescence on cooling are characteristics of compounds which 
can undergo singlet-triplet intercombinations. To resolve the question whether 
the observed phosphorescence was emitted from a triplet state or by a crystal phos- 
phor, we examined the decay kinetics and the activation energy of the emission. 
Emissions of light by crystal phosphors involve activation energies and are thus 
strongly temperature-dependent, while the emission decays from triplet excitation 
states require no activation energy and are temperature-independent. In addition, 
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the first-order decay kinetics characteristic of triplet emissions is not to be expected 
with crystal phosphor emissions. 

The temperature dependence of the emission lifetime for bullock lens was first 
qualitatively studied in the following manner. Two Dewar vacuum flasks were 
placed side by side, one containing dry-ice freezing mixture, the other liquid nitro- 
gen. A test tube holding the bullock lenses was immersed in each. After tempera- 
ture equilibrium was established, the two test tubes were illuminated side by side 
for 2-3 seconds with ultraviolet light and then returned to their respective Dewar 
flasks, and the gradual disappearance of their emissions was observed. The after- 
glow faded out in both tubes at about the 
same time, although the temperature dif- 
ference between the two systems exceeded 
110° ©. This temperature-independent 
lifetime strongly pleads for the assumption 
that the observed phosphorescence comes 
from a triplet excitation state. 

These qualitative observations were 
quantitatively corroborated from kinetic 
studies carried out as described under the 
technique for Criterion 3, decay kinetics. 
Figure 1 depicts the results. It is to be 
noted first that the curves fit a first-order 
kinetic equation, i.e., the decay rate is a 
function only of the concentration of the 
residual metastable, triplet, species. 
Second, since the decay-rate constant is 
temperature-independent, the activation 
energy is negligible, and we have further 
evidence for triplet-state involvement. 

The structural protein of the lens forms sec. 0 
a fairly solid gel in vivo and can thus be 
expected to be fixed in space. This ma- i , en 
trix should also impede the diffusion of gata for henselae: te pd con ll 
dissolved protein molecules. Thisassump- _ intact ox lens, 7 = 194° K., emissions unfil- 
tion must hold even more strongly for a tered &, iatict ox leat ®, homogeael 
frozen system. The change from fluores- lens (diluted 1 to 30 with water), T = 
cence to phosphorescence with cooling can a ee mees oF oe 
thus hardly be ascribed to dimerization or to higher polymerization of protein 
molecules but rather must be attributed to a change in the solvent and its relation 
to the protein. 

If the lenses are minced in a Latapie mincer and suspended in 30 volumes of water 
(our lens homogenate), the suspension shows no striking fluorescence and no phos- 
phorescence at 194° K. The same is true for less dilute suspensions. If, however, 
the suspensions are cooled further to 77° K., they emit a phosphorescence with an 
afterglow duration indistinguishable from that of the intact lens. 

It was noted earlier that the addition of glucose to aqueous solutions of fluorescent 
substances in many instances enhances the phosphorescence of the frozen state. 
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especially at 194° K. This effect was ascribed to a more intimate contact of the 
luminescent species with water, glucose disfavoring a sharp separation of water 
crystals from the solute. This effect was observed also with the lens suspension. 
If the suspension contained 0.5 per cent glucose, freezing to 194° K., with subse- 
quent illumination, caused the appearance of a strong blue emission resembling 
that of the intact lens and persisting for about half a minute after removal of the 
illumination. This afterglow was fairly independent of glucose concentration, 
being equally strong at 0.25 per cent and at 5 per cent. The'light intensities were, 
under the conditions of our experiment, also independent, within wide limits, of the 
concentration of the suspension, the afterglow being about equally strong whether 1 
part of lens was suspended in 30 or in 300 parts of water. 

The fact that the material of the lens loses its fluorescence and phosphorescence 
by being minced and suspended in water indicates that the proteins in the native 
state of the lens are kept in a specific condition which favors these light emissions 
and the underlying excitations. This condition in all probability involves some 
specific relation to the solvent water, which relation can be restored by the addition 
of glucose. 

The phosphorescence of the frozen (194° K.) lens suspension was strongly 
quenched if, in addition to the presence of 0.5 per cent glucose, the solution was 
made to 3 per cent (v/v) with glycerol or 0.25 per cent (v/v) with methanol prior 
to freezing. Acetone was found to be inactive even at a concentration of 8 volumes 
per cent. The phosphorescence was also insensitive to the action of oxygen. 
Equally inactive were classes of ‘quenchers’ such as KI (10-* M), NaSCN (107* M), 
2,4-dinitrophenol (10-° M), and 2,4-dinitronaphthol (5 K 10-* W). 

Figure | also gives the decay curve for homogenized lens phosphorescence at 77° 
K., which was no different at 194° K. In addition, it should be noted that at 77° K. 
the decay-rate constant for homogenized lens phosphorescence was the same with- 
out as with glucose. 

The question arises whether the phosphorescence of the lens is due to some super- 
structure or to smaller units like amino acids making up the molecule. The paper of 
Debye and Edwards‘ appears to be the first systematic attempt to examine protein 
phosphorescence. They observed that most proteins emitted a brilliant blue phos- 
phorescencé at liquid-nitrogen temperature. The visible spectrum was found to be 
pH-dependent, the protein phosphorescence from an alkaline medium exhibiting a 
minor maximum at 418 my and an intense maximum at 440 my. In an acid medium 
the results were less consistent, the main maximum appearing at 417 my. The 
phosphorescent spectra of tyrosine and tryptophan were found to lie in the same 
spectral region as the protein phosphorescence. Debye and Edwards found the 
amino acid phosphorescent emission decays to be exponential and suggested that the 
exponential decay phosphorescence of proteins was due to their content of tyrosine, 
tryptophan, and possibly phenylalanine. From their finding that the phosphores- 
cent decay was temperature-independent, they felt that an electronic transition was 
involved. Grossweiner,® using flash photolysis techniques, has complemented this 
work of Debye and Edwards and has presented spectra for the metastable states of 
ovalbumin, tyrosine, phenol, tryptophan, and indole. The wave-length shift, 
temperature independence, and metastable absorption spectra obtained by these 
workers suggested triplet-state involvement. Grossweiner further concluded that 
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the metastable absorption spectrum for ovalbumin principally represented the sum 
of contributions from tyrosine and tryptophan, with negligible effects from the 
peptide linkages. 

The lens phosphorescence which we have observed (see Fig. 1) appears to have a 
slightly longer mean lifetime (5 seconds) than the protein phosphorescence studied 
by Debye and Edwards (approximately 3 seconds). In addition, we have observed 
slight differences in the kinetic data for amino acid phosphorescences when com- 
pared with the data of Debye and Edwards, but these are probably differences in 
degree only and not in kind arising from differences in technique. 


2. AMINO ACIDS 


The following amino acids were tested: dl-alanine, /-glutamic acid, /-histidine, 
dl-phenylalanine, /-tyrosine, and dl-tryptophan. None of these in aqueous solution 
showed any striking luminescence on visual observation under the ultraviolet lamp, 
frozen or unfrozen. Of these amino acids it might be anticipated from absorption 
spectral data that tryptophan and possibly tyrosine would be the only compounds 
which could possibly emit a visible fluorescence. This fact becomes apparent when 
one examines the fluorescent spectral data for the aromatic amino acids reported in 
the recent publications of Bowman, Caulfield, and Udenfriend,* Sprince, Rowley, 
and Jameson,’ and Teale and Weber.’ Phenylalanine and tyrosine emissions lie 
exclusively in the ultraviolet, while tryptophan emission, with a maximum in the 
near-ultraviolet, has a long-wave-length tail which “spills” over into the visible 
out to about 440 my. For this reason, therefore, a tryptophan solution of proper 
dilution, upon exposure to an unfiltered ultraviolet lamp, emits a distinct, though 
weak, violet fluorescence easily observable with the unaided eye. The failure to 
observe any appreciable fluorescence in a 0.01 M tryptophan solution, while a solu- 
tion of 10-4 M emits a distinct violet light, indicates that the negative results at 
higher concentration may be due to self-quenching. 

This emission picture changed markedly for the aromatic acids when the solutions 
were made to 0.5 per cent in glucose prior to freezing. Tyrosine and tryptophan in 
10-* molar concentration gave short and prolonged afterglows, respectively. These 
emissions were of considerable intensity. The behavior of tryptophan appeared 
similar to that of the lens, being blue and lasting for about half a minute. The 
wave-length shifts for these two compounds in going from fluorescence to phos- 
phorescence were visually apparent, the colorless emission of tyrosine shifting to blue 
and the violet emission of tryptophan shifting to blue. 

The marked enhancement of phosphorescence in the frozen state displayed by 
many compounds in the presence of glucose (other sugars were equally as effective) 
motivated us to examine this effect in more detail. Studies were made with trypto- 
phan as the luminous species. 

It was found that the phosphorescent decay-rate constant for tryptophan-glucose 
systems in ice was practically independent of temperature® and tryptophan concen- 
tration. In addition, we have observed an influence of glucose concentration on the 
phosphorescent emission intensity. When the tryptophan concentration was held 
constant and the glucose concentration varied, the emission intensity was observed 
to pass through a maximum, while the decay-rate constant remained essentially 
unchanged. A preliminary examination of this effect made with a tryptophan con- 
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centration of 5 X 10-4 molar resulted in a maximum phosphorescent intensity at a 
glucose concentration of 0.27 per cent. Table 1 summarizes the results for the 
slight influence of temperature and tryptophan concentration on the phosphorescent 
decay kinetics. 

Substitution of the rate data from Table 1 for the 10~* molar tryptophan system 
into the equation, 


ke AH, (= a =) 
log 


ky 2.303R 


gave an activation energy of only 76 calories. In this equation kz and k; are the 
rate constants at the high (72) and low (7) temperatures, respectively; AH, is the 
activation energy; and RF is the gas constant. 

Figure 2 illustrates the influence of temperature on the phosphorescent decay- 
rate constant for a 10~‘ molar tryptophan concentration system containing 1 per 
cent glucose. It should be observed that the tryptophan phosphorescence was in- 
sensitive to oxygen. 


TABLE 1 


THE INFLUENCE OF TEMPERATURE AND TRYPTOPHAN CONCENTRATION ON THE DeEcay- 
CoNnsTANT FOR TRYPTOPHAN PHOSPHORESCENCE* 
TRYPTOPHAN TRYPTOPHAN 

CONCENTRATION -—Rate Constant, k- CONCENTRATION —Rate Constant, k— 
(Moves PER LITER) 194° K. x eee A (MOLES PER LITER) 194° K. i pee o 
5 X 1073 0.40 0.39 10-5 0.26 0. 
10-3 0.27 0.20 10-* 0.24 0. 
10-4 0.21 0.18 


* The glucose concentration was held constant at 0.25 per cent. Rate data were obtained with the spectro- 
photometer set at 420 my and with 2.0-mm. slit widths. Systems were contained in 25 X 200-mm. Pyrex test tubes 
and, following a minimum of 15 minutes at the proper temperature to allow temperature equilibration, were illu- 
minated with the unfiltered beam of a high-pressure Hanovia mercury are for 5 seconds. 


17 
17 


Figure 3 contains the phosphorescent decay kinetic data of the tyrosine-glucose 
system. The intensity of the emitted light in the absence of glucose was too low 
to permit its measurement by the technique used in this study. At 194° K., even 
in the presence of glucose, the intensity of the emitted light was so low that it was 
difficult to obtain a satisfactory decay curve, and the reported rate constant should 
be considered accurate within +5 per cent. We have made no effort to obtain 
kinetic data for phenylalanine. 

Wave-length shifts, long-lived afterglows, first-order decay kinetics, and negligible 
temperature dependence of the rate constants may be presented as presumptive 
evidence that tyrosine and tryptophan afterglows are emitted from a metastable 
triplet state. 


3. ADENINE, ADENOSINE, ADENOSINE PHOSPHATE, ADENOSINE DIPHOSPHATE, AND 
ADENOSINE TRIPHOSPHATE (ATP) 


Euler, Brandt, and Neumiiller’® and Stimson and Reuter" presented fluorescent 
data for several purines and pyrimidines and their derivatives and gave the color 
of their emissions as a function of pH. Since we have been interested in the energy 
configuration of ATP in relation to its possible functional participation in the energy 
transformations in which the molecule participates, we have limited our study to 
the purine adenine and its derivatives. 
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the phosphorescent decay rate constant, k, Fic. 3.—Phosphorescent decay kinetic 
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= 10-4 M, glucose = 1 percent. 1, system rateconstant. 1, 10~* M tyrosine + 0.5 per 
approaching temperature equilibrium at 194° cent glue ose, 7’ = 194° K., emission unfil- 
K.; 2, system at temperature equilibrium tered; 2, 4. 4X 10-*M tyrosine + 0.3 per 
at 194° K.; 3, system at temperature equilib- cent glucose, T =77° K., emission measured 
rium at 77° K., emission measured at 420 my, at 420 mu. 




















In view of the fact that the absorption bands of the adenine compounds lie rather 
deep in the ultraviolet, it might be anticipated that only slight or no visible fluores- 
cence would be observable from them upon proper excitation. Stimson and Reuter! 
reported a faint green fluorescence for adenine sulfate in strong alkali and a faint 
blue emission in weak alkali. They observed no emission for adenylic acid solutions 
at any pH. In our work we have never observed a visible emission from adenine, 
adenosine, or any of the phosphate derivatives at any pH at room temperature. 
Adenine compounds of 0.01—-0.0005 M concentrations, when frozen to 194° K., dis- 
played a very faint blue short-lived afterglow. This emission could usually be 
enhanced in intensity by the addition of glucose to 0.5 per cent prior to freezing. 
ATP, for example, with glucose added, displayed an intense afterglow lasting for 
about 15 seconds after the exciting light was removed. ATP, in contrast to trypto- 
phan, displayed no self-quenching in 0.01 M concentration. This may have bio- 
logical significance, especially in muscle, where ATP exists in rather high concen- 
tration. While the ATP phosphorescence was insensitive to KI, NaSCN, 2,4- 
dinitrophenol, 2,4-dinitronaphthol, and acetone (8 per cent v/v) it was abolished 
by 5 per cent (v/v) glycerol and 0.1 per cent (v/v) methanol at 194° K. These later 
quenching effects were not observed, however, at 77° K. 





486 BIOCHEMISTRY: STEELE AND SZENT-GYORGYI Proc. N. A. 8. 


Figure 4 summarizes the phosphorescent decay kinetics for all the adenine sys- 
tems studied. Kinetic data at 194° K. were obtained for ATP only and indicated 
that the rate constant was practically in- 
dependent of temperature. Though the 
emission intensities were found to be mark- 
edly sensitive to pH, the rate constants 
remained essentially independent of this 
variable. An important feature observed 
in studies of pH effects on the emission in- 
tensities was the finding that the amino 
nitrogen at position 6 on the adenine ring 
must be dissociated for luminescence to 
occur. Varying the pH of adenine-com- 
pound systems from 2 to7 by units of 1 gave 
a progressive increase in the phosphores- 
cent intensity from nil to brilliant. Since, 
according to the work of Alberty, Smith, 
and Bock,'? the amino group has a pK, of 
approximately 4, we visualize the increased 
alkalinity as dissociating the ‘adenin- 
ium” cation. This dependence of luminous 
intensity on the degree of dissociation 
of an amino group is not a unique find- 
ing. Ley and von Engelhardt!’ noted that, 
whereas the aniline molecule was fluores- 
, cent, the anilinium cation was nonfluores- 
pasted . sed cent. By contrast, Ley and Grife! and 

Fic. 4.—Phosphorescent decay kinetic Dickson® found that the blue fluores- 
data for, adnine (1, adenosine - (4) cence of alpha and. beta-naphthylamine 
phosphate (4), and adenosine-5’-triphos- was not lost by proton addition but, 
Date (6), i, somenurations COOL AM instead, was shifted into the ultraviolet, 
to freezing. 7 = 77° K. All emissions with, however, a considerable diminution 
measured at 420 mu. k = rate constant. in intensity. Lewis and Kasha'* have 

given phosphorescent emission data for ani- 
line and alpha- and beta-naphthylamine, but we have found no reports on the phos- 
phorescence of the anilinium cation or of the alpha- or beta-naphthylammonium 
cations. 

For the adenine derivatives studied, therefore, the wave-length shifts, i.e., from 
no visible emission to an intense blue afterglow, the long lifetimes, and the first- 
order kinetics plead for triplet-state participation in the emissions observed for 
these compounds. 





log Emission intensity 
ron 











4. NUCLEIC ACIDS 


Nucleic acids show a long-lived phosphorescence which might be due to the col- 
lective action of the column of cyclic structures constituting, according to the model 
of Watson and Crick,” the core of the molecule. The central cyclic compounds 
would thus, according to this assumption, be analogous to Scheibe’s pseudoiso- 
cyanine dyes,'* which polymerize to columns in the shape of ‘‘money-rolls” through 
which excitons are freely propagated by resonance. 
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The phosphorescent light emission from aqueous solutions of DNA and RNA 
preparations, 0.1—-1.0 per cent, frozen to 194° K., could be prolonged by the ad- 
dition of 0.254 per cent glucose. This prolongation was abolished by small con- 
centrations of glycerol (3 per cent, v/v) or methanol (0.25 per cent, v/v), while 
acetone (16 per cent, v/v) was found to be inactive. The phosphorescence was 
found to be insensitive to KI, 2,4-dinitrophenol, 2,4-dinitronaphthol, salicylate 
(0.002 M), or pyrrolazote (2 X 10-4 M). Morphia and aconitine (10-* M) were 
equally inactive, while colchicine in 2 X 10-4 M concentration quenched DNA 
strongly, the effect increasing somewhat on storage. Since the nucleic acids were 
present in rather high concentration, 0.5~-1 per cent, there was a disproportionality 
between the molarity of the alkaloid and the mass of the nucleic acids or the number 
of nucleotides building the nucleic acid molecule. This pleads for a collective 
action of the nucleotides in the production of phosphorescence. Scheibe!* found 
the energy transmission in his dye polymers decreased by relatively small concen- 
trations of hydroquinone, the conduction within the “‘money-rolls” being cut by 
the interposition of a “dud,” a hydroquinone molecule. Hydroquinone at 
0.001 M also strongly quenched the phosphorescence of DNA. The first pre- 
liminary experiments repeated with high-polymer DNA indicated that the quench- 
ing action of colchicine was about ten times stronger with this material than with 
the low-polymer DNA, the alkaloid being active in 10-* M concentration. This 
concentration is near the level which is effective in vivo. If the alkaloid acts by 
cutting the energy transmission through the row of central bases, then it should be 
more active the longer the molecule. 

The carcinogen benzpyrene also showed a fairly strong quenching activity. In 
these experiments the carcinogen was dissolved in acetone (0.1 per cent). One- 
tenth of a milliliter of this solution added to 5 ml. of 0.5 per cent DNA almost 
completely eliminated the long afterglow enhanced by glucose. In this instance, 
too, the action seemed to be more marked when high-polymer DNA was used in the 
experiment. Benzpyrene, as is generally known from the work of Brock, Druckrey, 
and Hamperl,?° emits a blue light when molecularly dispersed and a yellow-green 
light in colloidal emulsion when illuminated by ultraviolet light. This can be 
strikingly demonstrated by adding a drop of acetone solution of benzpyrene to 
water. The resulting colloidal solution emits a yellow-green fluorescence which 
changes dramatically to blue if the solution is shaken up with a little solid caffein. 
The purine forms a complex with the carcinogen, bringing it into molecular dis- 
persion. Analogous behavior can be demonstrated with DNA, indicating that the 
carcinogen complexes with the bases of the nucleic acids. This might have been 
anticipated also from the work of Booth and Boyland,®?! who demonstrated that 
the variably carcinogenic bicyclic dibenzacridines formed well-defined crystalline 
complexes with purines and nucleic acids. 

Figure 5 depicts the phosphorescent decay kinetic results as measured for DNA 
and RNA at 77° K. Decay kinetics are not reported for these compounds at 
194°K., for the emission intensities were too low to be recorded by the technique 
here employed. The rate constants given should be considered only as indicating 
the order of magnitude of the decay rate, for the purity of these preparations was 
uncertain. The rate constant for DNA appears to be complex, and until highly 
purified preparations of reasonably well-known polymer composition are avail- 
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Fic. 5.—Phosphorescent decay kinetic emission measured at 510 mu, 7’ = 77° K.; 
data for desoxy ribonue leic acid and ribonu- 2, progesterone, 7.3 X 10~* M in 95 per cent 
cleic acid. J’ = 77° K. k = rate constant. ethanol + 0.5 per cent glucose, blue emis- 
1, 1 per cent DNA + 0.25 per cent glucose, sion measured at 420 mu, T = 77° K.; 3, 
480 uw; 2, 1 per cent DNA, no glucose, 520 alpha-estradiol, 10~* M in 95 per cent eth- 
mu; 3, 1 per cent RNA + 0.25 per cent glu- anol, no glucose, violet emission measured 
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cose, 400 mu. at 410 my, 7’ = 77° K. k = rate constant. 














able, we will restrict our comments to one observation, namely, that while a 0.1 per 
cent DNA solution (glucose present) displayed a blue phosphorescence, a 1 per cent 
DNA solution (also in the presence of glucose) gave a green phosphorescence. We 
have noticed a somewhat analogous phenomenon with adenosine. Whereas ade- 
nosine in the absence of glucose emitted a distinct green phosphorescence, in the 
presence of glucose the emission reverted to the blue characteristic for adenine and 
the adenosine phosphate derivatives. The effect was most marked at low pH val- 
ues, and none of the other adenine compounds showed the effect. Therefore, the 
difference observed in the rate constant for DNA may reflect a mixture of these two 
emissions. The data predict that the green emitting system, if real, will be found 
to have the higher rate constant. No such wave-length shifts have been observed 
for RNA. 
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First-order kinetics and increased lifetimes are the only criteria we have obtained 
so far for the nucleic acids which suggest triplet-state involvement. 


5. MISCELLANEOUS 


Figure 6 contains rate data for quinidine sulfate, progesterone, and alpha- 
estradiol phosphorescences, all measured at 77° K. and at the wave lengths given. 
Quinine sulfate gave practically the same rate constant as quinidine sulfate, and 
glucose addition was necessary for both compounds in order to obtain sufficient 
luminous intensity to permit decay curves to be run. The emission was a distinct 
green in both instances and of rather intense brightness, though of very short du- 
ration. The blue emission of progesterone and the violet emission of alpha-estradiol, 
though of low intensities, could both be measured without glucose addition; in fact, 
glucose addition had little influence on their emission intensities. It should be 
noted that the solvent for the hormones was 95 per cent ethanol, not water. The 
very low emission intensities at 194° prevented rate-constant measurements by the 
technique used in this study. 

Quinine and quinidine sulfate fluoresced blue, while the phosphorescence was 
green. We have measured the fluorescence of alpha-estradiol (unpublished) and 
found it to be exclusively ultraviolet, with a maximum at approximately 310 u. 
We were unable to obtain any indication of progesterone fluorescence. With these 
compounds, therefore, wave-length shift (with the possible exception of proges- 
terone), first-order kinetics, and increased lifetimes represent the criteria suggestive 
of triplet-state involvement. 


6. DYES 


In Figure 7 are presented a few preliminary kinetic studies made for the dyes 
acridine yellow and acridine orange. The break in rate curves at the higher dye 
concentration and elevated water content suggest that polymer involvement may 
influence the rate constant. Zanker®? has made a detailed study on the influence of 
solvent, concentration, temperature, and pH on the absorption and emission char- 
acteristics of acridine orange that may prove to be pertinent here. 

We have measured the following emission shifts for these two dyes in going from 
fluorescence to the phosphorescent state. Acridine yellow fluoresced at 511 my and 
phosphoresced with a maximum at 614 my. Acridine orange gave a fluorescence 
band with a maximum at 530 my and a phosphorescence-band maximum at 660 mg. 
These wave-length shifts, first-order kinetics, and long-lived afterglow, we cite, 
therefore, as evidence for triplet-state participation in the phosphorescence of these 
dyes. 

7. WATER 

A finding which we can mention only briefly at this time was an observation con- 
cerning the emission from ice following ultraviolet excitation.2*> The ‘‘afterglow” 
was not visible in an undarkened room, and we have made no attempt to visualize 
it with dark-adapted eyes. The intensity appeared greatest when the emission 
from ice at 77° K. was permitted to impinge unfiltered directly onto the RCA 
1P28 photomultiplier tube set electronically at maximum sensitivity. Precau- 
tions were always taken to insure that this ice emission did not interfere with the 
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dine orange. 7’ = 77° K. k = rate con- 
stant. 1, 10-4 M acridine yellow in 90 per Fig. 8.—Phosphorescent decay kinetics 
cent ethanol + 10 per cent water (v/v), emis- for the unfiltered afterglow from ice at 77° 
sion measured at 550 mu; 2-3, 5 K 10~-* M and 194° K. The ice was frozen around the 
acridine yellow in 50 per cent ethanol + 50 outside of a hollow aluminum bar (freezing 
per cent water (v/v), emission measured at mixtures inside) and exposed for 5 seconds 
530 mu; 4-6,5 X 10~* M acridine orange in to the unfiltered beam of a high-pressure 
50 per cent ethanol + 50 per cent water mercury lamp prior to the kinetic run. k= 
(v/v), emission measured at 600 mu. rate constant. 1 = 77° K., 2 = 194° K. 
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kinetic measurements. Rate data for this emission at 77° and 194° K. are given 
in Figure 8. These results were fairly reproducible. It is interesting that the de- 
cay appears to obey first-order kinetics and has a rather low activation energy. 
Unfortunately, the intensities were too low to permit us to obtain reliable informa- 
tion about their spectral distribution. The rate constant, however, appeared to 
be about the same whether the emission was permitted to traverse a Corning 5860 
band-pass filter (maxinum wave length passed at 360 my) or a Corning 3486 cutoff 
filter, cutting out light below 510 mu. 


DISCUSSION AND SUMMARY 


Fluorescent and phosphorescent spectral data, kinetic studies (including order of 
reaction), lifetimes, and activation-energy considerations have been presented for 
a wide variety of biological compounds and the dyes acridine yellow and acridine 
orange. Results were obtained which we submit as evidence that electronic transi- 
tions involving singlet-triplet intercombinations are functional in these compounds. 

The marked optical changes observed in freezing aqueous solutions, namely, the 
shift from fluorescence to phosphorescence, appear to be a universal property of 
fluorescent materials. 
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REPRESSED. AND INDUCED ENZYME FORMATION: 
A UNIFIED HY POTHESIS* 
By Henry J. VoGe. 
DEPARTMENT OF MICROBIOLOGY, YALE UNIVERSITY 


Communicated by E. L. Tatum, April 2, 1957 


A specific antagonism, termed “enzyme repression,” of the formation of the 
biosynthetic enzyme acetylornithinase in Escherichia coli has recently been de- 
scribed.! The repressing agent was found to be arginine (or a substance related to 
arginine), which is the “end product” of the biosynthetic sequence involved. 
Enzyme repression by an exogenously available end product (combined with 
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anabolic utilization of this product) represents a highly effective regulatory device, 
since it permits the cell involved to curtail superfluous enzyme production. Enzyme 
repression by an endogenously produced substance can provide a feed-back ar- 
rangement that governs, in line with cellular demands, the flow of metabolites 
through the pathway leading to this substance. Accordingly, enzyme repression 
is a control mechanism that complements enzyme induction: in either case, the 
cell tends to form enzymes when they are needed and tends not to form enzymes 
when they are not needed. It has been pointed out that repressibility appears to 
be a widespread property of enzyme-forming systems, which presumably was posi- 
tively selected in the course of evolution. The possibility that enzyme repression 
plays a role in the development or differentiation of higher forms of life has also 
been briefly mentioned. The present communication is concerned with a hypothesis 
of enzyme induction and repression. 

““Small-M olecule”’ Control of Enzyme Formation.—As previously discussed, repres- 
sion and induction appear to be separable from the process of enzyme formation 
per se; in particular, there are good grounds for the view that induction is not a 
sine qua non of enzyme formation. It is for this reason that repression and induc- 
tion are best pictured as control mechanisms rather than as inherent features of 
enzyme biogenesis. 

The following small-molecule control situations with respect to enzyme forma- 
tion are conceivable: (i) indifferent enzyme synthesis (no control); (ii) induced 
enzyme synthesis; and (iii) repressed enzyme synthesis. One and the same enzyme- 
forming system may be both inducible and repressible,? or may be neither. Again, 
a given enzyme-forming system may be inducible but not repressible, or repressible 
but not inducible. 

Repression and induction can thus be considered to have a functional corre- 
spondence. In addition, these processes may well depend on corresponding molecu- 
lar mechanisms. The question then arises how a relatively small molecule, such as 
an inducer or a represser, can specifically modify the formation of a particular 
protein of large molecular weight. 

A Unified Hypothesis of Induction and Repression.—F rom their relative specificity 
it seems reasonable to assume that induction or repression depend, respectively, on 
an interaction between inducer or represser (or their ‘active’ derivatives), on the 
one hand, and a macromolecule, presumably of template nature, on the other. 
Such interaction either may involve the template and the “regulator” (inducer or 
represser) or may involve template, template product, and regulator. In line with 
current assumptions,’ this template may be one immediately involved in the pro- 
duction of enzyme protein. 

Since action of the regulator—in particular, of the inducer—does not appear to 
be a sine qua non of enzyme formation,' the cell is thought to have, prior to contact 
with the inducer, the essential information required for enzyme structure. Accord- 
ingly, the inducer would seem to act not by furnishing a prototype for the configura- 
tion of the enzyme molecule but, rather, by improving the performance of the 
template. Thus the inducer would assist the template catalyst, and may be con- 
sidered a “promoter.’’ In contrast, a represser may be compared to a catalyst 
“poison.” 

It is proposed that inducers and repressers act by affecting the rate of dissociation 
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of a template product from its template. Repression could then be the result of 
binding of newly formed enzyme protein to its site of synthesis through the agency 
of the represser involved. Induction may reflect the neutralization of a binding 
effect which, in the absence of the inducer, would tend to hold the nascent protein 
near its template. It is considered that rapid removal of template product from 
template will accelerate enzyme formation, while occupation of a template by its 
product will prevent the template from functioning in further enzyme synthesis. 
It is not necessary to assume that the represser causes binding of enzyme to tem- 
plate by acting as a “bridge”; it is possible that the represser causes a change in 
macromolecular configuration which results in binding. Similarly, an inducer may 
act either through the local “‘neutralization”’ of a bridge group or through an effect 
on macromolecule configuration. 

Regulators may be expected to show greater or lesser tendencies to depart from 
their locales of action. Thus, in line with Pollock’s observations,‘ the inducer of 
penicillinase (in contrast to inducers of 6-galactosidase) appears to be tightly held 
at the induction site. 

In the present view of regulator action, the separation of template product from 
template is regarded as a more or less well-defined, single event. However, this 
view is not necessarily in conflict with a (simultaneous) multiple partial functioning 
of a template, such as is envisaged by Dalgliesh.° In ‘‘two-dimensional” protein 
synthesis,® however, one wonders to what extent such partial functioning of ‘one- 
dimensional” templates is likely to occur. 

The hypothesis here proposed, which will be referred to as the ‘‘regulator hy- 
pothesis,” implies that all enzymes, be they constitutive or adaptive, inducible or 
repressible, are synthesized in the same general manner but may differ from case 
to case in susceptibility to regulators. This picture differs from Cohn and Monod’s’ 
“generalized induction hypothesis” (ef. Vogel’). 

The Regulator Hypothesis and the “Final Shaping” of Enzymes.—In its general 
form, the regulator hypothesis is concerned with the dissociation of a template 
product from its template. Let it now be assumed that the template product is 
the nascent enzyme protein and that the regulator acts in a region that corresponds 
to the dynamic site of the “finally shaped’’* (see below) enzyme molecule. If so, 
the regulator would be expected to be in close physical contact with the nascent 
dynamic site at the very moment that the template product separates from the 
template. Such contact may well be of considerable consequence for the fruitful 
final shaping of the enzyme molecule, particularly in the folding of the nascent 
enzyme protein or in the aggregation of template products. A directed folding 
process in protein biogenesis has been considered by previous authors, especially in 
connection with antibody formation (Pauling-Haurowitz-Mudd hypothesis; ef. 
Haurowitz*). In principle, the importance of the final shaping process has been 
anticipated—for example, in the specific-precursor postulates of Monod® and of 
Pollock.* 

Lipmann” has recently drawn attention to the final shaping of proteins as dis- 
tinct from the emergence of a defined sequential arrangement (‘‘patternization’’) of 
component amino acids. While, as suggested above, final shaping may be affected 
by regulators, this process is considered to be capable of occurring in the absence 
of regulators. For instance, appropriate folding of nascent enzyme protein could 
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Template + Activated patternized precursors be a spontaneous consequence 
of the patternization  in- 
volved.” The possible role of 
regulators in influencing (a) 
the separation of template 
Sens product from template and (b) 
| the final shaping is presented 

|—Active regulator . pects 

schematically in Figure 1." 
Final shaping] Comments on the Regulator 
res tatin 
Hypothesis—The regulator 
Template + Enzyme hypothesis is generally con- 
Fira. 1.—Scheme of regulator action sistent with and supported by 
the available data on enzyme 
induction and enzyme repression. This hypothesis readily accounts for the ob- 
servation (cf. Pollock*) that inducible enzymes usually show a “‘basal’’ enzyme level 
in the absence of added inducers. The basal level would correspond to the rate of 
enzyme production that the cell is capable of achieving without the benefit of the 
promoting effect of the inducer. 

The emergence of mutants in which an originally inducible enzyme has become 
constitutive (cf. Monod and Cohn'*) can be visualized, at least in some cases, as 
involving a template alteration such that the corresponding template product is 
enabled to dissociate rapidly enough from its site of formation without aid from a 
regulator substance; in other words, there would be a change from inducible to 
indifferent enzyme synthesis. 

The present hypothesis is also in harmony with other results, including the find- 
ing'* that there is no necessary relationship between the properties of a given sub- 
stance as an inducer (or a represser), on the one hand, and a substrate or a com- 
plexant, on the other. Again, the regulator hypothesis is in accord with the further 
findings that different inducers can evoke enzymes of indistinguishable specificity !* 
and that an enzyme produced under conditions of repression has been found not to 
differ detectably from the corresponding nonrepressed enzyme. ! 

A dual regulator function, such. as that illustrated in Figure 1, appears to be 
consistent with available results on induced enzyme formation, including those of 
Cohn and Torriani.'* These authors have considered the possibility that the 
8-galactosidase and the immunochemically related (but enzymatically inactive) Pz 
protein of /. coli are synthesized independently at the same site. In line with this 
possibility, the increased rate of 8-galactosidase formation and the decreased rate 
of Pz formation observed upon induction may well involve the above-mentioned 
dual regulator effect. However, since it has not been definitely established that 
8-galactosidase and Pz actually are produced at the same site, other explanations 
for the relation between these two proteins remain to be considered. 

The Regulator Hypothesis and Other Theories of Enzyme Formation.—In a recent 
review on theories of enzyme adaptation, Mandlestam™ has aptly commented on 
the unusually high ratio of theories to facts in this field; he has also pointed out, 
however, that it is unlikely that this subject would have advanced so rapidly if the 
theories had not been formulated. A comparison of the regulator hypothesis with 
some earlier concepts might therefore not be without some value. 


Template — Nascent enzyme complex 


Template + Nascent enzyme 
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The major theories of induction can be classified as (1) those that depend on some 
kind of combination of inducer with (finally shaped) enzyme (for either equilibrium 
displacement or enzyme stabilization or template-like action) and (II) those that 
do not. 

Group I includes Yudkin’s'* mass-action theory and Mandlestam’s” extended- 
mass-action theory; Spiegelman’s'® plasmagene theory and Monod’s® specific- 
precursor theory; and Campbell and Spiegelman’s template-inducer-enzyme com- 
plex theory (cited by Spiegelman and Campbell’), in which such a triple complex 
is considered to be the only structure that can effectively and rapidly lead to enzyme 
synthesis. Group II contains Pollock’s* organizer hypothesis (which has been 
extended by Cohn and Monod’), Monod and Cohn’s!? inducer-prototype hypothe- 
sis, and the present regulator hypothesis. 

Arguments against those theories in Group I which invoke a binary complex of 
(finally shaped) enzyme and inducer have been presented by Monod et al.'* These 
authors showed that, for 6-galactosidase of FZ. coli, inducer function is separable 
from substrate or complexant functions. Their arguments do not apply to a triple 
complex, such as is contemplated by Campbell and Spiegelman. The triple-com- 
plex hypothesis was advanced primarily to account for “long-term adaptation” in 
yeast.* Since, however, other possible explanations'® of long-term adaptation have 
not been excluded, one is not forced to assume that a given finally shaped enzyme 
participates in the synthesis of its own molecular species. It is also worthy of note 
that the above-mentioned results of Monod et al.'* do not provide evidence against 
a complex involving an inducer and the corresponding nascent enzyme (which may 
be similar to, but not identical with, the finally shaped enzyme in affinity for com- 
plexants and in enzymatic activity). 

Group II, which includes hypotheses that do not invoke an active participation 
of finally shaped enzyme in the production of its species, can be subdivided into 
the organizer and inducer-prototype hypotheses, on the one hand, and the regulator 
hypothesis, on the other. The inducer-prototype hypothesis assumes that the steric 
configuration of an enzyme’s active site is directly derived from the inducer.’ *° 
Similarly, the organizer hypothesis has been interpreted’ as implying that the 
characteristic pattern of an enzyme’s dynamic site is imposed by a co-organizer 
(“active”’ inducer). Arguments against such a prototype or co-organizer function 
of the inducer have recently been discussed.! In the regulator hypothesis, the 
assumption of prototype function is avoided. 

Summary.—A hypothesis (regulator hypothesis) has been proposed to account 
for induced and repressed enzyme formation, which are viewed as analogous 
phenomena in terms of function and of mechanism. Regulators (inducers or re- 
pressers) are thought to act, respectively, by promoting or impeding the separation 
of template product (nascent enzyme protein) from template. An effect of regu- 
lators on the final shaping of enzyme molecules is also contemplated. Regulator 
action is pictured as a control mechanism rather than as a sine qua non of enzyme 
biogenesis. 


* This work was aided by a contract between the Office of Naval Research, Department of the 
Navy, and Yale University, and by the Damon Runyon Memorial Fund. 
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THE SITES OF ORTHOPHOSPHATE UPTAKE BY BARLEY ROOTS 
By C. E. Hagen, J. E. Leacert, anp P. C. Jackson 


AGRICULTURAL RESEARCH SERVICE, UNITED STATES DEPARTMENT OF AGRICULTURE, BELTSVILLE, 
MARYLAND 


Communicated by Sterling Hendricks, April 3, 1957 


Kinetic analysis of measurements of steady-state uptake has established the ex- 
tent to which the HPO] and H,PO;j ions contribute to the total uptake of ortho- 
phosphate by excised barley roots.'. The binding compounds involved in the rate- 
limiting step of orthophosphate uptake are shown here to be components of the 
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respiratory chain. Further, the point where each phosphate species is coupled to 
the electron-transfer system is shown to be the site of energy conservation during 
oxidative phosphorylation by mammalian mitochondria.2~* The bearing of these 
results on theories of ion uptake is not developed. Discussion rather awaits results 
on the pathways of uptake of ions other than phosphate. 


KINETICS OF PHOSPHATE UPTAKE 
Uptake of ions by plant roots arises from an initial combination between an ion 
and a specific binding compound, as written in equation (1). Breakdown of this 
labile ion complex results in the uptake of ions into the system, as expressed by equa- 
tion (2). 


ky 
R + M cateide c————- MR, 
ke 


ke 


MR Cees R’ + M insiae, 
ke 


where F represents a binding site or compound; M, the ion; MR, the labile inter- 
mediate; k, the rate constant for each reaction. The rate constants of association 
and dissociation, k, and ks, of the intermediate, MR, are large with respect to the 
rate constant, k3, of the breakdown of the intermediate; k, is negligible. 

Uptake of the ion into the root may be described by the velocity, v, in the linear 
equation (3), when the concentration of the external ion is held constant and the 
concentration of R is constant during the period of absorption. Equation (3) is 
analogous to the steady-state expression describing the kinetics of enzyme reac- 
tions.® 6 


u “ , 
U ' [M] Ke + J max (3) 
The maximum uptake at infinite ion concentration is Vax, and K,, denotes the 
apparent dissociation constant of the intermediate ion complex, MR. When the 
steady-state uptake, v, by excised barley roots from single-salt solutions of alkali 
cations,’ as well as halides,® is plotted against v/[M], a straight line results over a 
wide range of concentrations. Thus a single first-order reaction is involved in the 
rate-limiting step of uptake for each of these monovalent ions, with the maximum 
uptake, Vinax, given directly by the ordinate intercept and the apparent K,, by the 
slope, according to equation (3). The steady-state uptake measurements of ortho- 
phosphate, when plotted in this manner, are curvilinear, because two first-order re- 
actions are acting simultaneously, but independently on different substrates which 

are in a non-rate-limiting equilibrium. 

Steady-state analysis of the measurements of orthophosphate uptake from solu- 
tions containing various concentrations of orthophosphate (10-*-10~-* M) and hy- 
drogen ion (pH 4-8) leads to the quantitative interpretation that excised barley 
roots absorb the ion species H,PO{ through one site and HPO through another.' 
The system is described by the following reactions: 


H;PO, 2 H+ + H.POZ @ H+ + HPO] @ H+ + POT, (4) 
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Ri + H:POZ @ R,H2PO,, (5) 
Ry} + HPO; @ R,HPO;. (6) 


Further, the hydroxyl ion, OH~, competitively inhibits the uptake of both ion 
species without affecting the concentration of either of the binding sites: 


R* + OH- @ R,OH, (7) 
Ri + OH- = R,OH. (8) 


The concentration of all the intermediate ion complexes, expressed in equations 
(5)—-(8), approximate equilibrium values, since their rate constants of association 
and discussion are large in comparison with the rate constant of breakdown, k;. 
Thus the observed velocity of steady-state uptake, v, is directly proportional to the 
concentration of the ion-complex intermediates at all values of orthophosphate con- 
centration. However, the total phosphate associated with the root, as measured in 
uptake studies, also includes the equilibrium value of the phosphate in combination 
with both sites. The observed uptake where these values contribute a significant 
amount to the total is denoted as v’, as expressed in equation (9): 


v= { [RaH2PO4 |ka + [RoHPO 5 Jhsn } t + [R,H2PO; ] + [R,HPO § }. (9) 


Sorption periods measured in seconds show a significant amount of phosphate in 
association with the roots because of the equilibrium value of the intermediate 
phosphate-site complexes (Fig. 1). Sixty-second absorption measurements must, 
therefore, be expressed by equation (9), where a plot of v’ against v’/(ZP] (Fig. 2) 
results in a numerical value of the ordinate intercept which describes Vimgx ¢ and 
Vmax » for each reaction. The numerical values of [R,H2PO,] and [R,HPOj{] are 
equal to [R,] and [R,] at infinite substrate concentration (v’/[ZP] = 0), since the 
sites exist completely in the combined form. Therefore, the maximum uptake for 
l-minute absorption periods is expressed as follows: 


View a + oe eg { [Ra lksa + [Ro kao} t + [Ra] + [R,]. (10) 


Inspection of equation (9) shows that the phosphate in combination with the sites 
becomes negligible with respect to the total phosphate absorbed with longer periods 
of uptake, ¢. The 3-hour sorption periods used in previous studies! were of suf- 
ficient duration to express uptake solely as a velocity. 


RELATIONSHIP OF PHOSPHATE UPTAKE TO RESPIRATION 


The hypothesis for electron transport in the oxidative phosphorylation system 
has been extended by the accurate measurement and rigorous interpretations of the 
effects of substrates and inhibitors on the oxidation-reduction levels of the compo- 
nents of the respiratory chain.* !° Certain inhibitors intercept electron transport 
at different points along the respiratory chain. These inhibitors cause those com- 
ponents on the oxygen side of the interaction to become more oxidized, while those 
on the substrate side become more reduced. Other types of inhibitors restrict or 
block electron flow and cause a difference in the steady-state oxidation-reduction 
level of the components of the respiratory chain. Each type of inhibitor has been 
studied with respect to the point of interruption of electron transport,’ as expressed 
by the following scheme: 
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The sites at which energy-conservation reactions occur in isolated mammalian 
liver mitochondria?~‘ and in intact baker’s yeast* * have been shown to be reduced 
cytochrome c, reduced cytochrome b, and DPNH. 

A kinetic interpretation of the effect of similar respiratory chain inhibitors on the 
absorption of orthophosphate shows that reduced cytochrome b and DPNH are 
the components of the respiratory chain involved in the rate-limiting steps of ortho- 
phosphate uptake by excised barley roots. 

The method for determining the components of the respiratory chain taking part 
in the rate-limiting steps of orthophosphate uptake is based on the following ap- 
plications of first-order reaction kinetics: 

a) A separation of the component reactions involved in both the presence and 
absence of an inhibitor, according to Figure 2, determines whether the effect is on 
one or both of the first-order reactions involved in orthophosphate uptake. 

b) A change in slope of either reaction in the presence of an inhibitor (Fig. 3) 
reflects an apparent change in K,,. A strictly competitive type of inhibitor is de- 
noted by a change in K,, without change in Fic: 

c) An inhibitor which changes the intercept with the ordinate, Vj.ax, of either 
reaction reflects an effect on either the respective concentration of the site, [R], 
k3, or both, according to equation (10). 

d) Measurements of uptake with time permit a further separation of the effect 
of an inhibitor which changes the Vox of either reaction. When absorption is 
plotted against time in the presence of an inhibitor, a change in the slope of the re- 
sulting curve, with no change in the ordinate intercept, reflects an effect solely on 
ks, according to equation (10). An effect on [R] is reflected by both the slope and 
the intercept changing by the same ratio in the presence of an inhibitor. 


EXPERIMENTAL METHODS 


Excised roots from 6-day-old seedlings of Hordeum vulgare, variety Atlas 46, were 
used as the experimental material. The seedlings were grown essentially as de- 
scribed by Epstein and Hagen.’ The seeds were soaked for 24 hours in continu- 
ously aerated demineralized water and then planted on a layer of cheesecloth sup- 
ported by a stainless-steel screen suspended over a solution of 2 X 10-* M CaSQ,. 
Roots were excised from the seedlings, which were grown for 5 days in a dark 
chamber maintained at 24° C. 
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Fics. 1-4.—The uptake of orthophosphate by 0.5 gm. of excised barley roots. Fig. 1, from 
1 X 10-* and 5 X 10~* M orthophosphate as a function of time. Fig. 2, as a function of ortho- 
phosphate uptake divided by total orthophosphate concentration with ‘separation of the com- 
ponent reactions, in the presence and absence of nembutal. Fig. 3, as a function of orthophosphate 
uptake divided by total orthophosphate concentration with separation of the component reac- 
tions, in the presence and absence of pyrophosphate. Fig. 4, from 2 X 10~¢ M orthophosphate in 


the presence and absence of azide. 


The concentration of orthophosphate remained essentially constant over an ab- 
sorption period of 5 minutes, using 25-ml. quantities of solution with each 0.5-gm. 
portion of excised roots. The roots, which were placed in test tubes, were rinsed 
three times with demineralized water. Single-salt solutions, at pH 4.0, of ortho- 
phosphate containing P**, both in the absence and in the presence of various con- 
centrations of substrate or inhibitor, were added to the roots at zero time. The 
pH value remained constant to within 0.1 unit during the absorption periods. At 
pH 4, the substrates and inhibitors were predominantly in the undissociated form, 
which has been shown to penetrate cells more readily than the charged ions.!!~"* 
At the end of the absorption period, the roots were rinsed four times with deminer- 
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alized water, which was sufficient to remove the phosphate not associated with ac- 
tive uptake. The absorption of phosphate was determined from the P*? associated 
with the roots.! 


EXPERIMENTAL RESULTS 


Location of Phosphate-binding Sites —Amytal intercepts electron transport be- 
tween DPNH and flavoprotein.® The presence of this inhibitor causes all the com- 
ponents in the respiratory chain on the substrate side of flavoprotein to become re- 
duced, while those on the oxygen side become oxidized. A decrease of orthophos- 
phate uptake by barley roots was observed on testing with either amytal or nem- 
butal. Separation of the component reactions showed a decrease solely of Vnax » 
(Fig. 2). Both the slope and the intercept decreased by the same ratio when the 
phosphate uptake was measured as a function of time; therefore, the decreased up- 
take of the HPO{ ion was due to a decrease solely of [R,]. The action of amytal 
and nembutal establishes that the rate-limiting step of HPO] uptake could be 
either an oxidized component of the respiratory chain on the substrate side of 
flavoprotein or a reduced component on the oxygen side. 

Paraphenylenediamine has been demonstrated to be a reductant of cytochrome ec. 
All components in the respiratory chain toward the substrate side of cytochrome c 
become reduced in its presence. An increase in uptake occurred when the effect of 
this inhibitor was tested on orthophosphate uptake. Similar results were obtained 
with naphthoquinones, methylene blue, and 2,3-dimercaptopropanol. Separation 
of the component reactions showed an increase in Vax », With no effect on Vinax a 
Measurement of the phosphate uptake as a function of time established that [R,] 
had been increased, since both the intercept and the slope increased by the same 
ratio. The deduction is that the rate-limiting step of HPO{ uptake is coupled to 
reduced cytochrome b, since only changes in [R,] are involved both in the increase 
of HPO{ uptake in the presence of this type of inhibitor and in the decrease of 
HPOj{ uptake in the presence of nembutal and amytal. Further, the uptake of 
H.PO{ is coupled to a reduced component of the electron-transfer system toward 
the substrate side of flavoprotein. 

Addition of acetoacetate, the product of 6-hydroxybutyrate oxidation, inhibited 
both HPO{ and H.PO{ uptake. The inhibition was a function of a decrease in 
Vinax » 20d Vinay a The inhibition was further resolved by time studies to show that 
a decrease solely in [R,] and [R,] had occurred. Thus H,PO;{ uptake is coupled to 
a reduced component of the respiratory chain which lies between the substrate and 
flavoprotein: the DPNH. 

Systems in the steady state of aerobic metabolism generally maintain cytochrome 
b in a considerably less reduced state than DPNH, which is largely reduced. A 
measurable effect on HPO 4 uptake could occur with a negligible change in H,PO; 
uptake in the presence of inhibitors which cause both cytochrome b and DPNH to 
become more reduced. This is the case, apparently, in the nonmeasurable effect of 
p-phenylenediamine, nembutal, and amytal on the DPNH which is coupled to 
H.PO,4. 

Site Concentrations —Chance and Williams have shown that the concentration 
of DPNH exceeds many fold that of cytochrome b in intact mammalian mitochon- 
dria. The relative concentrations of the two sites involved in the rate-limiting 
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step of orthophosphate uptake by barley roots may be approximated without ref- 
erence to the chemical identity of the sites. The ratio of the concentrations pro- 
vides a criterion for corroboration of cytochrome b and DPNH as the two sites of 
uptake. Most of the total phosphate associated with the roots from solutions of 
1 X 10-* M orthophosphate (where v’,,,/[ZP] = 5 X 10~‘) is contributed by the 
HPO; ion (Fig. 2). Thus an approximation of [R,HPOj ] can be determined from 
the value of the intercept at zero time obtained from extrapolation of uptake 
measurements with time (Fig. 1). A value of ky can be calculated by substitution 
of a value of v,’ and [R,HPOZ] at 1 X 10-* M[ZP], obtained from Figure 1, into 
the expression for HPO4 uptake, equation (12): 


», = [RHPOG |knt + [R,HPO; I. (12) 


Solution of the equation gave a value for ky, of 2.3 K 10~? moles P/mole R, X sec. 
The steady-state uptake of HPO; is directly proportional to the concentration of 
the R,HPO at all values of orthophosphate concentration. Thus a value for 
[R,] of 11.4 & 10-'° moles R,/0.5 gm roots is obtained when the value for ky, ‘s 
substituted in equation (13), which expresses the maximum absorption where [F, ] is 
completely saturated with HPO{; the [R,] is equal to [R,HPO; J: 


ao. = [Ro Ika» t + [R,]. (13) 


Essentially total saturation of R, is obtained at 5 X 10-4 M orthophosphate 
concentrations (where V’,,,/{2P] = 0.2 X 10~*) (Fig. 2) and [R,HPOj ] is approxi- 
mately equal to [R,]. The value of [R,H»PO,] is obtained by subtraction of [R,], 
since the intercept of the time curve for concentrations of orthophosphate of 5 X 
10-4 M is a measure of phosphate in combination with both sites (Fig. 1). Total 
uptake at 5 X 10-4 M [ZP] may therefore be expressed as 


Varo = | [RaH2PO, ha + [Ry liso} t + [R,H2PO,] + [Ry], (14) 


where numerical values of [R,H2PO,], [R,], ks are known approximately and the 
value of v’,4, is obtained from Figure 1. Solution of equation (14) gave a value of 
kz, of 1.4 & 10-* moles P/mole R, X sec. Substitution of this value into an ex- 
pression for Vinax @ results in a value for [R,| of 78 K 10-'° moles R,/0.5 gm roots. 
The relative concentration of R, to R, must be at least 1:7, a ratio which is of the 
same magnitude as the ratio of cytochrome b to DPNH in mammalian mitochron- 
dria. 

Actions of Competitive Inhibitors.—Adenosine diphosphate and pyrophosphate 
similarly inhibited orthophosphate uptake by barley roots. In the presence of a 
5 X 10-* M concentration of either inhibitor, a separation of the component reac- 
tions showed a change in K,» with no change in Vinax » (Fig. 3), with no measur- 
able effect on the H,PO{ uptake at this concentration. Within the range of in- 
hibitor concentrations from 1 XK 107% to 1 & 10~-! M, ADP and pyrophosphate 
decreased the uptake of both HPO{ and H,.POy. The inhibition was shown to be 
a function solely of a change in apparent Ky». and K,,». The inhibition of HePO; 
uptake became apparent only at the high concentrations of the inhibitor, because of 
a lesser affinity for the H,PO{ site than for the HPO site. The effect of either 
ADP or pyrophosphate was competitive to the uptake of both H,.POZ and HPO;. 
since the inhibition was abolished by increasing concentrations of orthophosphate. 
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Several schemes describing oxidative phosphorylation implicate an acceptor 
which forms an inorganic phosphate-ADP complex. Since R, and R, exist com- 
pletely in the forms R,H2PO, and R,HPO, at infinite orthophosphate concentra- 
tions, additions of ADP would be expected to increase the maximum uptake, 
Fins a+o, if such a complex were essential to uptake. However, aa «and Vinexd 
are unchanged by the additions of ADP. Either the endogenous ADP differs 
from the ADP in solution, or a phosphate-ADP complex is not required in the rate- 
limiting step of orthophosphate uptake by barley roots. 

The competitive inhibition by hydroxyl ion,t DNP," ADP, and pyrophosphate 
implies that these competitive inhibitors combine with the same sites involved in 
the rate-limiting steps of orthophosphate uptake: the reduced cytochrome b and 
DPNH. The relatively high affinity of the hydroxy] ion' suggests that the sites of 
phosphorylation may exist primarily in combination with the hydroxyl ion. Com- 
petitive inhibition of orthophosphate uptake implies a direct reaction of the ortho- 
phosphate species and the competitive inhibitors with the sites involved. First- 
order reaction kinetics cannot exclude the possibility that the actual competition 
occurs at some intermediates which are coupled to cytochrome b and DPNH. 

Inhibitors Which Block Electron Flow.—Azide is believed to block electron flow 
by combining with the oxidized form of the terminal oxidase, cytochrome a3, which 
résults in the components of the respiratory chain becoming more reduced.* The 
presence of a low concentration of azide (10~*-10-* M) caused a small increase in 
uptake from 10~-* M orthophosphate. The increased phosphate uptake was a func- 
tion primarily of an increased [R,], with no measurable change in the value of kz». 
A decrease in the apparent rate of breakdown of the phosphate-site complexes, 
ksa+», is obligatory upon an electron restriction, since oxidative phosphorylation is a 
direct function of electron flow. The apparent absence of a change in the value of 
kz is in accordance with Chance’s observation that large changes in the oxidation- 
reduction levels of components involved in electron transfer are observed with only 
a small inhibition of respiration over a considerable range of azide concentrations. 
With greater azide concentrations (10~* /), the rate of phosphate uptake decreased 
with increasing times of sorption (Fig. 4). This deviation from steady-state uptake 
was peculiar to this type of inhibitor. Inspection of the initial slope of the time 
curve, where steady-state uptake was approximated (Fig. 4), showed that a de- 
crease in [2A] had occurred. Further, a separation of the component reactions 
after 1-minute sorption periods ascribed the inhibition to a decrease of [R,]. In- 
creasing concentrations of orthophosphate did not overcome the inhibition. Thus, 
to be consistent with the previous interpretations, an oxidation of cytochrome b, 
R,, must be assumed at these concentrations of azide. In the presence of even 
greater concentrations of azide (10~-? M), a drastic decrease in the apparent k3a4» 
was evident, as well as an even further decrease in [ZR]. Separation of the com- 
ponent reactions showed that the phosphorylations coupled to both R, and R, are 
decreased. Increasing orthophosphate concentrations did not overcome the in- 
hibition in either case, and an oxidation of both cytochrome b and DPNH must have 
occurred. The contribution of HPO4 uptake is small in comparison with total up- 
take in the presence of 5 X 10-4 M orthophosphate. At this phosphate level, up- 
take is inhibited over the whole range of azide concentrations, with a significant 
decrease in both [R,] and ks,. 
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The substrate-inhibitor method for the examination of a respiratory chain is 
based on the premise that components of the respiratory chain on the oxygen side of 
the point of action of an inhibitor become more oxidized and those on the substrate 
side more reduced. The increased uptake of orthophosphate at low azide and phos- 
phate concentrations is in accordance with the reduction of the components of the 
respiratory chain induced by an azide block of electron flow at cytochrome az. 
However, higher concentrations of azide inhibit uptake by decreasing both the 
kay» and [R,4,]. A secondary site of action of azide with a respiratory component 
toward the substrate side of DPNH must be assumed on the basis of the substrate- 
inhibitor method. Further, the deviation from steady-state uptake with time sug- 
gests that azide combines with a state of the respiratory component which is a 
product arising from electron transfer. The conclusion is that azide combines 
with the oxidized form, DPN, as well as with cytochrome a;, on the basis that the 
current hypothesis adequately describes the components of the respiratory chain 
involved in oxidative phosphorylation (eq. [11]). 

The recent finding that cyanide combines with the oxidized form, DPN,'* sug- 
gested a study of the effect of cyanide on orthophosphate uptake. The charac- 
teristics of the cyanide inhibition were identical with those of azide, although slightly 
greater effective concentrations of cyanide were required to effect the same degree of 
inhibition as azide. 

Action of Substrates.—The respiratory chain responds differentially to the type of 
substrate which is dehydrogenated. The respiration rates of intact mammalian 
mitochondria, for instance, have been shown to be about three times larger with 
succinate as a substrate than with 8-hydroxybutyrate.'’? Further, a greater reduc- 
tion of components of the respiratory chain occurred with the greater rate of elec- 
tron flow effected by succinate.'* 

Mitochondrial preparations, after starvation periods of relatively short duration, 
essentially are depleted of endogenous substrates. However, excised barley roots 
are not amenable to a quantitative depletion of endogenous substrates, and addition 
of a high concentration of substrate is required for assurance that the measured 
response is a function primarily of the prevailing substrate. 

Succinate addition to an effective concentration of 1 X 10~! M decreased the 
uptake of orthophosphate by excised barley roots. Separation of the component 
reactions showed the inhibition to be essentially a function of a decrease in Vinax a- 
Measurement of phosphate uptake as a function of time showed that the slope de- 
creased to a lesser extent than the intercept (Fig. 5). Thus the effects of a high 
concentration of succinate on phosphate uptake were resolved as a decrease in [R, | 
and an increase in the apparent k;,4). The prevailing concentration of succinate 
severely restricts electron flow through the competing pathway of 8-hydroxybutyr- 
ate. The decrease in [R,], therefore, results from the DPNH becoming more 
oxidized with the slower rate of electron flow through the 8-hydroxybutyrate path- 
way. 

A significantly larger decrease in uptake was effected by a 1 K 107! M concen- 
tration of 8-hydroxybutyrate. An analysis of the effect, identical with that used 
for succinate, resulted in the conclusion that 6-hydroxybutyrate had decreased 
significantly both [R,] and the apparent k3.,,. A severe restriction of electron flow 
through the succinate pathway is effected by the competing pathway of 8-hydroxy 
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butyrate. This effect alone would not necessarily decrease [R,], since both cyto- 
chrome b and DPNH are involved in the pathway of the reducing reaction. How- 
ever, the decrease in k3,,, showed that a slower rate of electron flow was effected by 
this substrate than by succinate. The slower electron flux results in cytochrome b 
becoming more oxidized. 

These experimental observations are in qualitative agreement with the current 
hypothesis of oxidative phosphorylation in the electron transport system.® Inspec- 
tion of the chemical representation of the respiratory chain (eq. [11]) shows that 
the electrons from both the succinate and the 6-hydroxybutyrate pathways of reduc- 
ing reactions converge on cytochrome b. The dehydrogenation of a substrate such 
as 6-hydroxybutyrate, which requires DPN in the pathways of reducing reactions, 
results in the uptake of both the H:PO, and HPO, ions which are coupled to DPNH 
and reduced cytochrome b, respectively. The values of k;, and ky» must be equal, 
since the rate of breakdown of both R,H,PO, and R,HPO, is a direct function of 
electron flow, and the breakdown of both occurs with the transfer of an electron 
pair. The magnitude of their numerical value is a measure of the rate of electron 
flow effected by the 8-hydroxybutyrate pathway. 

; The succinate pathway does not involve 
DPNH but is linked to cytochrome b. In 
this instance, however, the numerical value 
of kz is a measure of the rate of electron flow 
effected by the succinate pathway. The 
value of k3,,, was shown to be increased by 
succinate and decreased by §-hydroxy- 
butyrate. A greater value of kz must have 
occurred in the presence of succinate than 
with 6-hydroxybutyrate, since kz, is a func- 
tion of 8-hydroxybutyrate dehydrogenation. 
Therefore, a more concise description must 
be accorded ky»—namely, ky» succinate and 
kx» g-onn- Lhe apparent k3, of control roots, 

2-08 , 107M utilizing an endogenous level of substrate, is 
e—e- 3 numerical value dependent on _ the 
6020180240800 Michaelis affinities, concentration, and ratio 


Fia. 5.—The uptake of orthophosphate of the oe spe 
from 0.5 gm. of excised barley roots from Approximate values for — kp succinate, 
1 X 10~° M orthophosphate in the pres- — f;,, , 4, and k;, may be calculated on the 


ence and absence of succinate or 6-hy- et ; é 
droxybutyrate. supposition that a high concentration of an 
added substrate effects a single pathway 
of electron flow. Values of [R,] and »,’, obtained from the succinate curve of Figure 
5, are substituted in equation (12), which expresses the uptake of the HPO, ion. 
Solution of the equation resulted in a value of 2.1 X 10-* moles P/mole R, X sec. 
Similar treatment of the measurements of uptake from 1 X 10~! M 6-hydroxybutyr- 
ate resulted in a value of 8.3.X 10-* moles P/mole R X sec. This value expresses 
either ks_ or ks» g-onn, Since both are equal in the 6-hydroxybutyrate pathway. 
The 2.5 times larger uptake of the HPO, ion with succinate than with 6-hydroxy- 
butyrate as a substrate is in close agreement with the 3 times larger respiration rate 
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shown with the same substrates by Chance and Williams for mammalian mito- 
chondria. 

Also, the value of 2.1 X 10-? moles P/mole R X sec for kp succinate 18 in essential 
agreement with the value of 2.3 X 10~-* moles P/mole R X sec for the kz, calculated 
from kinetic analysis of uptake measurements of control roots utilizing an endoge- 
nous level of substrate. However, the value of 8.3 K 10-* moles P/mole R X sec 
for either kz or ks g-onn is much lower than the value for k;, of 1.4 X 107? moles 
P/mole R X sec which was calculated from uptake measurements of control roots. 
A study of the effect on orthophosphate uptake of numerous substrates which re- 
quire DPN in the pathways of reducing reactions has resulted in the conclusion 
that the k;, of glutamate corresponds most closely to the k;, calculated from up- 
take measurements of excised barley roots utilizing an endogenous level of sub- 
strate. 


CONCLUSIONS 


Orthophosphate uptake by barley roots is coupled to sites identical with those for 
oxidative phosphorylation during electron transport in mammalian mitochondria. 
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THE NORM FORM OF A RATIONAL DIVISION ALGEBRA* 
By A. A. ALBERT 
YALE UNIVERSITY 
Communicated April 4, 1957 


Let D be a finite-dimensional associative division algebra over the field % of all 
rational numbers. Then the center § of D is an algebraic extension of finite degree 
over 8, and D is a cyclic algebra of degree n and dimension n? over §. We may 
then write D = (3, S, y), where 8 is a cyclic field of degree n over §, y is a nonzero 
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element of §, and S generates the cyclic automorphism group of 3 over §. The 
general element x of D now has the form x = xo + ry +--+ + 2,—1y"~!, where the 
vc, are in 3, and y in D has the property y" = y, yz = (2S)y. 
The algebra D now has an isomorphic representation as an algebra of n-rowed 
square matrices with elements in 8. It is the mapping x ~ X = Xo + X,Y + 
-- + X,_,Y"—'!, where X, = diag {x., ae aes a1} and 


Then we define the norm form of D to be the determinant 


A(x) = |X 


of the representation matrix. It is a homogeneous polynomial of degree n in 
the n? co-ordinates in § of the general element x of D, and it permits composition, 
that is 


A(xy) = A(a)-A(y). 


The purpose of this note is to provide a proof of the following result. 

THEOREM. Let D be a division algebra of odd prime degree over its center §, where 
& zs an algebraic number field. Then the norm form of D is universal, that is, every 
nonzero element 6 of § is the norm 6 = A(x) of an element x in D. 

Let us begin by recalling some of the properties of rational division algebras. ' 
We let & be the ring of all integers of §, so that every prime ideal $ of $ determines 
a valuation dy of §, and a completion §y of §. The scalar extension Dg of D by 
yg will then be a total matric algebra for all but a finite number* of valuations 
¢y. Let us call the primes $ for which Dg is a division algebra the critical primes 
for D. Then it is known? that a field &, of the given odd prime degree n over §, 
is a splitting field for D if and only if the degree ng of each composite (R, Fg) over 
& is n for every critical prime. Thus, if R = §[£], where & is a root of 


fw) = wo + aw"! + +++ + a, 


with a, in §, and f(w) is irreducible in ¥g for every critical prime §, it will follow 
that & splits D. Then & will be isomorphic to a subfield [x] of D, f(x) = 0, 
a, = (—1)" A(x). It should now be clear that our theorem is equivalent to the 
property that if a, is any nonzero element of §, there exist elements ay, +++, Gn 1n § 
such that f(w) is irreducible in &y for any prescribed finite set of prime ideals J of §. 
However, it will be convenient for us to remain closer to the theory of algebras 
than this remark might suggest. 

There is no loss of generality if we assume that a, is an integer of §, since a,= 
— A(x) if and only if \"a, = —A(Ax). We can clearly select \ so that \"a, is in- 
tegral. We shall also select the elements a, --- , @,—; to be algebraic integers, so 
that any root é of f(w) = 0 is always integral. Every element of § has the form 
a = mu, where m can be taken to be any element of $ not in $2, and yu is a unit 
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of the ring Sg of all integers of §g. Then A > 0 for every element of Sy. We 
take $ to be a critical prime for D and write 


On = 1 Br, h=nk+t (0<t<n), 


where 8, is a unit of Xy. If t > 0, we select the a; so that 
a; = 0 (mod $”"* +). (1) 


Then a; = x” +8, for integers 8; in Sg, and f(w) = w" + 2” (B,)9w"—! + Bow*w"-* + 
-->+ 6, :9"~!w + 7'B,). As we have said, f(w) has leading coefficient 1 and other 
coefficients in Yg, and so every root of f(w) is integral over $y. But then f(w) of 
prime degree is reducible if and only if it has a root p in Jy, and we may write 
p = n'a fora unit o of Sy. Clearly r”**! divides p”, and so vy = k + s, where we 
must haves > 0. Then x” is a factor of p” as well as of 8;2‘p”—', and so x” divides 
x'B,, Which is impossible. Thus the congruences (1) make f(w) irreducible in the 
field §g for every critical prime divisor $ of a, such that the exact power of $ divid- 
ing a, is $*"*+, where 0 < t < n. 

There remains the case where ¢ = 0, so that a, = 7*"8,, with 6, a unit of Yg. 
We assume that 


a, == w"*—"+H4y, (mod P"—"+#+1), (2) 


where the elements y; # 0 (mod §8) will be selected later. Then a; = 2™*-‘+#@,, 
where 8; is a unit of Yg, and may be determined modulo 7. But then f(w) = 
wo" + Byr™*®—M—D atl 4. +. + Bie —lw + 2B, has a root in §y if and only if 
a "*f(war*) = g(w) = w* + Bw"! + +++ + Brio + B, has a root in Fg. It is 
known‘ that the algebra Dg = (QW, 7’, +‘), where YW is a cyclic unramified field over 
ge with generating automorphism 7’, YW = Fg(h) for a root of unity h. The resi- 
due-class field Hg = Sg — (7) is a finite field of g elements, and Hg(h) has degree 
n over x, A is a primitive (¢" — 1)st root of unity, kT = h*. Every unit of Sg is 
known’ to be a norm in YW, and we observe that, if 


a= a = he-), 

then the norm of a as an element of Wis 1. However, a is never in §g, since other- 
wise the residue class of a has period gq — 1, h®-@-” = 1 (mod z), whereas the 
period of h is q” — 1 > (q — 1), since g”-!' + ---+q+1>q— 1 wheng > 2. 
But if 8, = —N(b) for b in BW, and 6 is in Fg, then 8, = —N(ab), where ab is not in 
dy. It follows that B, = —N(b), where b generates the cyclic field YW of degree n 
over fy. Then b is a root of go(w) = w" + aw"! + +++ + opsw + Bn, and we can 
take 


8; = a; (mod §). (3) 


Since YW is unramified, the polynomial [go(w)] = w" + [oi]w"~! + +++ + [onulw + 
[8,], over the residue-class ring Og, is irreducible. Then (3) implies that g(w) is 
irreducible, for a reducible polynomial cannot have an irreducible residue class. 
Thus the selection (3) makes f(w) irreducible in the cases of critical primes $ with 
a, exactly divisible by $B". Since we have a finite number of congruences to be 
satisfied by the a,;, we can select them so that f(w) is irreducible, as desired. 
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This completes our proof in the case where § is an extension of finite degree over 
®. The result is true, however, even when § is infinite-dimensional over R. For 
D has a basis uw; = 1, ue, +++, Um over § and a multiplication table ujuj = yijxnu, for 
algebraic numbers y;, where m = n?. If 6 ¥ 0 is in §, we let f be the field ob- 
tained by adjoining the y;, and éto®. Then D = (Do)g, where Do = who + --- 
+ Um ¥o is a central division algebra of degree n over %. The field ¥ is a finite 
algebraic extension of R, and soé = A(x), where x isin Do. But A(z) is the norm 
form of D, and our result is proved. 


*This paper was sponsored by the Office of Ordnance Research, United States Army, under 
Contract DA-19-059-OR D-2329. 

1 An exposition of this theory is given as chap. IX of the author’s Structure of Algebras. 

2 See Theorem 9.30. 

3 This is a consequence of Theorems 9.31 and 9.23. 

4See Theorem 9.21. 

5 See Theorem 9.19. 


THE RANDOM FUNCTIONS OF COSMIC-RAY CASCADES 
By T. E. Harris 


RAND CORPORATION 
Communicated by S. Bochner, April 2, 1956 


1. Introduction.—We consider the electron-photon cascades of cosmic radiation 
in the formulation of Bhabha and Heitler! and of Carlson and Oppenheimer.? We 
make the simplifying assumptions (used in some of the physical literature) appro- 
priate to the so-called ‘complete screening” case. Hence the asymptotic results 
given below cannot be strictly true, since our assumptions hold only over limited 
ranges of the parameters. Presumably, detailed analytical study of more com- 
plicated models would be necessary to determine whether the results have some 
qualitative validity. 

A photon of positive energy ¢«, moving through homogeneous material, has 
probability \ dt + o(dt) (henceforth we shall omit the o-terms in such expressions) 
of being transformed in the thickness interval (¢, ¢ + dt) into two electrons, posi- 
tive and negative, which receive respective energies ew and «(1 — u),0< u <1, 
with probability density q(u). We assume that g has a bounded derivative for 
0 < wu < 1 and is symmetric about '/z. An electron loses (by “‘collision” or “ioniza- 
tion’’) the deterministic amount of energy 6¢ in any interval of length ¢, provided 
that the energy is not thereby reduced below 0; 8 and X are constants independent 
of t and of the energy. Also, an electron radiates photons (‘““Bremsstrahlung”’) : 
the probability that an electron of energy ¢ emits a photon of energy between 
eu and e(u + du) in the interval dt is k(u) du dt, the energy which goes to the photon 
being subtracted from that of the electron. The usual simplifying assumptions are 
covered if we take 


Ik) 
ku) ="+h(u), i <el—w, O<u<il, (1.1) 
u du | 


' | 
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1 
where c, u are positive constants and b < 2. Since JS k(u) du = © (infinite cross- 
0 


section), with probability 1 infinitely many photons are produced by an electron 
in each ¢-interval in which it has positive energy. 

2. The Straggling’ and Generating Functions When 8 = 0.—Let e(t) be the 
energy at ¢ of an electron which had initially energy 1, in the case 8 = 0. Putting 
X)(t) = —log e(t), it can be seen that Xo(¢) has for t > 0 a distribution of the 
infinitely divisible type with a probability density h,(x) whose characteristic func- 
tion is, for x > 0, 


Sf hlzje de = exp J, S (e'™ — I)k(1 — ee du. 2.1) 
é Uo ) 
The case k(u) = —p/log (1 — u), where h,(x) is then x“!~'e-*/T (ut), was given by 
Bethe and Heitler.® 

Let N(e, ¢) be the total number of electrons at ¢ whose energies are greater than 
e, and let f,(s, ¢, t) and f2(s, ¢, t) be the generating functions, )>P(N(«, t) = r)s’, if 


initially (i.e., at £ = 0) there is, respectively, one photon or one electron of energy 1. 
To derive equations for f; and fo, the usual method of regeneration points must be 
modified, since there is no first photon emission. Since differential equations 
for f; and fy have not previously been justified in the present case of infinite cross- 
sections, we mention that differentiability in ¢ and ¢ of f,; and f. can be established 
by direct consideration of the random process, and we get equation (2.3) below; 
the more simply derived companion equation (2.2) is also given: 


Ofi(s, ¢, t : 
fi(s, ¢, t) snide 44) > VR (s, ott ) fr (s, 4 8 -, ) q(u) du, (2.2) 
di é u l1—wu 


. 1 ( 
ols, ae es S ih (s, a : fo (s, i ’ : aR fo(s, €, ot k(u) du, (2.3) 
ot 0 u l—u 
with fi(s, «, 0) = fi(s, 1,0 = fo(s, 1, t) = 1, t > 0, and fo(s, «, 0) = s for e < 1. 
Uniqueness of the solution among a wide class of functions follows from the fact 
that the factorial moments (derivatives of the f; with respect to s at s = 1) have 
Mellin transforms on ¢ which satisfy linear differential equations. 
1 


In the case of a finite total cross-section (f° k(u) du < ), integral equations 
0 


were given for the characteristic functional of N(e, t) by Bartlett and Kendall® and 

for the generating functions f; and f. by Janossy’ and others. Integral equations 

can be obtained in the present case using a method similar to the one employed by 
Moyal* to get his equation (2.22). 

3. The Expectation Process When 8 = 0.—Define a vector Markov process 

1 


(I(t), ¢(t)), t > 0, with J = 1 or 2 and ¢ > 0, by the following scheme (let K = { 
0 


uk(u) du). I(t) is itself a two-state temporally homogeneous Markov process (we 
can consider ¢ as a “time” rather than a “thickness” parameter) with respective 
probabilities \ dt and K dt for the transitions 1 — 2 and 2— 1indt. Whenever the 
transition 1 — 2 occurs, ¢(t) is multiplied by u, 0 < u < 1, with probability density 
2uq(u), and whenever 2 > 1 occurs, ¢(t) is multiplied by u with density uk(u)/K; 
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furthermore, if /(t) = 2 fort; < t < th, then ¢(t)/¢(4) decreases, in the interval 
t; < t < t, in the manner of the straggling process of Section 2, except that the 
function k(u) is replaced by (1 — u)k(u); ¢(¢) remains constant in any ¢t-interval in 
which /(t) = 1. We may think of /(t) and ¢(¢) as describing the condition of a 
single particle which can be sometimes a photon (J = 1) and sometimes an elec- 
tron (J = 2) and which has an energy ¢(t). We call (J(d), ¢()) the “expectation 
process,”’ for reasons to be seen below. 

If we take ¢(0) = 1 and J(0) = 7,7 = 1 or 2, then if 0 < « < 1, (/(0), ¢(#)) has 
for t > 0 a probability density p,,(e, t) de for the probability that 7(t) = 7 and 
e<f(t)}<e+de. Let m,,(e, t) de be the expected number, for the cascade process 
with 8 = 0, of particles of type j at ¢ with energies between ¢ and ¢ + de, starting 
with one particle of type i and energy 1; 7, 7 = 1 for photons and 2 for electrons. 
(These expectation densities are easily seen to exist for t > 0,0 < e< 1.) 

TuHeoreM. We have* p,,(¢,t) = emi(e, t). 

This can be seen by comparing differential equations in ¢ for p,; with those for 
m;;; the latter can be obtained from equations (2.2) and (2.3) and the corresponding 
equations for the number of photons. Similarly related Markov processes can, 
of course, be found for other branching processes. 

Since for the expectation process there is, with probability 1, a smallest ¢ where J 
changes, we have a natural means of using the regeneration-point method to get 
integral equations for the p,; and hence for the m,,; the equations obtained in this 
way are different in form from those mentioned at the end of Section 2. Here we 
merely mention that these integral equations are convenient for proving differenti- 
ability properties in ¢ of the mi or pi. The correspondence between the two 
processes might be useful in artificial sampling (Monte Carlo) experiments; using 
the expectation process to determine p and hence m, one wouid deal with less com- 
plicated histories than if the original process were sampled directly. 

4. Steady States and Limits for the Case 8 = 0.—The process (/(t), ¢(8)) has the 
stationary “distribution’’’® which assigns the weights K(a, — a;) and A(az — a;) 
to the events {/(t) = 1, a < —log ¢(t) < a>} and (T(t) = 2,a, < —log ¢(t) < ar}. 
Hence the original cascade has the stationary densities K de/e? and \ de/e? for the 
expected number of photons or electrons in the energy range e, « + de. For par- 
ticles of a single type, a corresponding result was given by Janossy."! 

Next let z(t) be the total energy at ¢ in all electrons, for the cascade process with 
8=0. Since E[z(t)] = P[(t) = 2], and since P[/(t) = 2] has from the theory of 
Markov processes the limit \/(A + K) ast— ©, we have (the symbol EF denotes 
expected value) 


lim £[z (t)] = 


41 
ps + K’ (4-1) 


a result given by Bhabha and Chakrabarty.’ We can strengthen equation (4.1) to 
get convergence in probability also.'? In fact, 


wm af) - 2] =o « 
im |} z ie naan tas = U. 2 
t+ @ r + K ( 

* 5. Cascades with 8 > 0.—Let e(t) be a random function defined as in Section 2; 
for any 8 > 0, define*® 
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€,(t) = Max {o, a(0| 1 - B S ale (5.1) 


0 (8) 


This formula, which has apparently not been noticed in the literature, gives a 
representation for the straggling when @ is the collision factor. It may be of use in 
sampling experiments, since it gives a convenient way of constructing sample «- 
functions corresponding to various values of 8, when a single sample ¢-function is 
available.'4 

Now let V(t) be the total number of electrons of all energies at t, electrons whose 
energies have reached 0 not being counted. It must be recalled that we are using 
the usual simplifying assumption that the functions q and k of Section 1 do not 
depend on energy, so that no lower limit is assumed for the energy of an electron, 
whereas in actuality an electron has a certain minimum (rest) energy. However, 
even for this simplified model Bhabha and Chakrabarty* obtained from the moment 
equations the result (deducible also from energy considerations) that EH [N(t)] is 
finite when 8 > 0. We can strengthen this to the result that #[N?(t)] is also finite 
when 8 >0. However, a direct analysis of the successive generations of the process 
shows that if there is one initial electron, then, with probability 1 N(t) reaches © in 
every t-interval. Although this result must be regarded as a. peculiarity of the 
mathematical model (the so-called “Approximation B”’), it suggests that in reality 
the distribution of the maximum of N(¢) over a ¢t-interval may differ markedly from 
the distribution for fixed t. 


1H. J. Bhabha and W. Heitler, ‘Fast Electrons and the Theory of Cosmic Showers,’ Proc. 
Roy. Soc. London, A, 159, 432-458, 1937. 

2 J. F. Carlson and J. R. Oppenheimer, “On Multiplicative Showers,’’ Phys. Rev., 51, 220-231, 
1937. 

3H. J. Bhabha and 8. K. Chakrabarty, ‘“The Cascade Theory with Collision Loss,’’ Proc. 
Roy. Soc. London, A, 181, 267-302, 1943. 

4 “Straggling’’ is the name applied to the random decrease in energy of an electron. 

5 H. Bethe and W. Heitler, ‘On the Stopping of Fast Particles and on the Creation of Positive 
Electrons,’’ Proc. Roy. Soc. London, A, 146, 83-112, 1934. 

6 M.S. Bartlett and D. G. Kendall, “On the Use of the Characteristic Functional in the Anal- 
ysis of Some Stochastic Processes Occurring in Physics and Biology,’’ Proc. Cambridge Phil. Soc., 
47, 65-76, 1951. 

7 L. Jdnossy, “Studies on the Theory of Caseades,’’ Acta phys. Acad. sci. Hung., 2, 289-333, 
1952; this article mentions earlier work of JAnossy. 

8 J. KE. Moyal, ‘Theory of the Ionization Casecade,’’ Nuclear Physics, 1, 180-195, 1956. 

® Both p and m have versions which are continuous in e; taking these, the equality holds for 
all, not just almost all, e. 

© On stationary ‘‘distributions’’ see T. E. Harris, ““The Existence of Stationary Measures for 
Certain Markov Processes,’’ Proc. Third Berkeley Symposium Math. Statistics and Probability, 2, 
(University of California Press, 1956), 113-124. 

|. Jdnossy, ‘On the Absorption of a Nucleon Cascade,’’ Proc. Roy. Irish Acad., A, 53, 181- 
188, 1950. 

2 The random process e‘X+)*{2(t) — »/(K + 2)] is a martingale but possibly not a con- 
vergent one; at least its mean square is unbounded at t— ©. Equation (4.2) can be shown by 
considering a functional equation for the characteristic function of z(t). 

18 The e-process is defined in such a way that with probability 1 the sample functions are posi- 
tive nonincreasing step functions (although with infinitely many jumps in each ¢-interval). 

14 For a direct study of the distribution of « see A. Ramakrishnan and P. M. Mathews, “Strag- 
gling of the Range of Fast Particles as a Stochastic Process,’’ Proc. Indian Acad. Sci., 41, 202-209, 
1956. 





THE IDENTITY OF THE WEAK AND STRONG EXTENSIONS OF 
A LINEAR ELLIPTIC DIFFERENTIAL OPERATOR 


By M.S. NARASIMHAN 
TATA INSTITUTE OF FUNDAMENTAL RESEARCH, BOMBAY, INDIA 
Communicated by S. Bochner, March 12, 1957 


1. Let Q be an open subset of the Euclidean space R". Let 
ot: re | 


D = » > G;,.. jn (2) 


+ £ in 
jit...tjnsm Oa". le 


be a linear differential operator in Q with indefinitely differentiable complex-valued 
coefficients aj, jn(x), = (a1, ... , Xn) €Q. D operates on the space of distribu- 
tions in Q. Let D, denote the restriction of D to the space of C” functions with 
compact supports. The weak extension Dw of D, in L*(Q) is defined as follows: the 
domain M of Dy consists of all those elements f in L* for which Df, formed in the 
sense of distributions, is in Z?, and for such f, Dwf is defined to be Df. The strong 
extension Ds of D, in L? is the closure in L? of the operator Do, where D, is defined 
as follows: its domain consists of all fin L? which are C” and for which Df, formed 
in the usual sense, belongs to L*, and for such f, Dof = Df. It is easily seen that 
Dy is an extension of Ds. It is not known whether, in general, the weak and strong 
extensions coincide. In this note we prove their identity for a class of differential 
operators. ! 

2. TueroreM 1. Let Q be an arbitrary open subset of R", and let D be formally 
self-adjoint and elliptic.2 The weak and strong extensions of D, in L?(Q) coincide. 

Proof: By assumption, D, is a symmetric operator in L*; also, Dy is the adjoint 
of D, in L?. We have, by a theorem of von Neumann,’ 


Domain of Dw = domain of the closure of D,; @ E(i) @ E(—2), (1) 


where #(z) and E(—7) denote the spaces of all square-integrable solutions of the 
equations DT = iT and DT = —iT, respectively. From the results of F. John 
and L. Schwartz,? the elements of E(i) and E(—7) are C” functions. It follows 
immediately from this and from decomposition (1) that the weak and strong ex- 
tensions coincide. 

3. THEOREM 2. Assume that Q is bounded. Let D be uniformly elliptic, * not 
necessarily formally self-adjoint. Then the weak and strong extensions of D, coincide. 

Proof: Since the identity of (D + A)s and (D + X)w, for A real, implies that of 
Ds and Dy, it is sufficient to prove that (D + A)s = (D+ A)wfor some X. Let 
D+X=D’. The domain M of Dy is a Hilbert space with the scalar product 


(fi, fo)) = (fi, fo)uz + (D’fi, D’ fo) x1, 


and we have to show that M n & is dense in M in this topology, § denoting the space 
of C” functions in. GAarding’s inequality” ‘ allows one to define the generalized 
Friedrich’s extension‘~* of D,; we denote its domain by NV. We consider N as a 
Hilbert space with the scalar product induced from M. Then, for X sufficiently 
large, D’ | N:N — L’ is a topological isomorphism of N onto L?.5* Since C* func- 
tions in L? are dense in L? and since (D’| N)-!: L?-—>» N maps C” functions into C~ 
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functions by the elliptic character of D, it follows therefore that N n & is dense 
inN. Nowiff eM, we can find a unique f’ « N such that D’f = D’f’. It follows 
immediately that 


M=N O01, 


where L is the space of all square-integrable solutions of the equation D’f = 0. 
Again, by the elliptic character of D, L ¢ & So, from decomposition (2), we see 
that M n &is dense in M. 

Remark 1: A similar argument shows that if, in Theorem 2, the coefficients are 
analytic, Dy is the closure of its restriction to analytic functions in the domain of 
Dy. 

Remark 2: Actually, for the proof of Theorem 2, it is sufficient to know that D’ 
maps M onto L*. For, by Banach’s open mapping theorem, the map D’| M:M—> 
L? will then be an open mapping, and the inverse image of the set & n L?, every- 
where dense in L?, will be everywhere dense in M. 


1 For the case of second-order elliptic operators see M.S. Birman, Doklady Akad. Nauk S.S.S.R., 
N.S., 92, 205-208, 1953; and for that of a class of operators with constant coefficients see L. 
Hérmander, Acta Math., 94, 241, 1955. 

2 L. Garding, Math. Scand., 1, 55-72, 1953. 

3 J. von Neumann, Math. Ann., 102, 85-86, 1930. 

4L. Nirenberg, Communs. Pure and Appl. Math., 8, 648-684, 1955, in particular pp. 657-659. 

5 J. L. Lions, Acta Math., 94, 26, 1955. 

6 P. D. Laxand A. N. Milgram, ‘Contributions to the Theory of Partial Differential Equations,” 


Ann. Math. Studies, 33, 167-190, 1954. 


ON THE PRINCIPLE OF INVARIANT IMBEDDING AND DIFFUSE 
REFLECTION FROM CYLINDRICAL REGIONS 


By RicHarp BELLMAN AND RoBert KALABA 
THE RAND CORPORATION, SANTA MONICA, CALIFORNIA 
Communicated by S. Chandrasekhar, April 11, 1957 


1. Introduction.—The importance of various principles of invariance in treat- 
ing problems of radiative transfer has been discussed by Ambarzumian” *? and 
Chandrasekhar,* who applied these principles to a variety of problems connected 
with plane-parallel regions. However, the application of these concepts to the 
corresponding problems for cylindrical and spherical regions is not immediate. 
Chandrasekhar has observed: ‘The extension of the approximate methods of 
solution to other geometries is straightforward; but it is not equally apparent what 
the general principles are which will play the same unique role as the principle of 
invariance in plane-parallel atmospheres.’’4 

In our first paper® we introduced the principle of invariant imbedding as an ex- 
tension of the ideas of Ambarzumian and Chandrasekhar and utilized this principle 
to treat the case of inhomogeneous plane-parallel regions. The regions are inhomo- 
geneous in the sense that the phase function, the absorption coefficient and the 
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emission coefficient, to use the terminology of Chandrasekhar, are all allowed to be 
functions of position. In this paper we propose to show that this new concept 
enables us to treat cylindrical and, analogously, spherical geometries in a unitary 
fashion and thus provides an answer to the query of Chandrasekhar. Analytic and 
computational questions arising from this treatment will be discussed in greater de- 
tail in a subsequent publication. 

2. The Physical Model.—Let us now consider the problem of determining the 
intensity of diffusely reflected radiation from the surface of a circular cylinder. 
We assume that the symmetry properties of the medium are such that it suffices to 
consider interactions in a cross-section perpendicular to the axis, that the absorp- 
tion and scattering parameters of elemental volumes depend only upon distance 
from the axis, and, finally, that the scattering is isotropic. 

Consider an incoming flux of radiation making a constant angle @ with the radius 
of the circle, as indicated in Figure 1, and taken to have intensity 1 per unit length 
normal to the direction of propagation. This is the appropriate definition of flux 
for circular regions corresponding to the usual one for plane regions. As in the 
plane case, the ratio of scattering coefficient to absorption coefficient is denoted by 
d(x) /(1 — A(x)), where z is the distance from the center of the circle. 

We now consider the intensity of flux of reflected radiation making a fixed angle 
y with the extended radius of a cylinder of radius x, due to an incoming flux at 
angle @ (Fig. 1). This intensity is considered to be a function r(v, u, x) of u = cos 8, 
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v =cosy,andz. In this way we imbed the particular process within a continuous 
family of processes. 

3. Application of the Principle of Invariant Imbedding.—In order to obtain a 
functional equation for r(v, u, x), we consider the scattering from a cylinder of 
radius x + dz as the result of scattering from the infinitesimal annulus (x, « + dz) 
and the cylinder of radius x (see Fig. 2). Although we are definitely considering 
multiple scattering, our functional equations are considerably simplified by the 
fact that, to higher-order terms, only one scattering can occur within the annulus. 

Proceeding as in our first paper,’ which in turn follows the procedures of Am- 
barzumian and Chandrasekhar, we obtain the following nonlinear integro-dif- 
ferential equation :® 


1 1 1 — v? 1 — wu? d(x) 
Ft -+}- dy, Cammecaren Fett !. = — + 
U v xv ru 4v 


A(x) 


WA) 


= 0 


lz 


A(x 1 0 
r(w, u, x)dw + \ ' fir, ia + 
2 0 Zz 


1 1 dz 
d(x) f r(w, u, x)dw f r(v, x, 2) - (3.1) 
0 0 


together with the boundary condition, r(v, u, 0) = 0. 

This differs from the corresponding equation for the inhomogeneous parallel- 
plane case, derived in our first paper,® by the presence of the partial derivatives 
with respect tou and v. These terms arise from the fact that the radiation incident 
at angle @ on a circle of radius x + dz is incident, to first-order terms, at an angle 
x + (tan 6) dx/x on a circle of radius x. 

4. Reduction in Dimensionality.—As the equation stands, it is an equation for a 
function of three variables, r(v, u, x), two spacelike variables, wu and v, and one time- 
like variable, x. Let us now show that this equation may be replaced by two equa- 
tions for the functions P(v, x) and Q(u, x), defined as 


R d. U 
P(v, x) =So r(v, z, 2) = (4.1) 


Q(u, x) =Sh r(w, u, x) dw. 
Equation (3.1) has the form 
L(r) = F(PQ, x), Q(v, x)), (4.2) 


where L is a linear partial differential operator. Solving for r in terms of F, using 
characteristic theory or otherwise, we obtain a relation of the form 


T= A(F), (4.3) 


where A is a linear operator. Using this relation, we form P(v, x) and Q(u, x) as 
defined above. 


1 
Pv, 2) f ACP) ©, 
0 4 


Q(u, x) = {2 A(P) dw. 


(4.4) 
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This last set of equations constitutes a system of two nonlinear integral equations 
for the functions P(v, x) and Q(u, 2). 

Computational and analytic aspects of these equations will be presented subse- 
quently. 

1V. A. Ambarzumian, “Diffuse Reflection of Light by a Foggy Medium,’’ Compt. rend. Acad. 
Sct. U.R.S.S., 38, 229-232, 1943. 

2V. A. Ambarzumian, E. R. Mustel’, A. B. Severny, and V. V. Sobolev, T'eoreticheskaia Astro- 
fizika (Moscow, 1952), chap. viii. 

8 §. Chandrasekhar, Radiative Transfer (Oxford, 1950). 

4 Thid., p. 364. 

5 R. Beliman and R. Kalaba, “On the Principle of Invariant Imbedding and Propagation 
through Inhomogeneous Media,’’ these PRocEEDINGS, 42, 629-632, 1956. 

6 To simplify the analysis, we are considering only the case of spherinal scattering here. 


ON THE PRINCIPLE OF INVARIANT IMBEDDING AND 
ONE-DIMENSIONAL NEUTRON MULTIPLICATION 
By Ricuarp BELLMAN, Rospert KALABA, AND G. Mitton WING 


RAND CORPORATION, SANTA MONICA, CALIFORNIA, AND LOS ALAMOS SCIENTIFIC LABORATORY, LOS 
ALAMOS, NEW MEXICO 


Communicated by S. Chandrasekhar, April 11, 1957 


1: Introduction.—In this paper we wish to introduce a new method of treating 
problems involving neutron multiplication by fission, with special regard to ques- 


tions of critical mass and distribution of neutrons. Our results derive from ap- 
plications of the principle of invariant imbedding, which we have introduced and 
discussed in two previous papers," ? the first of which contains a brief discussion 
of the close interrelation between the method we employ here and the point-of- 
regeneration technique employed in Bellman and Harris.* 

In order to present the essentials of the technique, as free as possible of analytic 
difficulties, we shall consider first only one-dimensional flux of neutrons. This 
case is of interest in itself, from both the physical and the mathematical points of 
view. 

We shall consider first a one-velocity case, deriving partial differential equations 
for the generating functions of the probability distribution of ‘‘reflected’” and 
“‘transmitted” neutrons (those terms will be defined precisely below) and ordinary 
differential equations for the moments. Although the moment equations can be 
used to determine critical length, an inspection of the first and second moments 
shows that these yield only a very imperfect picture of the actual distribution of 
neutrons. Classical treatment is based on expected values. 

We then turn to a two-velocity case, to illustrate the applicability of our methods 
to the more realistic situations where fission and fission products are energy- 
dependent, as well as situations where several types of particles are involved. 

It will be clear from our discussion and from previous papers how our procedure 
can be modified to take into account absorption of neutrons and inhomogeneity of 
the medium. 

In subsequent papers we shall first consider the dependence on time of the dis- 
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tribution of neutrons and the case of continuous dependence on energy for the one- 
dimensional case and then turn to a discussion of the two-dimensional plane slab, 
the cylindrical case, and the spherical case. 

2. Description of the Process.—Consider a finite one-dimensional “rod,” the 
interval [x, 0], constituted of a material with the property that a neutron traversing 
an infinitesimal interval, [y + dy, y], has, to first-order terms, the probability 
a dy of splitting into two neutrons, of the same type as the original, one going to the 
left and one to the right. 

Assuming that a neutron, the “‘trigger’’ neutron, enters from the left, as dia- 
gramed in Figure 1, we wish to determine the number of neutrons emerging from 
the left (“reflected’”’ neutrons) and the number of neutrons emerging from the 
right (“‘transmitted”’ neutrons). 





> | | 
x O 





X+A 
Fig. 2 


3. Mathematical Formulation.—Let us define the two sequences of probabilities 
; | Pn(x)}, }n(x)} , as follows: 


p,(x) = the probability that n particles will be reflected from an 
interval of length 2, over all time, as a result of one neutron 
entering at x at time 0. 

n(x) the probability that n particles will be transmitted through 
an interval of length x, over all time, as a result of one neu- 
tron entering at x at time 0. 


The essence of the technique of invariant imbedding is to regard the process 
described above as a member of a family of processes of the same type. The 
mathematical statement of this invariance will yield a set of equations governing 
the process. This is accomplished readily in this case by regarding x as a parameter 
capable of assuming all positive values. 

4. Basic Functional Equations.—Let us begin by deriving relations connecting 
members of the sequences {p,(z + A)} and {p,(x)} for small A. Referring to 
Figure 2, we see that, to first-order terms, 


Po(x + A) = (1 = aA) po(x). (4.1) 


To obtain the equations for p,(2 + A), we consider reflection of n neutrons from 
[x + A, 0] as composed of the following occurrences: 
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a) Transmission of the neutron through [xz + A, x], reflection of n 
neutrons from [z, 0], transmission of n neutrons through [z, « + Al]. 

b) Transmission of the neutron through [* + A, x], reflection of k 
neutrons from [x, 0], k = 1, 2,..., n, fission of one of these k (4.2) 
neutrons in [r, 2 + A], reflection of (n — k) additional neutrons 
from [x, 0]. 

c) Fission of the neutron in [zx + A, x], reflection of (n — 1) neutrons 
from [z, 0]. 


This verbal statement is equivalent to the mathematical relation 


Pala + A) = (1 — aA) [p,(x) (1 — aA)" + aA YS kp, (x)pp_x(x)] + aAp, (2), 
k=1 
(4.3) 


to first-order terms in A.4 
Letting A — 0, the limiting equations are 


po (x) = —ap(x), 
Pn (x) = —(n + lap, (x) + apya(x) +a Yo kp,(x)pr«(2), n>. (4.4) 
k=1 


Similarly, we obtain the recurrence relations 


Qn (x) = —angn(x) +a D> kqx(x)pnx(x), n SS ae 
k=1 


with go(x) = 0. 
5. Partial Differential Equation for Generating Functions.—Let us define 


9) 2 


u(x, r) = Pr(x)r’, v(x, r) 5 » Qn(x)r”. 
0 1 


n= r= 
The recurrence relations of equations (4.4) and (4.5) yield the equations 


ur = au(r — 1) + aru,(u — 1), 


RO 
v, = ar(u — 1)p,. (5.2) 


These equations can be explicitly resolved using the method of characteristics. 
The details are not trivial, however. We shall discuss these matters in detail 
subsequently. 

6. Critical Length and Moment Equation.—The concept of “critical length” can 
be introduced as the smallest value of x for which the expected number of reflected 
and transmitted neutrons is infinite. 

Setting 


! 


U(z) = uj, <1, V(x) = 2,|, =1, (6.1) 
we obtain from equations (5.2), via differentiation and setting r = 1, the equations 


U’(x) = a(l + U®), V’ (x) = aU, (6.2) 


with U(O) = 0, V(O) = 1. From these equations we obtain U(r) = tan az, 
V(x) = sec ax, showing that the critical length is r/2a. 
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In a similar fashion, we find the differential equation for the second moment, 


qs _ oat + 2aU* + 3aUUr, U0) = 0, (6.3) 
dx 
which shows that U, = 0(1/(cos az)*) asx — w/2a. It follows that the expected 
number of neutrons is an unreliable measure of the nature of the distribution for z 
near /2a. 

7. Two-Velocity Case——Following the usual approximation to the physical 
situation where a neutron possesses a direction and an energy varying continuously 
over certain limits, let us assume that there are just two types of neutrons, “fast”’ 
and “slow,” and that in the fission process either can give rise to the other. 

To simplify the equations, let us assume that fast neutrons have a probability 
a,dy of splitting in |y + dy, y} and that, when they do split, there is a probability 
'/, of a fast neutron produced going to the right and a slow neutron produced going 
to the left, and probability '/. of a slow neutron going to the right and a fast neutron 
to the left. Let the same situation prevail with slow neutrons with a, replaced 
by Qs. 

Let 


Pn (a) = the probability that m fast neutrons and n slow neu- 
trons will be reflected from an interval of length z, over- (7.1) 
all time, as a result of one fast neutron entering at x at 
time 0, with pp,°>() defined similarly for a slow trigger 
neutron. 


Then, as above, we have 


Pmn (a + A) = (1 — aA) [Pm (x) (L — ad)" + Spe (x) ['/ePm—a, nt? (2) 
+ kay! 2Pm—k. n- 1 (x) }] : 
+ aA[! 2[Pm—1, n(x) + Pm. n—l (x) ]], (7.2) 


with a similar equation for Pm’ (x, + A). 

Passing to the limit, we obtain recurrence relations which permit us to obtain 
partial differential equations for the generating functions, u(x, r, t) = 
Lpmn(ayr™"t” and u(x, r,t) = Lpmn°?(x)r"t". From these we obtain ordinary 
differential equations for the moments and thus can determine critical length. 
At present we have not been able to solve these equations explicitly, though nu- 
merical solution would be straightforward. 


1 R. Bellman and R. Kalaba, “On the Principle of Invariant Imbedding and Propagation 
through Inhomogeneous Media,’’ these PROCEEDINGS, 42, 629-632, 1956. 

2 R. Bellman and R. Kalaba, ‘On the Principle of Invariant Imbedding and Diffuse Reflection 
from Cylindrical Regions’’ (Paper P-976 [Rand Corporation, 1956} ). 

’ R. Bellman and T. E. Harris, “On Age-dependent Binary Branching Processes,’’ Ann. Math. 
55, 280-295, 1952. For a thorough discussion of branching processes see T. E. Harris, “Some 
Mathematical Models for Branching Processes,’’ Proc. Second Berkeley Symposium Math. Statis- 
tics and Probability (Berkeley: University of California Press, 1951), pp. 305-328. 

‘ This is equivalent to the fact that in an infinitesimal interval it is sufficient to consider only 
first-order processes. 





ON THE EXPANSION OF FUNCTIONS WHICH SATISFY 
FOUR BOUNDARY CONDITIONS 


By S. CHANDRASEKHAR AND W. H. Retp 
UNIVERSITY OF CHICAGO 


Communicated April 27, 1957 


1. Introduction.—In the solution of many problems in applied mathematics 
it is often convenient to have expansions for functions which vanish, together 
with their derivatives, at the end points. In some recent work on hydrodynamic 
and hydromagnetic stability we have found certain types of expansions for such 
functions which have proved to be very useful and which appear to be novel in this 
connection. In this paper two types of such expansions will be considered, and the 
principles underlying them will be described. 

2. Fourier-Type Expansions.—Consider the characteristic-value problem de- 
fined by the equation 

BF Lig (1) 
dat 
and the boundary conditions 


d 
y=0 and w=6. et 8 (2) 
ax y 
Let ,, and \,, denote two distinct characteristic values; then, by using the boundary 
conditions (2). it can readily be shown that 


+1/2 +'/2 q2 1? ’ 
Am! f YmYn AX = f vd = a dex. (3) 
—'/, dx? dx? 


By subtracting from this equation the corresponding equation with m and n inter- 
changed, we obtain 


+1/2 
On* ~ Ax‘) f YmYn dx = 0; (4) 


and hence 


+1/s 
{ YmYn dx = 0 if m # n. (5) 


The solutions y, form, therefore, a complete orthogonal set of functions. We shall 
now obtain the explicit form of these functions. 

The solutions of equation (1) fall into two noncombining groups: the even 
and the odd solutions, with cosine- and sine-like behavior, respectively. Consider- 
ing first the even solutions of equation (1), we obtain the eigenfunctions 

cosh \,2 COs Av 


Ca(z) = cosh (An/2) cos (A»/2)’ = 


where the eigenvalues \, are the positive roots of the equation 
tanh (A/2) + tan (A/2) = 0; 
521 
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the corresponding odd solutions of equation (1) lead to the eigenfunctions 


iy = sinh p,x _ sin Ment 
A t)=- 
sinh (u,/2) sin (u,/2 


2 


where the eigenvalues u, are the positive roots of the equation 
tanh (u/2) — tan (u/2) = 0. (9) 


We shall presently see that, when chosen in this manner, the functions C, and S,, 
are normalized to unity. The first few roots of the characteristic equations (7) 
and (9) have been computed by Rayleigh! and are given in Table | for convenience. 
For large values of n, \, —> (2n — 1/2) and uw, — (2n + 1/2) r 


TABLE 1 
An Mn 
4.7300408 7 .8532046 
10.9956078 14.1371655 
17 .2787596 


The functions C,,(x) and S,(x) defined by equations (6) and (8) have the proper- 


ties 
+'/2 +'/2 
J C..(z) C,(z) dx = f S,.(2) S,(2) dz = bax (10) 


—t/, —I/, 


and 


+1/2 
f Cn(x) Sp(x) dx = 0. (11) 


They form, therefore, a complete set of orthonormal functions in the interval 
(—*/s, +1/2). Thus, if F(x) satisfies the boundary conditions (2), then it can be 
expanded in the form 


z; [A,C,(x) + B,S,(2)], 


n=1 


+'/2 +'/2 
An F(x) C,(x) dx and Ros { F(x) S,(x) dx. 


=x» i 


ae i / 


If F(x) is an even function, then B, = 0, corresponding to the familiar Fourier co- 
sine expansion; conversely, if (x) is odd, then A, = 0, corresponding to the Fourier 
sine expansion. 

3. An Application to Hydrodynamic Stability.—To illustrate the use of these 
functions in problems of hydrodynamic stability, consider the classical Bénard 
problem of thermal instability in a layer of fluid heated from below, for which the 
disturbance equation is 


(D? — a*)*» = —a*Rv, (13) 


where D = d/dx and R is the Rayleigh number. - For two rigid bounding surfaces 
the boundary conditions which » must satisfy are 


= (D? — a*)v = D(D? — a’)v = 0 at xz = oe (14) 


~ 
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Now Chandrasekhar? has shown that two variational principles exist for this prob- 
lem, one of which can be described briefly as follows. 
By letting 
(D? — a*)v = G, (15) 
we can rewrite equation (13) in the form 
(D? — a’)°G = —a?Rr, (16) 


while the boundary conditions (14) require that 


3 ] 
v=G= DG =0 ata. = +. (17) 


The characteristic value, R, can then be expressed as the ratio of two positive def- 


inite integrals in the form 
+1/s 
f [((D? — a*)G}? dx 


l/s 


Dati +1/2 ; Peet i 
of [(Dv)? + a*v?] dx 


Thus, if one assumes a form for G which satisfies the boundary conditions G = DG = 
0 at x = +'/2 and solves equation (15) subject to the boundary conditions v = 0 
at z = +'/., then all the boundary conditions of the problem are satisfied and the 
minimum value of R (with respect to a) can then be obtained from equation (18). 
Now it is known that instability first arises in this problem with respect to a sym- 
metric disturbance. Accordingly, to test the method, we assume the very simple 
trial function G = C(x), with no variational parameters, for which the solution of 


(18) 


equation (15) is 
2d? cosh ax 1 cosh \\x l COS Ayw 


—- - . (19 
Ait — atcosh (2/2) =A? — a? cosh (A;/2) Ay? + a? cos (A, /2) Sk 


{quation (18) then leads to the result 
R = I,/a*he, (20) 
where I, = \i* + 2A,2a? [tanh? (A,/2) — (2/A,) tanh (A,/2)] + a‘ (20’) 


and 


8r,4 Ay 
L= = d, tanh (A,/2) — a tanh (a/2)]+ 
Ot — a’)? (A; tanh (A;/2) (a/2)] 
9 a? 


A 5 , 
ie [tanh? (A,/2) — (2/d,) tanh (A,/2)] — = (20”) 


Ai! ,* — a” 
which gives R = 1,719.1 fora = 3.10. The other variational method,? which re- 
quires the solution of a fourth-order equation, gives R = 1,715.1 for a = 3.12. 
Compared with the exact value, R = 1,707.8 for a = 3.12, both methods give very 
good results. The particular advantage of the present method lies in its greater 
simplicity with little sacrifice in accuracy. 

Expansion in terms of the functions C, and S, can also be used for the solution of 
the problem of the stability of viscous flow between rotating cylinders (cf. Chan- 
drasekhar*). Thus, instead of transforming equation (10) (in the paper just re- 
ferred to) to the form (17) and expanding W in a sine series, we can seek a solution 
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of equation (10) directly by expanding (D? — a?)v in C, and S,, solving for v, and 
then proceeding to the solution of the characteristic-value problem; for many 
problems, a direct expansion in terms of C,, and S, is preferable to one which requires 
a solution of a higher-order equation. 
4. Fourier-Bessel-Type Expansions.—For problems involving circular bounda- 
ries one must consider the solutions of the equation 
a’ ld y?\2 
(z + - =e *) y = aty (21) 
era ° + 
which satisfy the boundary conditions 
dy 
y =0 and = = at r = 1 and n (<\l). (22) 
. dr 


To show that these solutions form an orthogonal set, we proceed as follows. Letting 


d? ld yp? 
p= — -——— 23) 
dr? r dr r2’ ( 
we can show that 
1 1 
f y Of rdr = f f Dy rdr (24) 
i] uv] 


for an arbitrary f, if y and dy/dr vanish at the end points. Thus, from equation 
(21), we have 


1 1 1 
Qm! f Yn, f ar = f 4, Dyj,ra = f DYn DYm 1 dr. (25) 
” " n 


By subtracting from this equation the corresponding equation with m and n inter- 
changed, we obtain 


1 
kin, ass f iy, dr = 0; (26) 
n 
and hence 


1 
f Ymynt dr = 0 ifm # n. (27) 
” 


The solutions y, form, therefore, a complete orthogonal set. We shall proceed to 
obtain the associated characteristic equation and the explicit form of the eigen- 
functions. 
The general solution of equation (21) is 
y = AJ,(ar) + BY,(ar) + CI,(ar) + DK,(ar), (28) 


where A, B, C, and D are arbitrary constants and the definitions of the Bessel func- 
tions are those of Watson.‘ In writing the solution in the foregoing form, we have 
assumed that » is an integer; if vy should not be an integer, we need only write J_, in 
place of Y, to obtain the corresponding solution. 

For purposes of applying the boundary conditions, it is convenient to note that 





7 (r’y) = AJ,-1 (ar) + BY,-1(ar) + CI,-,(ar) — DK,-,(ar). (29) 
a 


The boundary conditions (18) then yield the equations 
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AJ (a) + BY (a) + Cl,(a) + DK,(a) = 0, 
AJ,(an) + BY,(an) + CI,(an) + DK,(an) = 0, 


Ad, -se):-+ BE ox s(a). + Cl ,caledors Dpcake) > @ on 
AJ, _ (an) + BY, _,(an) + CT, —1(an) 8g DK, (an) 7” 0, 
so that the characteristic values a, are the roots of the equation 
| J,(a) Y,(a) I,(a) K,(a@) 
} J,(an) Y,(an) I,(an) K,(an) ee 0 (31) 
J,~1(a) Y,~1(a) T,-1(a) — K,_,(a@) 
J,—1(an) Y,—1(an) I,-1(an) — K,_;(am) 


The corresponding eigenfunctions then become 
Yn = J (anr) + BrY(anr) + Cyl (arr) + D,K,(anr), (32) 
where we have taken A = 1 for convenience. 

The first few roots of equation (31) for »y = 1 are given in Table 2 for three values 
of n, together with the corresponding values of B,, C,, and D,, and of the normal- 
izing integral, V,, the explicit form of which we shall now derive. For this latter 
purpose, consider the quantities 

u(r) = AJ,(ar) + BY, (ar) (33) 
and 
v(r) Cl (ar) + DK,(ar), (34) 


which satisfy the equations 


ld 1 2 Ss 
a (rat) + (a - S)u=0 (35) 
rdr\ ar r? 


and 
ld dv y? p 
~—(r—) —(at+=)v=0. (36) 
rdr\ dr r* 
TABLE 2 
Tue Roots or EquaTION (31) FoR vy = 1 AND RELATED QUANTITIES 
n= 1/4 n = '/2 n= 4/5 
ay +6.39034 +9.49896 +1.89344 X 10! 
Ni +1.34093 +1.42107 +4.5610 xX 107} 
B, —5.87014 —9.07642 —1.02780 X 10} 
C —1.49706 X 107? +1.00635 x 10-8 —8.91136 X 1078 
D, —1.07799 X 10° +4.39637 X 10? —6.65593 X 10° 
a +1.05517 X 10! +1.57397 x 10! 
N2 +2.68344 X 107? +3 .0630 
B, —4.28052 X 107! —1.73711 X 10! 
C2 +4.16730 X 107 +3.67889 K 10~° 
D, — 5.98695 —1.95527 x 104 
a3 +1.47290 X 10? +2.20178 X 10! 
N; +2.45037 X 107? +3.65915 
B; +7.12511 X 107 — 2.24667 X 10! 
C; —7.14451 X 1077 +8.86448 X 107° 
D, —1.99016 X 10! +5.91172 X 105 


Consider now the integral 
1 


1 1 
f r(u + v)? dr = f r(u? + v’) dr +2 f ruv dr. (37) 


” 
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Multiplying equation (35) by rv and equation (36) by ru and subtracting the results, 
we obtain 

1 1 

at f ruv dr = { E : (: 4 —v 4 (- a] dr 
” n dr \ dr dr \ dr 
| ( dv all 
= r u a i ye . 
dr dr] J, 


Since u(r) + v(r) and u’(r) + v’(r) vanish at the end points, the right-hand side 
of this equation vanishes also; equation (37) therefore reduces to 


(38) 


1 1 
J r(u + v)? dr = f r(u? + v?) dr. (39) 


To evaluate the integrals on the right-hand side of this equation, consider first the 
differential equation satisfied by v. If this equation is multiplied by 2rv’(r), then 
one obtains a result that can be written in the form 


d ( ~ ao get Se (40) 


dr dr dr dr 


and when this result is integrated between 7 and 1, we have 


1 dv\? : 
2at [ rv’ dr = | -" ( ) + (ar? + v?)v?] . 
” dr i 


Similarly, 


9 9 $ du : 98 9 ; 
2a? ru? dr = jr* + (a*r? — v?)u? |; 
n dr 


hence 
' f(du\? dv\*\ 1 
2at [ r(u? + v*) dr = E (=) ~ (*) ' + a*r?(u? + v?) + (Ww? — | ‘ 
(43) 


Since u = —vand u’ = —v’ at the end points, this reduces to the very simple result 


1 
i r(u2 + v?) dr = [r’u?]} = [r’v?]) (44) 
= u?(1) — n?u?(n). 


Thus the normalizing integral for y, can be written in the form 


1 
4 ry,? dr = N,, (45) 
n 


where 
N, = [J (an) + BrY (an) |? wns n° [J,( ann) + B,Y (ann) |’, (46) 
the numerical values of which are given in Table 2 for » = 1. 
A function F(r) which satisfies the boundary conditions (22) can therefore be 
expanded in the form 


wo 1 
F(r) = > An,y,(r), where Aj = 1 f F(r) y,(r) rdr, (47) 
n=1 Vn ” 
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and y,(7) are the eigenfunctions defined by equation (32) in close analogy with the 
usual Fourier-Bessel expansions. 

5. Concluding Remarks.—As we have indicated in the Introduction, expansions 
of the type considered here should prove useful in the solution of many hydrody- 
namic and hydromagnetic problems. In particular, the classical problems of the 
stability of viscous flow between rotating cylinders® and in a curved channel,® which 
have never been solved when the spacing between the cylinders is comparable with 
their radii, can be treated by this method. 
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ELECTROMAGNETIC DETERMINATION OF REGIONAL BLOOD FLOW 
IN UNANESTHETIZED ANIMALS* 


By ALEXANDER Ko.uin, NicHoLas ASSALI, GEORGE HERROLD, AND RODERICK 
JENSEN 


UNIVERSITY OF CALIFORNIA MEDICAL SCHOOL, LOS ANGELES 24, CALIFORNIA 
Communicated by D. B. Lindsley, March 26, 1957 


The interest in a method for continuous registration of blood flow in essentially 
normal conscious animals goes beyond the objective of studying flow as a circulatory 
parameter. By far the greatest utility of such a method would be in permitting 
the study of the activity of a great variety of organs by one single method. The 
rate of blood supply to an organ is an important parameter, revealing its rate of 
activity, which determines the required rate of supply of oxygen, hormones, and 
nutrient materials as well as the rate of removal of products of metabolism. The 
blood flow through an organ, when correlated with its state of activity, could be used 
to follow the variations in the activity of the organ. Thus it may become possible, 
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in studies aiming at initial detection of the effectiveness of various stimulating and 
inhibiting factors affecting various organs, to replace discontinuous and often 
cumbersome chemical tests of secretory or excretory functions of these organs, used 
as an index of their activity, by a continuous record of blood flow. The applicability 
of such studies to the solution of pharmacological and physiological problems as well 
as to investigations in experimental pathology is obvious. ; 

An ideal method for the determination of regional blood flow in normal animals 
should have the following characteristics: it should measure flow through intact 
blood vessels; the blood vessel should be left in its normal place and not be exposed ; 
the measuring device and procedure should be such that drugs (sth as heparin or 
anesthetics) would not be needed, the animal being as nearly as possible in its normal 
state; the inserted instruments should be suitable for chronic implantation; the 
relation between blood flow and instrument reading should be linear and independent 
of variations in phyiscal factors like temperature, viscosity, velocity profile (laminar 
or turbulent flow), etc.; the response of the instrument to velocity changes should be 
instantaneous; and frequent calibrations of the instrument should not be necessary. 
It is believed that the method to be described below is the nearest approach so far 
to these desired ideal conditions. 

1. Principle and Development of the Method.—The principle of the method is 
based on the phenomenon of electromagnetic induction discovered by Faraday. 
Faraday demonstrated the generation of electromotive forces (emf’s) in closed 
electric circuits traversed by a varying magnetic flux, as well as in solid conductors 
moving at right angles to a magnetic field (Faraday’s unipolar inductor). Faraday 
could easily have demonstrated, by means available in his laboratory, the induction 
of an emf in a liquid traversing a magnetic field, by allowing mercury to flow through 
a pipe disposed at right angles to the magnetic field and by picking up the voltage 
induced in the moving liquid by means of two electrodes disposed at the ends of the 
pipe diameter perpendicular to the magnetic lines of force (Fig. 1A). Strangely 
enough, Faraday did not perform such an experiment but, instead, spent a consider- 
able amount of time conducting experiments at the Waterloo Bridge in London, 
trying to demonstrate the induction of an electromotive force in the River Thames 
as its water moves across the magnetic field of the earth. Owing to difficulties caused 
by polarization of the electrodes, Faraday finally gave up these experiments. This 
negative result is mentioned by him in a Bakerian lecture.! 

The qualitative idea of the possibility of inducing an emf in a moving fluid, which 
Faraday evidently conceived, is, owing to complications introduced by the velocity 
distribution in the flowing fluid, still considerably removed from the realization of the 
possibility of quantitative measurement of flow by electromagnetic induction. The 
importance of this factor was realized by Williams,” who, in fact, used electromag- 
netic induction for an indirect determination of the velocity distribution in pipes, 
using a special nonpolarizable system (copper electrodes in a copper sulfate solution). 
Although Williams’ results imply the possibility of development of a flowmeter 
based on electromagnetic induction, there is no such suggestion made in his paper. 
That paper was not followed by further publications on this subject. 

Unaware of the work of Williams and Faraday, Fabre proposed in a letter to 
the French Academy of Sciences’ a qualitative biological use (“Cun hémodromograph 
sans palette’’) for electromagnetic induction. His note suggested the possibility of 
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Fig. 1A4.—Scheme of electromagnetic flowmeter. N and S, pole pieces of magnet; + and 
—, electrodes which pass through the pipe wall to make contact with the fluid. 

Fig. 1B.—Scheme for determination of blood flow in intact artery. H, magnetic field vector; 
A, artery; EH, and EF», nonpolarizable electrodes touching the artery wall. 

Fig. 1C.—Pickup sleeve S with artery and nonpolarizable electrodes. W, wicks soaked in agar- 
saline; B, porous clay boots filled with zine sulfate; Z, amalgamated Zn electrodes. 

Fig. 1D.—Pickup sleeve with metal electrodes. Si, slit through which artery is slipped into 
the sleeve channel; E; and £2, pickup electrodes contacting the artery; W; and W2, electrode lead 
wires. 


detecting variations in flow of blood through a pipe equipped with electrodes and 
tied between the two parts of a cut blood vessel. His discussion was qualitative, 
and no experimental evidence was presented. The instrumentation described by 
him could not have permitted measurement of the rate of flow, which evidently 
was not attempted. This note was not followed by further publications on this 
subject. 

The above publications have had no effect om the development of the electro- 
magnetic flowmeter. The method was developed independently by investigators 
who have been unaware of the close approaches of the preceding authors. In 1936 
it was demonstrated‘ that electromagnetic induction could be used for quantitative 
determination of fluid flow in pipes, and the anticipated linear relation between the 
volumetric rate of flow and induced voltage was demonstrated. Furthermore, it 
was shown that flow of blood through arteries could be recorded without injuring 
the artery wall, merely by applying nonpolarizable electrodes to the outside of the 
walls of the artery at the ends of its diameter perpendicular to the magnetic field 
(H) (Fig. 1B). At the same time, the instantaneous response of the apparatus to 
rapidly pulsating flow was demonstrated, and the average rate of flow of blood 
through a dog’s carotid artery was determined. The instrumentation used in this 
communication was very simple, consisting of a d’Arsonval galvanometer for 
average-flow measurement and a string galvanometer for recording of instantaneous 
flow. In the subsequent communications’? the use of a chopper amplifier was 
described, and a sleeve was introduced for convenient application of nonpolarizable 
electrodes (Fig. 1C). 
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Shortly after Kolin’s 1936 publication,‘ an independent paper by Wetterer 
appeared,® describing the same method in a form very similar to that presented 
by Kolin. Experimental evidence in the form of calibration curves and blood-flow 
records was presented. The use of a d.c. amplifier was described in a subsequent 
paper.’ 

The initial version of the method, as described by Kolin* * 7 and Katz and Kolin,® 
utilizes a constant magnetic field, so that a d.c. emf is induced with unidirectional 
flow. This necessitates the use of nonpolarizable electrodes and of a suitable 
system of d.c. amplification. The nonpolarizable electrodes are quite cumber- 
some in use,’ and d.c. amplifiers are, as a rule, less stable than a.c. amplifiers. The 
next logical step thus consisted of replacing the constant magnetic field by an alter- 
nating one.'! In this case, an a.c. emf is induced in the case of unidirectional flow, 
so that nonpolarizable electrodes become unnecessary, ordinary metal electrodes 
touching the artery wall being quite adequate (Fig. 1D). An a.c. amplifier can be 
used to detect the signal induced due to flow of blood. It is, however, necessary to 
use a special compensator,'? in order to adjust for a zero reading at zero flow. Such 
zero adjustment requires temporary occlusion of the blood vessel by special occlusion 
clamps or pneumatic devices.'! '% 14 

In the initial experiments*~'! the blood vessel, 
usually an artery, had to remain exposed during the 
measurement of blood flow. Thus anesthesia had to 
be maintained. A step toward elimination of an- 
esthesia was taken by implanting the lucite sleeves 
containing the pickup electrodes into the animal, 
leaving the blood vessel in its normal place.'*!* In 
addition, a way was found to establish a zero base 
line without occluding the blood vessel to stop the 
blood flow.'* This was accomplished by means of 
pickup sleeves with four electrodes. After the im- 

z plantation of the pickup unit, the animal was placed 

Fic. 2.—External magnet 7 ’ 
used in conjunction with im- in the gap of a large electromagnet (Fig. 2), and the 


bedded sleeve of type shown flow was recorded as in previous modifications of the 
in Fig. 1D. P, and P2, mag- hod.! 16 

net pole pieces; C, magnet method." 

coils; S, grounded shield. The drawback of the above modification of the 


ie 23 “ae nes st method lies in the necessity of maintaining the ani- 

mal immobile in the magnetic field to avoid unbalanc- 
ing the zero adjustment. The possibility of such a disturbance is avoided in the 
modification presented in this paper. This is achieved by a rigid coupling between 
the magnet, the electrodes, and the blood vessel. This coupling is accomplished by 
using a miniature magnet unit (Figs. 3A [front view] and 3B [side view)) with a 
pickup sleeve (S) fixed in its gap. This entire unit is implanted into the animal, 
the blood vessel under study being inserted into the pickup sleeve channel through 
the slit Sl (Fig. 3C). 

2. Design of the Miniature Flowmeter Implant.—Figure 3C shows a photograph 
of a magnet used for flow measurement in the carotid artery of a dog. The lucite 
sleeve S is inserted in the gap of the magnet. It was found most practical to use 
miniature magnets manufactured for erase heads and recording heads of magnetic 
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Fic. 3B.—Side view of magnet M and sleeve S. 
W, suspension wires; L, electrode lead wires 
Fig. 3A.—Front view of magnet scale in inches). 
£ 


tape recorders. The magnet shown is from a Brush Electronics Company magnetic- 


tape erase head, Model BK944. Its dimensions are indicated by the inch scale. 
Its weight was 5.5gm. The coils had to be rewound on a nylon or teflon bobbin, since 
the commercially provided one dissolves in the monomer used in the polymerization 
process (see below). The total number of turns of both coils (439 formvar insulated 
copper magnet wire) was 1,600, and the current used in continuous operation was 
65 milliamperes, which resulted in a magnetic field of 800 oersteds in the gap. The 
electrodes were made of */,-mm.-diameter brass rod. The tips of the brass elec- 
trodes, flush with the inner surface of the sleeve S, were gold-plated after the com- 
pletion of the flowmeter unit. Figure 4 shows the sleeve S in the gap of the magnet 
M (greatly magnified). P; and P, are the magnet pole pieces, H, and FZ; are the 
electrodes, S/ is the slit, and W, and W, are the lead wires connecting the electrodes 
to the amplifier. _W, has been shown for clarity inside the sleeve body, but actually 
it runs through a shallow groove in the lucite which is located between the sleeve 
and the pole piece P; coplanar with the line joining £, to Z2. The two wires W, and 
W, are twisted to minimize the emf induced in the input leads. 

In order to secure adequate electrical insulation between the magnet wire, the 
electrode leads, and the animal, the combination of magnet and sleeve is cast into 
a block of plastic (methyl methacrylate). The limits of the plastic are indicated by 
dashed lines in Figure 4. The magnet lead wires and the electrode leads are in- 
sulated by separate flexible polyvinyl tubes (°/3:-inch i.d.; '/s-inch o.d.) through 
which they run (7; and 7:). A ground wire connected to the laminated magnet 
core made of sheets of silicon steel (*/; mm. thick) is imbedded in the same tubing as 
the magnet leads. The polyvinyl tubing is polymerized within the methyl meth- 
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Fic. 3C.—View of flowmeter encapsulated in pastic. SI, 
slit; C, artery channel. 


acrylate plastic, so that a very secure watertight seal is achieved. The adequate 
moistureproof electrical insulation between the above-mentioned components of the 
system was the main experimental difficulty which had to be overcome. Miniature 
magnets, small enough to be inserted into the neck of a cat, were prepared by one of 
the authors (A. K.) in 1952 (unpublished), and successful tracings of blood flow 
in the cat’s carotid artery have been obtained. But the flowmeter units were of 
limited value, owing to an inadequate seal between the polyethylene tubing which 
was used to insulate the leads and the insulating material encapsulating the magnet. 
The present method of insulation results in flowmeter units which remain in perfect 
condition even after prolonged implantation into an animal. 

Figure 5 illustrates the casting procedure. The magnet M (seen in side view) is 
suspended with the gap (space behind P;) facing downward in a glass vial V. The 
lucite sleeve S, with the pickup electrode EZ, showing at the bottom, is inserted in the 
gap. The polyvinyl tubing (7 and 72) insulating the lead wires dips into the vial. 
The methy! methacrylate is then poured into the vial. The vial is now immersed 
for 6-8 hours in a constant-temperature bath kept at 120° F. One milligram of 
benzol peroxide is then added per cubic centimeter of monomer, to act as an accelera- 
tor for the polymerization process. The casting is accomplished by polymerizing 
several layers of the monomer successively. Upon completion of the polymeriza- 
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CROSS SECTION OF FINISHED FLOW METER 


F ia. 4.—Front view of electromagnetic flowmeter. 7, and 7, 
insulating polyvinyl tubing; M, magnet core; C; and C2, magnet 
coils; W; and W2, electrode lead wires; FE, and E2, pickup elec- 
aeomts P; and P», pole pieces of magnet; S, pickup sleeve; Si, 
slit. 


tion process, the glass vial is broken and excess lucite (polymethyl methacrylate) is 
removed by abrasive action of a wet sanding-machine belt. 

Prior to the polymerization, the cylindrical channel of the sleeve S (Figs. 4 and 5) 
is filled with a low-melting-point metal alloy (Cerrolow 136). This secures the align- 
ment of the electrode pin (a */,-mm. brass rod of which F, and £; are initially the 
ends) during the process of polymerization, during which the lucite sleeve, which 
supports the electrodes, may dissolve. After solidification of the plastic, the 
Cerrolow core is drilled out, and with it the portion of the brass rod which lies be- 
tween FL, and £,. The surfaces of the electrodes EZ; and E; thus produced are flush 
with the inside wall of the sleeve. The last step consists in gold-plating the elec- 
trode tips electrolytically. After removal of all sharp edges and polishing, the 
flowmeter unit is ready for implantation. A more detailed description of the 
manufacturing process of the flowmeter units will be published elsewhere. 

3. Calibration and Response to Pulsating Flow.—The emf (in volts) induced in a 
conducting liquid flowing through a pipe at right angles to a magnetic field is given 
by the expression 


V = 10° pds, (1) 


where yu is the permeability, H is the magnetic field strength in oersteds, d is the 
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be 21", 
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FLOW METER SUSPENDED IN VIAL Salt Concentration in ¥% 


Fic, 6.—Independence of the flowmeter sensi- 
tivity of the electrical conductivity of the solu- 
tion. A constant flow is maintained through an 
a.c. pipe flowmeter. The conductivity of the 
wires; W, end Wes, seaeaalin liquid is varied by varying the concentration of 
wires: TP, top plate (aluminum); NaCl solution (abscissa). The instrument read- 

dave ? ing is measured by the ordinate. (From A. 


M, suspended magnet; C), one Oy. ney we OQ 104% 
of the magnet calls: ¥, ines Kolin, Rev. Sci. Instr., 16, 109, 1945.) 


vial; P;, one of the magnet pole 
pieces; S, pickup sleeve; F,, 
protruding pickup _ electrode; 
BP, base plate (aluminum); SP; 
and SP2, support pins. 


Fia. 5.—Magnet suspended in 
vial for encapsulation in plastic 
(side view). F; and F2, clamps 
for fastening the suspension 


diameter of the pipe in centimeters, and v is the average velocity of the liquid in 
centimeters per second.‘ This emf is independent of the electrical conductivity of 
the liquid (Fig. 6, reproduced from an earlier paper") as well as of the velocity 
distribution, provided that the latter is axially symmetrical.'© Contrary to intuitive 
expectation, the emf induced in the moving fluid for a given volume rate of flow is the 
same regardless of whether the fluid flows through a nonconducting pipe or 
through an artery whose electrical conductivity is nearly the same as that of the 
fluid. This is true regardless of the thickness of the artery wall.'‘ This paradox 
is a consequence of the independence of the induced emf of the velocity distribution, 
as pointed out in an earlier paper.'‘ It is thus unnecessary to calibrate the instru- 
ment with each individual artery. One single calibration, once performed without 
an artery, standardizes the instrument. These results do not depend on the character 
of the magnetic field used. They are valid for constant as well as for alternating 
fields; for sinusoidal as well as for nonsinusoidal ones, such as the square wave 
shape used by Dennison, Spencer, and Green.!” 

The linearity of the calibration is apparent from equation (1) and is verified by 
experimental results using an excised carotid artery (Fig. 7A, reproduced from a 
previous paper‘). The same result is obtained when calibration is carried out with 
a pulsating flow in vivo (Fig. 7B). The instantaneity of the response of the flow- 
meter to variations in flow has been confirmed by Inouye and Kuga" through care- 
ful experimental studies. 

It is sufficient to establish the calibration for one single size of pipe. The instru- 
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ment sensitivity for other diameters follows easily from the following consideration: 
The emf induced for the rate of discharge Q = rR*# is given by the equation 
2R-Q  10-*2uHQ 


V = 10° wH(2R) - 6 = 10 pH — 


Stee Fe 2 
wh? rR 2) 


V is thus inversely proportional to R at a constant rate of discharge Q. The sensi- 
tivity S = V/Q can thus be expressed as follows: 
ee 7m 
iy = (9) 
R 
where K = 10~* (2uH/7) is an instrument constant. Having determined S for one 
particular diameter, we have calibrated the instrument for all artery diameters. 
In the case of the miniature-magnet method, one has the choice of perfusing the 
sleeve channel with saline at a known rate prior to cutting the slit S/ or of inserting 
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‘ F ‘ Fi — Fic. 7B.—Calibration is carried out with a pul- 

\ 2 3 4 5 sating blood flow through a dog’s carotid artery 

RATE OF FLOW IN CC/sEC in vwo. A constant-field flowmeter is used. 
: ‘ (From A. Kolin, Am. J. Phystol., 122, 797, 

Fic. 7A.—Calibration curve for a.c. 1938. ) 

blood flowmeter. An excised carotid 

artery is perfused with blood at con- 

stant rates. (From A. Kolin, Rev. 

Sci. Instr., 23, 235, 1952.) 





an excised artery through the slit and perfusing it with blood. Practically, the 
perfusion consists of ejection of a known volume of saline or blood through the flow- 
meter unit from a syringe while recording the flow signal. The record thus obtained 
is shown in Figure 10, B,c. The injection time is given by the extent of the record 
along the horizontal time axis. The ratio of the injected volume to this time is the 
average flow. The average-flow signal is obtained by planimetry of the record of the 
instantaneous flow and division of the area obtained by the horizontal extent of 
the curve. The ratio of the average-flow signal thus obtained to the average flow 
gives the instrument sensitivity in terms of the instrument reading in centimeters per 
unit flow of 1 ce/sec. In the experiments described below, the sensitivity was 4 mm. 
per ec/sec, but higher values can easily be obtained by increasing the amplification. 

4. Apparatus.—The amplifying and recording apparatus used was essentially 
the same as described in a previous paper.‘ The magnet was energized by a 60- 
cycle alternating current. Figure 8 shows the block diagram of the apparatus. The 
signal voltage derived from the flowmeter FM is put in series with the output of the 
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compensator (¢omp.). The compensator voltage can be varied independently in 
phase and in amplitude,* by one single control knob. This permits rapid com- 
pensation of stray 60-cycle pickup in the input circuit which may mask the flow 
signal to be recorded. This compensation must be accomplished at zero flow. In 
order to be able to stop the flow through the artery, a pneumatic device was used. It 
is an improvement upon the occlusion bags described earlier.‘ Figure 9 shows the 
occlusion unit and its components. A tuning-fork-shaped form figure (a) is pre- 
pared from a low-melting-point alloy (Wood’s metal). By repeated dipping of this 
form in a solution of polyvinyl in a suitable solvent (“Cadco 201”’), a sufficiently 
strong skin of polyvinyl sheeting is formed on the surface of the metal after the 
evaporation of the solvent. The Y-shaped form is then dipped into water at 100° 
C. The metal form melts and flows out of the polyvinyl bag. The open 
central ‘‘tail” of the Y is cemented by polyvinyl cement (solution of polyvinyl in 
“Cadeo 201”) to a long polyvinyl 
tubing which is attached to a 
. syringe, while the two remaining 
@ closed legs of the Y are imbedded in 
the animal, being placed about the 
artery as shown in Figure 9, 6. A 
hull of tygon tubing, shown in Fig- 
Yr ti ure 9, c, is then pushed over the in- 
sert as shown in Figure 9, d, to pre- 
vent outward expansion of the bags 
~ aaa and to insure perfect occlusion of 
the vessel when the syringe plunger 
_ Fie. 8.—Block diagram of the flowmeter circuit. is pushed in. The tygon hull is 
FM, the flowmeter transducer comp. COMEDY finally closed by stitches near the 
end of the split section. Figure 10 
illustrates the effectiveness of this occlusion device. Short base lines, like those 
shown in Figure 10, are obtained frequently throughout the measurements, so that a 
continuous base line can be obtained by joining them. It is during such occlusion 
periods that the base line is occasionally rapidly adjusted to zero. Once the ad- 
justment is made, readjustment is rarely required. 

The combined signals of the flowmeter and of the compensator are conveyed in 
series to the amplifier (RC-coupled single-ended 3-stage) input. The output signal 
of the amplifier is applied to the oscilloscope (osc.) (Dumont type 304), and the 
oscilloscope reading is recorded by a continuous automatically developing camera, 
as described in a previous paper. '4 

5. Examples of Determination of Blood Flow in Unanesthetized Animals.—Figures 
10-12 show records of blood flow obtained with an implanted miniature flowmeter 
unit of the type described above. These records are not primarily intended to 
demonstrate physiological results but, rather, are presented as an illustration, ob- 
tained with one particular animal, of the possibilities of application of the new 
modification of the electromagnetic blood flowmeter. 

The experiments presented below were obtained with a mongrel dog (weighing 
14 kg.) two days after the implantation of the measuring unit. The dog was un- 
anesthetized and was strapped in standing position to avoid moving too far from the 
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Fia. 9.—Components of the occlusion bag in relation to the artery 
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Fia. 10.—Effect of ether anesthesia upon blood flow in dog’s left carotid 
artery and effect of inhalation of smoke. A, onset of ether anesthesia (mask 
applied at a). B, recovery from ether anesthesia (mask removed at a; dog 
regains consciousness at b; c, calibration record). C, effect of inhalation of 
smoke upon carotid blood flow; inhalation begins at a. The base line is ob- 
tained by pneumatic occlusion of the artery. 
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Fig. 11.—Effects of a visual stimulus and of a pharmacological agent upon 

blood flow through left carotid artery in dog. A, reaction to flame (a torch is 

presented to the dog at the interruption of the record). B, effect of amyl 

nitrite (inhalation begins ata). C, recovery from effect of amyl! nitrite. The 
inhalation ends at a. 


‘ehooogggrgtiabe 


4 


? : (NAMA WY 


the tageyberagg tt! 


AAAANOANAAAANOANAAY! 


ee Mat ee ee ee Mee ae rte Mae 


° 
=< 





py Petee gette 


4 prota; cee asghesen tc PATIL eee eee Ee 4 
B 5 err Tae RON OTIC 
° emer 


i | 





FLOW IN CC/SEC.® 


TIME IN SEC. 


Fia. 12.—A, effect of occlusion of carotid artery interrupting the flow 
signal upon the subsequent carotid blood flow. Record taken in left carotid 
artery. The other carotid was tied. B, reaction to food. Blood flow in left 
carotid artery. Food is presented at a. Chewing and swallowing of the food 
begin at b. 


recording apparatus. Blood flow was recorded in all the subsequent tracings in the 
left carotid artery; the right carotid artery was ligated. In Figure 10, A and B, 
the profound effect of anesthesia upon the blood flow in the carotid artery is illus- 
trated. The first section of Figure 10, A, represents the normal control tracing. It 
is followed by a brief occlusion of the carotid artery to establish the base line at zero 
flow. The ether mask is then applied to the dog at a, which results in a marked 
reduction in blood flow. Figure 10, B, shows the increase in carotid blood flow dur- 
ing recovery from ether anesthesia. The mask isremoved ata. The sharp increase 
in flow upon regaining consciousness at b is noteworthy (c is a calibration record). 
Figure 10, C, shows a marked effect of inhalation of cigarette smoke (inhaled at a) 
upon the carotid blood flow. 
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In Figure 11, A, the effect of excitement by presentation of a flame (torch) upon 
the carotid blood flow is illustrated. In Figure 11, B and C, a pharmacological 
effect is depicted. At B, a, the animal begins inhaling amy] nitrite. The reversal 
of the normal direction of diastolic flow during the initial diastoles is noteworthy. 
At C, a, the inhalation of amy] nitrite is discontinued and the blood flow goes up. 
In Figure 12 the filmstrip A shows the effect of a brief occlusion of the carotid artery 
upon the subsequent blood flow through it. As can be seen from this record, the 
initial brief increase in flow following release of the artery is substantial (approxi- 
mately 45 per cent). Figure 12, B. shows the effect of stimulation by food upon 
carotid blood flow. At B, a, food is presented to the dog. There is an immediate 
noticeable increase in flow; at b, as the animal begins to chew the food (ground 
meat), the flow increases and becomes somewhat less regular. 

Other types of stimuli have been tried, such as frightening the animal by a light 
flash. It was interesting to note the very rapid conditioning which occurred. The 
carotid blood flow was reduced upon the presentation of the light bulb to the animal 
prior to flashing it for the second time. 

These examples illustrate the possibility of studying the effects of physiological 
and psychological stimuli as well as of pharmacological agents upon essentially 
normal conscious animals. Chronic research in experimental pathology also appears 
to be within the realm of possibility. Until now, the implants have not been left in 
the animal for more than a week, because of thrombosis caused probably by excessive 
constriction of the artery. There are good reasons to expect that, after elimination 
of this condition, the implants could be left in the animal indefinitely. The demon- 
stration of the absence of objectionable biological effects of the implants, when prop- 
erly used, as well as of physical artifacts affecting the measurement of flow, will be 
presented in a subsequent publication. 

The preceding examples demonstrate some of the main advantages of the new 
modification of the electromagnetic method over the traditional methods of blood- 
flow determination. The possibility of obtaining a continuous flow record showing 
instantaneous responses without interference with the flow enables one to detect 
phenomena which are not observable by means of conventional methods. The 
possibility of making chronic implantations of sensing elements suggests the applica- 
bility of this method to research in experimental pathology and therapeutics. The 
use of conscious animals promises pharmacological results of greater validity than 
those obtained in earlier studies on anesthetized animals and permits the applica- 
tion of this method to problems in physiological psychology. Finally, the possi- 
bility of implanting miniature flowmeters simultaneously about several different 
arteries will make it possible to study the correlation between the functions of 
several different organs and their simultaneous response to a variety of stimuli. 

The next step in this work will consist of improving the apparatus and procedure 
so as to make it adaptable to the measurement of blood flow in human beings and in 
animals in any accessible blood vessel of sufficiently large diameter. It will be pos- 
sible to measure blood flow in arteries '/,mm. in diameter. This is not necessarily 
the lowest attainable limit. 
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Communicated by Horace W. Magoun, April 8, 1957 


Several attempts'~* have been made to study the combination between adenosine 
triphosphate (ATP) and actomyosin quantitatively, by determination of the vis- 
cosimetrically or turbidimetrically measurable change induced in the protein by its 
combination with ATP. The first efforts + * were preliminary, in that the avail- 
able methods did not permit sufficiently rapid measurement of these transient 
phenomena. But even the measurements performed with a rapid technique‘ are 
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not quantitative, since, by the time the response has become maximal, an unevalu- 
ated fraction of the ATP has been hydrolyzed and replaced by ADP, which in- 
hibits competitively.’ Hence all previous estimates of the quantity of actomyosin 
reacting with 1 mole of ATP are too low. 

A new method of determining the stoichiometry of this interaction is provided 
by the experiments of Mommaerts and Hanson,‘ who measured the turbidity 
change Ar caused by small amounts of ATP maintained by a rephosphorylating 
enzyme (pyruvate kinase) and a reservoir substrate (phosphoenol pyruvate), at 
pH 6.9, 0.01 M MgCle, ionic strength about 0.5 (mainly KCl). From the ATP- 
actomyosin association curve so obtained (Mommaerts and Hanson,® Fig. 6), we 
shall derive the total number n of moles of ATP bound by 10® grams of actomyosin 
necessary to produce the total physical change, Armax, and, furthermore, calculate 
the intrinsic association constant, k, of this combination reaction. 

It is assumed that the relative turbidity drop is proportional to the amount 
of bound, ATP, r being the bound ATP in moles per 10° grams of actomyosin: 

: Ar 

Dis ghetaay reais 

n ATmax 
The concentration of bound ATP is r(./)-10~* moles per liter for an actomyosin 
concentration of (M) grams per liter and becomes n(M/)-10~* = x for the case of 
maximal response. If c’ is the concentration of free ATP, and c that of the total 
ATP, we can write’ 


and hence 


r r/n a 


~ (n — roe’ ie 








(1 — r/n)(e — ra /n) Pe (1 — a)(e emis ax) 


This is valid for any combination of points 7, 7, on the association curve, and z is 
obtained from 


ai aj; 


(l—ad(e— aa) (1 — a)(c; — aye)’ 


In this manner we obtain n = 1.98 + 0.12 moles per 10° grams, so that 1 mole of 
ATP is bound per 504,000 grams of actomyosin. This equivalent of protein may 
still be somewhat underestimated, since it has been implied that all the nucleotide 
is maintained in the form of ATP. This condition was aimed at in the experiments 
by the use of a large excess of pyruvate kinase, but it is unknown as yet to what 
extent this was successful. Any imperfection in this regard would tend to reduce c 
below the accepted values and hence to increase the equivalent weight of the pro- 
tein. We decline, for the moment, to interpret the result in terms of the molecu- 
lar weight of myosin, in view of the uncertainty of the latter. 

Further, we find k = (6.7 + 0.15)- 10°, and hence AF° 20 = RT Ink = 9,450 eal. 
The fact that a constant value of k applies throughout the curve corroborates the 
underlying assumption that Ar is proportional to a, unless there chances to be a 
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compensation of errors due to this factor and to differences in the completeness by 
which ¢ is maintained by rephosphorylation at different ATP levels. It is of 
importance to emphasize the high value of AF° which, having been estimated be- 
fore,* * is now obtained with more satisfactory precision and is found to exceed 
even the AF° and AH values of the hydrolysis of ATP.*-" It should be kept in 
mind that the value of k, notwithstanding its seemingly different mode of deter- 
mination, is still essentially identical with a Michaelis-Menten constant. Hence 
it is a true association constant only if the velocity constants for combination and 
dissociation are considerably greater than that for the enzymatic decomposition. 
If this were not the case, k and AF° would have been underestimated. 

Considering the excess of Mg over ATP in these experiments, it should be kept 
in mind that practically all ATP is in the form of the Mg complex! '?, so that the 
above derivation applies to this rather than to ATP as such. 
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Angeles, Calif. 
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RESULTS OF INVERSION OF NEURAL PLATE MATERIAL* 


By J. S. NicHouas 


OSBORN ZOOLOGICAL LABORATORY, YALE UNIVERSITY 
Communicated April 13, 1957 


Capacities for differentiation have frequently been tested by transplanting or ex- 
planting parts of the embryonic nervous system. Experiments of this nature have 
contributed much information concerning the time of tissue determination, of 
regional determination, and specificity.'~* In the present investigation the ability 
of the medullary plate to form nervous system under varied conditions of auto- 
transplantation has been studied. 

It is well known that the chordamesoderm influences the overlying ectoderm to 
form nervous system. The present experiments were made to’study development 
of the nervous system under conditions of inversion. A strip of medullary plate 
together with underlying archenteric roof was cut out of neurulae of the spotted 
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salamander Amblystoma punctatum and replaced in all possible varying axial rela- 
tions. The following questions were asked: (1) How do these reversals affect the 
nervous system and the mesodermal derivatives? (2) Will the structures formed 
adjust themselves to the level of the embryo as a whole? (3) Will the reim- 
planted structures be otherwise adjusted to their surroundings? 

The operations performed are diagramed in Figures 1-3. Axes of the strips 
were reversed as follows: anteroposterior and mesiolateral (Fig. 1); mesiolateral 
and dorsoventral (Fig. 2); and anteroposterior and dorsoventral (Fig. 3). Stages 
used varied from Harrison’s Stage 14 to Stage 16. The level of the strip was varied 
either anteriorly or posteriorly in the embryo to test its regional specificity. 

Healing in all series occurred remarkably rapidly. The graft maintained its 
shape. In the operations involving dorsoventral reversal the graft remained 
sharply outlined because of the contrast between the lightly pigmented chordameso- 
derm and the fully pigmented neural plate surrounding it. 


1 3 


Fias. 1-3.—Diagrams of the early neurula of A. punctatum, showing the field of operation and 
indicating the orientation of the strip of the prospective neural plate which was modified experi- 
mentally in its axial relations. Fig. 1, anteroposterior and mesiolateral axes reversed. Fig. 
mesiolateral and dorsoventral axes reversed. Fig. 3, anteroposterior and dorsoventral axes re- 
versed. 


In the most favorable cases in the first series (as in Fig. 1), reversal without in- 
version, the closing of the neural folds was not retarded. The folds arched over the 
graft and obscured it from view. Subsequent development was fairly normal dur- 
ing the early stages, except for a thickening of the neural tube in the vicinity of the 
region of operation. As development progressed, the larvae were seen to have 
shortened bodies; it was suspected that this was due to defective development of 
the notochord, a suspicion later justified by study of the sections. 

The nervous system in the first series (Fig. 1) has always developed regularly, 
irrespective of the level of operation. Marked deficits, however, have been found 
in the principal axial mesodermal structure, the notochord. In the less favorable 
cases, particularly those performed in the late neurula, the irregularities were chiefly 
those of blockage and malformation of the posterior head and trunk. The chorda- 
mesoderm evidently served as a partial block to one or the other of the enveloping 
neural folds. 

In the second and third series (Figs. 2 and 3), where the dorsoventral axis was 
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reversed, an entirely different course of development resulted. The nervous sys- 
tem from the graft was inverted. The only animals in which the transplanted 
nervous system showed any regularity were among the group operated at Stage 
14. The inverted strip acted as a block to the nerve tube of the host; at the level 
of the graft it can be seen in sections that three nerve tubes are present (Figs. 4 and 
5). Two lateral ones derive from the host medullary folds, a median one from the 
graft. In the case of the series (Fig. 2) in which the anteroposterior axis of the 
transplant was not reversed, the nervous structures induced in the graft proved to 
be of the appropriate level, as would be expected. 

The chordamesoderm in inverted position is in most cases prevented from attaining 
its characteristic organization. It forms dissonant elements of unique cellular 
structure, as shown in Figures 4 and 5. In the third series, where the anteroposte- 
rior axis was reversed, the notochord appeared to form a more extensive indeter- 
minate mass than in the second. 

The chordamesoderm had already gastrulated normally before the operation and 


4, 
en me Pl, : 

Fias. 4, 5.-The photomicrographs show the development of three typical nervous systems re- 
sulting from the operation shown in Fig. 3. The median nervous system has developed from the 
grafted strip; the two lateral ones, from regularly developing neural folds which were separated by 
cordamesoderm during the period of closing of the neural folds. The inchoate mass between the 
three nervous systems is composed of chordamesoderm which forms an indeterminate cellular mass 
lacking characteristic notochordal structure. 


had performed its primary induction. Marx,‘ in Spemann’s® laboratory, had 
shown that, under other conditions of transplantation, the notochord-forming region 
was the one region of the early gastrula able to realize its own potentialities. The 
most striking result of the present experiment has been to show that at a later 
stage, after the neural plate and the archenteric roof have started on the course of 
their differentiation, the notochord is far more dependent on its surroundings than 
is the nervous system. 

Thus, in answer to the questions posed above, the following would appear to be 
true: (1) The nervous system maintains its identity under the various sorts of re- 
versal imposed. The notochord fares less well. Even if left beneath the ectoderm, 
reversal of anteroposterior and transverse axes appears to disturb its elongation and 
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its integrity. If differentiating notochord is placed dorsal to the medullary plate, 
it evidently ceases to develop as a structural entity, although cellular differentiation 
progresses. (2) The nervous system at the time of operation is capable of regulat- 
ing, with regard to level, reversals of the anteroposterior and transverse axes but 
not reversals of the dorsoventral and anteroposterior axes imposed simultaneously. 
The notochord seems not to regulate anteroposterior reversal. (3) The archenteric 
roof and medullary plate retain their respective identities under conditions of re- 
versal and do not adjust to their new position when the dorsoventral axis is re- 
versed. In such a graft the chordamesoderm acts as an intrusive developmental 
block and separates locally the developing neural folds each of which forms a 
complete reduced nervous system characteristic of its original localization. 

* Partially supported by the Anna Fuller Fund and by a grant from the American Cancer 
Society administered through the Yale University Committee on Atypical Growth. This paper 
is dedicated to Professor Alexander Petrunkevitch, my teacher, colleague, and friend, on the 
occasion of his eightieth birthday. 
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